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THE EFFICACY OF FLOOR SILLS UNDER DROWNED HYDRAULIC JUMPS 


Ahmed Shukry,! M. ASCE 
(Proc. Paper 1260) 


SYNOPSIS 


In this paper are presented the results of tests on the hydraulic perform- 
ance of various types of floor sills. The study is mainly concerned with low- 
head river barrages which are generally operated under conditions of 
drowned jumps. Example of this type of barrages is the Edfina Barrage 
which controls the escapage discharges of the Rosetta Branch of the River 
Nile into the sea. A scale model of this barrage was constructed with a fixed 
bed. The distribution of velocities, for various types and locations of sills, 
were recorded by a Pitot-static tube. The efficiency of any sill against bed- 
scour is indicated by the rate of adjustment of the flow to the normal distri- 
bution in the downstream channel. 


INTRODUCTION 


The hydraulic jump which forms within or downstream of an hydraulic te 
structure may be described as a: (a) perfect jump, (b) repelled jump that . 
takes place downstream of the location of a perfect jump, and (c) drowned "3 
jump that occurs upstream of the location of a perfect jump. ee 

On horizontal floors, the type of jump may be predicted from the well- 3 
known formula: 


24 Y 
4 2 


in which Yq is the depth of flow before the jump, YQ is the conjugate depth of 


Note: Discussion open until November 1, 1957. Paper 1260 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. aS 
83, No. HY 3, June, 1957. 3 
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flow after the jump, q is the discharge per unit width and g is the acceleration 
of gravity. 

A perfect jump is formed when the depth of water in the downstream chan- 
nel is equal to the conjugate depth Y2g. On the other hand, the repelled and the 
drowned jumps occur, respectively, with a downstream depth smaller or larg- 
er than the conjugate depth Yo. 

The perfect jump is a powerful energy dissipator. However, if it occurs 
on the unprotected bed, intensive erosive action may take place. In tail- 
water design it is necessary, therefore, not only to create favourable condi- 
tions for the occurrence of a perfect jump, but also to ensure that the zone of 
its erosive effect is well confined within the edge of the solid floor. If the 
jump is repelled out of the floor edge, its location on the floor should be con- 
trolled by the use of baffle piers, sills or stilling basins.(1, 2, 3) These con- 
trol works are generally designed so as to secure a perfect jump. As a safe- 
ty measure, some authors recommend a slight submergence of the jump to 
accommodate for differences between the actual and the assumed hydraulic 
coefficients which may cause the jump to shift downstream from its estimated 
position.(4) On the other hand, low-head river barrages are generally oper- 
ated under conditions of deeply drowned jumps. As energy destroyers, the 
efficiency of such jumps is very low. To make use of the full energy dissipa- 
tion associated with a perfect jump, the floor-level of the barrage is generally 
raised above the bed-level of the downstream channel. Furthermore, a slop- 
ing or stepping floor transition is so designed that a perfect jump is secured 
at all downstream stages.(5) However, such a solution is not always feasible, 
especially in silty streams where high elevations of the floor may hinder the 
passage of the silt loads. 


Models with Movable Beds 


Tests on models with movable beds have often been made to investigate 
the efficiency of anti-scour designs downstream of hydraulic structures. The 
model is usually constructed according to Froude’s criterion of similarity, 
whereas the movable bed is generally kept out of scale. In spite of this par- 
tial distortion in the model, experience indicates that models with movable 
beds, representing high-head structures, have supplied successful designs of 
various stilling basins.(2) Perfect hydraulic jumps are apt to form down- 
stream of these structures and gravitational forces are the most predominat- 
ing of all other involved parameters. 

In a low-head river barrage, where a deeply drowned jump is liable to 
form downstream of the piers, viscous forces may be effectively called into 
action. This is due to the fact that a considerable part of the excess kinetic 
energy in water escapes from the head of a drowned jump and wanders through 
the low-velocity regions, where it is dissipated partly by erosive action and 
partly by viscous shear. Accordingly, both Reynolds and Froude parameters 
may be equally effective. Adding to this combination the non-scaled parame- 
ters which define the bed material, the model results may become uncertain, 
even so far as qualitative similarity is concerned. The geometry of the scour 


2. Notation: The letter symbols in this paper are defined where they first ap- 
pear in the text and are assembled for convenience of reference in the 
Appendix. 
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hole may be described, in this case, by the following general relationship: 


3 a 


a 
in which = is a ratio describing the geometry of the scour hole, F is the Wari 
s a 
Froude number, R is the Reynolds number, length ratios 
describing the geometry of the floor, piers, sills . . . etc., > a > ae 
bm’ Bm 


are dimensionless parameters describing the physical and the hydraulic prop- 
erties of the bed material, C is a constant and @ denotes a function. 

The erratic results obtained from bed-scour tests on open barrages oper- 
ated under deeply drowned jumps have been noticed in several instances in 
Egypt. The scour-protection works in most of the Nile barrages were de- 
signed according to the results of tests on models with movable beds. In all 
these tests, the scales of the models varied between 1 : 100 and 10: 100. Al- 
though the tests were carried out at different dates since 1928, the testing 
technique was almost the same. For any barrage, a model representing two 
or more vents was fixed in a long flume with a glass side. The model was 
kept well above the floor of the flume so that a bed of movable material might 
be provided. The downstream water-level was adjusted by means of a sliding 
plain sill at the lower end of the flume. The results of many tests indicate 
that in spite of the fact that the general features and hydraulic conditions of 
ali the Nile Barrages are much the same, floor-sill locations do not follow 
any general rule. Figures la through f show longitudinal sections of these 
barrages and the dates of original construction and reconstruction of each. 
The various points marked by letters in Figures la through f are identified 
as follows: 


Description 


eee Upstream emergency stop-log grooves 
rere Main grooves fitted with two sliding gates ‘i 
Downstream emergency stop-log grooves a 
Main solid floor 
Upstream bed-protection works 

a eee Lip wall or floor sill 


An example which raises questions as to the validity of the results ob- 

tained from movable-bed tests, under conditions of deeply drowned jumps, 

are those tests carried out by Butcher and Atkinson(6) on a model of the old 
Esna Barrage, Fig. (1f). From the results of those tests, the authors con- 
cluded that, “Under equal conditions of discharge and downstream level, the 
model showed that erosion beyond the end of the solid floor was almost always 
most severe when the water was passed over the gates, rather less so when 
passed under the bottom gate, and very much less when passed between the 
gates.” As a matter of fact, experience indicates that serious tail-erosion 
problems are generally associated with undergate rather than with over-top 
regulation. 
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Before the construction of Edfina Barrage, Fig. (la), earlier bed-scour 
tests showed that the best location of the floor sill was at the middle of the 
solid floor. The same tests were repeated later under the same hydraulic 
conditions, but possibly with different compaction of the non-scaled bed ma- 
terial. The results of these latter tests are shown in Fig. (2). It may be 
noted that the scour holes for some locations of the sill are deeper than the 
hole formed with plain floor. Furthermore, the smallest scour pool is ob- 
tained when the sill is placed at the edge of the floor. It is believed that the 
disagreement between the two test results may be due to the fact that qualita- 
tive similarity of the model is not feasible. The distorting effects of the bed 
material are magnified under the two predominating flow parameters. Ac- 
cordingly, it was decided to perform velocity-distribution tests on a fixed- 
bed model of the barrage. It is believed that more precise conclusions could 
be obtained by comparing the relative effects of various sills on the flow 
distribution. 


Determination of Shear Stresses Along the Bed 


Laboratory determination of the magnitudes of the shearing stresses along 
the bed, which may cause bed erosion downstream of a barrage, requires a 
complete survey of the velocity distribution. The bed of the flume should be 
fixed in order to eliminate the distorting effects of the bed material. Due to 
the separation wakes of the piers, the problem should be dealt with in three 
dimensions. Referring to Fig. (3), and assuming a horizontal bed, the dif- 
ferential equation may be expressed as follows:(7) 


a F, 


3 
+ Sx + = = 0 (3) 


in which Fy, is the pressure difference, b is the width of the flume, T, is the 
shearing stress along the bed, 8 is the momentum corrective factor, P is the 
density of water, Q is the rate of discharge and V, is the forward mean 
velocity. 


The momentum corrective factor at any cross-section is given by the fol- 


lowing equation: 
Y b > 
Vy dz dy 
(4) 
2 


in which v, is the forward velocity at any point and Y is the mean depth of 
flow in the cross-section. 


Assuming hydrostatic distribution of pressures, and re-writing equation 
(3) in a finite form, 


2 


(yg - x2) + 6x (8, - (5) 
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For a qualitative analysis of the relative performance of various sills, it 
is sufficient to investigate their effects on the flow distribution in the vertical 
plane passing by the centerline of the model. The efficiency of any sill 
against bed-scour is measured by the rate of adjustment of the flow to its 
normal distribution in the downstream channel. 


The Model Tests and the Results Obtained 


Fig. (4a) shows a model representing three vents of the Edfina Barrage 
which has 46 vents. The model was constructed of metal to a scale of 1:66.7. 
Fig. (4b) shows the general arrangement of the test facilities. The tail gate 
which controlled the water-level downstream of the model, consisted of multi- 
butterfly vertical ribs which could be rotated by a common header. This de- 
sign eliminated the retarding effect on the bed velocities which might be 
caused by the commonly-used plain sill. A calibrated Pitot-static tube was 
used to obtain velocity measurements at seven vertical sections downstream 
of the piers. The corresponding locations of these sections in the prototype 
are shown in Fig. (la), sections I through VII. An over-head flushing tank 
was used to expel any air bubbles in the rubber-tubing connections. Tests 
were made without floor sills and with sills of various shapes and locations. 
Fig. (5) shows the types of sills used in tests. Figs. (6), (8) and (1la) show 
velocity-distribution diagrams for some of the tests. In each diagram, the 
upper thick arrow represents the maximum forward velocity vmax in the 
section, and the lower thick arrow represents the bed velocity vp. The latter 
velocity was always measured at a height above the flume bed of i th the 
normal depth of flow. It may be noted that the positive areas in the velocity 
diagrams for any one test are not equal. This may be attributed to the wake 
effects of the piers, unsteady flow conditions at the surface, flow vortices 
originating at the bed and to the variation in specific gravity of the water due 
to the entrained air.(8) 

For each test, the following dimensionless plots were drawn: 

(i) and versus , where and Yo are the mean velocity and the 
Vo Yo 


mean depth of flow, respectively, when the flow assumes its normal distribu- 
tion in the downstream channel, and Ly is the distance along the floor mea- 
sured from the pier edges. Examples of these plots are given in Figs. (7b), 


y 
(9) and (11b). (ii) —* and -“M&% versus an where Y is the mean depth of 
Yo Yo Yo 


flow at any section and Wissen the height of the filament of maximum veloci- 


ty above the bed-level of the flume. Examples of these plots are given in 
Fig. (7a). 
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When the flow assumes its normal distribution in the downstream channel, 
the general trends of the functions given in (i) and (ii) become horizontal. Ac- 
cordingly, the plots indicate the relative flow disturbances as well as the ef- 
fects of the different sills. Since there were no appreciable differences be- 
tween recorded and hydrostatic pressures at any section, comparisons were 
not made of pressure distributions. 


Discussion of Results 


. The Methods of Passing the Discharge through the Double Sliding Gates of 
the Vents 


These initial tests were carried out to investigate the effect of the method 
of regulation on the erosive forces. When the discharge is passed under the 
bottom gate, a huge surface vortex is formed, Fig. (6a). Under a constant 
discharge and a variable downstream depth of flow, the extension of this vor- 
tex beyond the pier edges is nearly constant. However, when the discharge is 
either passed between the bottom and the top gates, Fig. (6b), or over the top 
gate, Fig. (6c), the vortex does not appear downstream of the piers. The be- 
havior of the downstream flow for the three methods of regulation is shown 

in Figs. (7a) and (7b). It may be noted that the erosive action is most severe 
when the discharge is passed under the bottom gate and very much less when 
passed either between the gates or over the top gate. It may be also noted 
that the least erosive action is associated with over-top regulation. These 
conclusions are in disagreement with those obtained by Butcher and Atkinson 
from bed-scour tests. (6) 

Fig. (7a) shows the test results for undergate regulation associated with a 
deeply drowned jump. The erratic manner of the flow may be noted from the 
fluctuation of the path of maximum-velocity sheet. Therefore, in all other 
tests the discharge was always passed under the bottom gate in order to mag- 
nify the effects of the floor sills. The working head h on the barrage, (differ- 
ence between upstream and downstream water-levels), was kept constant at 
6.0 cms. (2.46 inches) and the discharge at 8.3 liters per sec. (131.5 gal. per 
min.). These hydraulic conditions correspond, respectively, to a hypothetical 
working head on the prototype of 4.0 ms. (13.12 feet) and to a flood discharge 


through the three barrage vents of 302 cubic ms. per sec. (106,700 cubic feet 
per sec.). 


b. Location of the Floor Sill 


When a solid sill of the type shown in Fig. (5a) is placed between the piers 


or near the pier edges, a bed vortex of a considerable length is formed, Figs. #8 
(8a) and (8b). The strong reversed currents in this vortex may thus take i, 
place over the movable bed and will then participate in the formation of the . 
scour pool. On the other hand, when the sill is placed far from the pier edges, i 


the length of the surface vortex is increased, Fig. (8d). In this case, the de- 
velopment of the normal velocity distribution may be retarded until the flow 
has discharged over the movable bed, sections I through V, Fig. (8d). Since 
the shearing forces along the bed increase with abrupt changes in the velocity 
distribution, the development of the normal distribution of flow on the loose 
bed may give rise to serious local scour. 

It seems reasonable, therefore, that the best location of the sill is indicated 
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by the smallest sizes of both the surface and the bed vortices. In Edfina 
Barrage, this condition is satisfied with a sill located at the middle of the 
solid floor, a distance of 2.5 h from the pier edges, Fig. (8c). 

Fig. (9) illustrates the relative disturbances in the bed velocities for vari- 
ous locations of sill on the floor. The disadvantage of the extreme locations 


of sill may be evidenced from the trends of functions - , sections I, II and V. 


c. Efficiency of Different Types of Sills 


The three types of sills shown in Fig. (5) were tested alternately at sec- 
tions II and V, which correspond to the middle and the downstream edge of 
the floor of Edfina Barrage, respectively. Fig. (10) illustrates the relative 
effects of these sills on the velocity distribution when they are placed at sec- 
tion III. The dentated sill shows slightly better qualities than a solid sill of 
the same peripheral dimensions. This is true for both positions II and V. 
The broad sloping sill, which is sometimes recommended, causes larger 
variations in bed velocities and has no advantage over ordinary types of sills. 


d. The Combined Action of Two Floor Sills 


In order to ensure safety against bed-scour, some of the River Nile Bar- 
rages are fitted with two floor sills placed opposite each other, Figs. (1a), 
(1d) and (1f). The upstream sill, called the baffle sill, is placed between piers 
or near pier edges. Its function is to interrupt flow issuing from the gates 
and to form a jump which dissipates part of the hydraulic energy. The down- 
stream sill, called the lip wall or the floor sill, is intended to deflect the high- 
velocity current upward from the river bed. 

Fig. (1la) shows velocity-distribution diagrams for the baffle sill fixed at 
the downstream emergency grooves and the floor sill placed alternately at 
the middle and the edge of the solid floor, sections III and V, respectively. 
In the first position, a bed vortex occupies the space between the two sills, 
(upper diagrams in Fig. 1la). The water cushion formed thereby is so effec- 
tive in energy dissipation that the flow is adjusted to its normal distribution 
before it leaves the floor edge. In the second position, the size of the bed 
vortex does not change, (lower diagrams in Fig. 11a), but the sill creates re- 
disturbance in the flow over the loose bed. The advantage of a middle-floor 
sill in stabilizing the flow is clear in Fig. (11b). Functions 7 and 2 m 

(trends in full lines), become almost parallel upstream of section V, the edge 
of the solid floor. 


GENERAL CONCLUSIONS 


The behavior of anti-scour sills below open river barrages, operated under 
conditions of deeply drowned jumps, may be clearly investigated by means of 
velocity-measurement tests on models with fixed beds. The results of such 
tests indicate that the best efficiency of sill is obtained when it is placed 
neither too close to the downstream edges of piers nor at the edge of the pro- 
jecting solid floor. This is true whether the sill is used alone or in combina- 
tion with a baffle sill. It is also true in case of dentated sills. The hydraulic 
advantage of these latter sills seems, however, to be over-estimated. 
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The following letter symbols adopted for use in this paper conform essen- 
tially with American Standard Letter Symbols for Hydraulics (ASA-ZIO.2- 
1942) prepared by a Committee of the American Standard Association, with 
Society participation and approved by the Association in 1942. 
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form, with subscripts s, f or m referring to the scour hole, 
the floor or the bed material, respectively; 


= width of the flume; 


= a constant; 


= the Froude number; 


= force component along centerline of the channel; 


= acceleration of gravity; 


= working head on a barrage, the difference between upstream 
and downstream water-levels; 


= length; 


= longitudinal distance along the floor measured from the down- 
stream edges of piers; 


rate of discharge; 


rate of discharge per unit width of channel; 


the Reynolds number; 


= forward mean velocity in a given section; 


= mean velocity in the section where the flow attains its normal 
distribution; 


= forward velocity at any point; 


= forward velocity near the bed of the flume measured at a depth 
0.1 Y, above the fixed bed; 


= maximum forward velocity in a vertical section; 


= mean depth of flow, with subscripts generally denoting the sec- 
tion of measurement; 


= mean depth in the section where the flow attains its normal 
distribution; 


= momentum corrective factor, with subscripts generally denot- 
ing the section of measurement; 


= specific weight of water; 


= density; 


= shearing stress along the flume bed; and 


= a function. 
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PIG.(4a) THE MODEL FIG.(4b) GENERAL ARRANGE- 
MENT OF TEST FACILITIES 


L752) 
(a) Solid sill (bd) Dentated sill (c) Broad sloping sill 


PIG.(5) FLOOR SILLS USED IN TESTS 
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yic.(7) THE BEHAVIOUR OF THE DOWNSTREAM PLOW IN THE THREE 
METHODS OF REGULATION. 
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LOSSES DUE TO ICE STORAGE IN HEART RIVER, NORTH DAKOTA! 


H. M. Erskine,2 M. ASCE 
(Proc. Paper 1261) 


SYNOPSIS 


Winter flows in the lower reaches of the Heart River in western North 
Dakota have been regulated by the Heart Butte Dam since 1949. Flow 
records indicate water losses due to storage as ice is relatively low for 
smaller discharges, particularly if the rate of release is decreased some- 
what as the winter progresses; however, the opposite is true for higher flows 
and increasing rates of release. 


During recent years the flow in the lower reaches of the Heart River in 
western North Dakota has been controlled by the Heart Butte Dam. Winter 
flows have generally been maintained at a rate considerably higher than pre- 
vailed under uncontrolled conditions. Gaging station records at three well- 
spaced locations between the reservoir and the mouth frequently show a de- 
creasing rate of flow from the upper to the lower stations during the winter 
months. Storage in the form of ice in the channel probably accounts for near- 
ly all of the decrease in flow at the downstream points. 

The Heart River, one of the principal tributaries to the Missouri River, 
rises about 30 miles east of the western boundary of North Dakota and flows 
in an easterly direction for about 120 miles to its confluence with the Missouri 
River near Bismarck, N. Dak. The drainage basin of approximately 3,400 
square miles lies between latitudes 46°22’ and 47°11'. (See Fig. 1.) 

Most of the basin is composed of gently rolling prairie although there are 
some relatively large areas with quite steep slopes, particularly in regions 
adjacent to the main stream and its larger tributaries. Altitudes range from 
about 1,650 feet at the mouth to about 2,800 feet at the divide at the upper end 
of the basin. Water courses are well defined and have moderate slopes. It is 


Note: Discussion open until November 1, 1957. Paper 1261 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
83, No. HY 3, June, 1957. 
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2. Dist. Engr., Geological Survey, U.S. Dept. of the Interior, Bismarck, N. Dak. 
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estimated that about one-fifth of the area is under cultivation, principally for 
the production of small grains. The remainder is range land with a firm 
cover of native grass, 

The soils range from shallow to moderate depth and are largely derived 
from the weathering of the underlying sandy clays and shales of the Fort 
Union formation. Exceptions are the bottom lands adjacent to the river and 
its larger tributaries which are composed of alluvium containing much fine 
sand and ranging in width from a few hundred feet to about a mile, relatively 
small areas in the lower part of the basin where some glacial till is found, 
and scattered sandstone capped buttes near the edges of the basin. Figures 2 
to 5 are typical of the terrain along Heart River. 

The average annual precipitation for the basin is approximately 16 inches. 
About three-fourths of this normally occurs during the growing season. Dur- 
ing the winter season precipitation, which is nearly always in the form of 
snow, is usually light. The averages during the winter months range from 
about half an inch of water content per month during December, January and 
February, to about three-fourths inch during March. The average 


- 


Figure 2.--Gaging Station on Heart River below Heart Butte Dam. 
February 28, 1956. 


Figure 3.--Heart River near Flasher. February 28, 1956. 
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Figure 4.--Heart River just below the Lark Gaging Station. 
February 28, 1956. 


Figure 5.--Heart River near Carson. June 17, 1956. Note Ice Scars 
on Cottonwood Trees on Left Bank. 


temperatures during these months are approximately 15, 10, 13 and 26°F. 
respectively. Periods with minimum temperatures of -30°F. or colder usual- 
ly occur each winter. The average annual temperature is about 41°F, 

The average annual runoff for the basin has been about 1-1/4 inch during 
the past few decades. The volume of runoff is very erratic from year to year 
and the bulk of it generally occurs during the spring breakup period. 

Near the end of 1949 the Heart Butte Dam was completed by the Bureau of 
Reclamation. This dam is an earth-fill structure located about 104 miles 
above the mouth. The drainage area above the dam is 1,760 square miles, 
slightly over one-half of the total area of the basin. It is a multi-purpose 
project to provide water for irrigation, flood control, low-flow augmentation 
for stock and other water supplies, waste dilution, and for recreation. The 
gross capacity of the reservoir is 428,000 acre-feet; however, only 68,700 
acre-feet of this is controlled storage held below the crest of the uncontrolled 
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glory hole type spillway. The releases from the controlled storage are through 
a gate controlled conduit. The reservoir first filled to above the spillway level 
in the spring of 1950, 

Ice thickness ranging from 2 to 3 feet are normal on the water bodies in 
this area by midwinter. The volume of ice accumulating on the Heart River in 
the reach below the Heart Butte Reservoir is of considerable importance par- 
ticularly to those in the Mandan area where serious flooding sometimes 
occurs, 

The spring breakup on the Heart River traditionally precedes the Missouri 
River breakup in the Mandan vicinity by about two weeks. Thus in years when 
there are moderate or high flows during the breakup period, the Heart River 
ice is prevented from entering the Missouri River by the heavy ice cover on 
that stream. As a result the Heart River ice accumulates in the channel and 
on the bottom lands for several miles immediately above its mouth obstructing 
the flow of the stream and causing flooding over rather large areas. Although 
levees now protect the city of Mandan against all but extremely severe floods, 
a large area in the vicinity is subject to fairly frequent flooding and consider- 
able damage. It is apparent that any reduction in the volume of ice forming in 
the river would reduce the flood hazard in its lower reaches. 

Three stream-gaging stations are operated by the U.S. Geological Survey 
on the Heart River between Heart Butte Dam and Mandan, and one is operated 
reasonably near the mouth of each of the three principal tributaries to this 
reach of the river (see Fig. 1). Pertinent data relative to these gaging sta- 
tions is given in Table 1. The flow records collected at these stations offer 
an excellent opportunity to evaluate the volume of water that is retained in the 
river in the form of ice under various winter conditions and for different rates 
of stream flow. 


Table 1. --Gaging stations in lower Heart River basin 


Miles 


Stream and location above 
mouth 


Datum of Gage Drainage area 
above in square miles 


sea level 
Heart River below Heart Butte 
Dam ------------------- 103 1,998. 87 1,710 
Antelope Creek near Carson- 4 1,974 221 4 
Muddy Creek near Almont--- 12 1,864 456 J 
(Intervening Area - ungaged)- - - 313 id 
Heart River near Lark -------- 62 1,802. 83 2,750 7 
Sweetbriar Creek near Judson 16 1,886 157 ” 
(Intervening Area - ungaged) - . 403 
Heart River near Mandan ------ 14 1,638.70 3, 310 


Note, --All stations have recording gages, except Antelope Creek 
near Carson, which has a wire-weight gage. 


An examination of the records indicated that portions of the winter periods 
during the four water years 1953 to 1956 were suitable for use in a study of 
the losses due to ice formation. The winter records for each of the six gag- 
ing stations are given graphically in Figures 6 to 9. U.S. Weather Bureau 
maximum and minimum daily temperatures and daily precipitation recorded 
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at the Mandan Experiment Station are also shown for each year. 

The periods selected for study begin as soon as possible after the freezing 
weather started at a time of fairly steady releases from the reservoir and 
with base flow conditions prevailing on the tributaries. The periods were 
ended just prior to the first sign of surface water inflow from snow melt as 
indicated by the hydrographs for the tributaries provided reasonably steady 
flow from the reservoir prevailed for several days subsequent to the ending 
date. This procedure eliminated making adjustments for changes in channel 
storage and snow-melt runoff inflow which could not be adequately evaluated 
with the data available. 

December 16 to March 10 was chosen as the study period for the 1953 
water year. Figure 6 shows that the reservoir outflow was quite low ranging 
from 13 cfs (cubic feet per second) at the beginning to 10 cfs near the end of 
the period. Tributary inflow was reasonably constant. The winter was mild 
with no extremely cold periods and each month several degrees above normal. 

December 6 to January 31 was selected as the study period for the 1954 
water year. Figure 7 shows the reservoir release was steady at 100 cfs at 
the beginning of the period and decreased rather gradually to 90 cfs at the 
end. Tributary inflow declined steadily during the period. The period was 
rather short due to snow-melt inflow developing on February 1. The average 
December temperature was several degrees above normal but the January 
average was about the same amount below normal. 

The period of study for the 1955 water year extended from November 27 to 
March 2, Figure 8 shows the reservoir release was steady at 31 cfs until 
December 5 when it was reduced to from 6 to 10 cfs for about 30 days then 
increased to from 29 to 31 cfs where it was held the remainder of the winter. 
Tributary inflow had a generally declining tendency during the period. Aver- 
age temperatures were within a few degrees of normal during each month. 

Particularly noteworthy during this period is the hydrograph for the Heart 
River near Mandan following the increase in release from the reservoir on 
January 7. It will be noted that the flow at Mandan shows no response to the 
increase in flow from the reservoir until 40 days later and then the increase 
in flow is less than half of the increment at the dam. Apparently all but a 
very small portion of the increase in the release was absorbed in storage in 
the form of ice. 

November 11 to February 29 was used as the study period for the 1956 
water year. The data for this period are shown in Figure 9. Release from 
the reservoir was reasonably steady, varying from 73 to 79 cfs, until January 
16 when it was reduced to 54 cfs for about 12 days, then fluctuated down and 
then up moderately until February 11 when it was held steady at 29 cfs for the 
remainder of the period. This was a cold and prolonged winter with the aver- 
age temperature several degrees below normal every month although there 
were no periods of unusually intense cold. 

Of particular significance on the hydrograph is the position of the Lark 
graph with respect to the Heart Butte Dam graph during the latter half of 
February. This was the first time during this winter period that the flow at 
the lower station exceeded that at the upper. Apparently this happened be- 
cause the reduction in release at the reservoir allowed the water to drop 
away from the ice cover or at least reduced the under pressure on the ice 
cover so that there was no further freezing of the flowing water. 

The runoff in acre-feet for each gaging station for each month or partial 
month during the periods of study are Tabulated in Table 2. An estimate of 
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the inflow from the ungaged tributary area in the reach from Heart Butte Dam 
to the Lark gage, and from the Lark to the Mandan gage was made. This in- 
flow is believed to be extremely small during the winter months and has been 
estimated on the assumption that the runoff rate per square mile from the 
ungaged area was about one-half of the rate from the gaged tributaries. 

The total flow entering each reach is computed and compared with the flow 
at the gaging station at the lower end of the reach for each month or partial 
month and for each complete period. The indicated gain or loss for each 
reach and the total gain or loss between Heart Butte Dam and Mandan is listed 
in Table 2. 

Efforts were made to develop relationships between the indicated water 
losses and some of the other variables but no satisfactory relationship was 
found. 

Figure 10 has been prepared to illustrate the indicated losses in the reach 
from Heart Butte Dam to Mandan. The mean temperature at the Mandan Ex- 
periment Station and the mean discharge at Heart Butte Dam has been shown 
with the plotting of the indicated water loss for each month or partial month. 

The plotting for 1952-53 indicates a relatively small accumulation of ice, 
with the water loss for the 85-day period totaling 1,015 acre-feet. The aver- 
age temperature was about 8 degrees above normal. The rate of release from 
the reservoir was quite low at the beginning of the period and was decreased 
slightly at intervals duringthe period. A comparison with the graphs for 
other years indicates that this combination of conditions was very favorable 
for minimizing ice formation. 

The 1953-54 graph shows a much more rapid accumulation of ice during 
December than in the previous year. The average temperature was several 
degrees higher and the rate of release from the reservoir nearly eight times 
as great. Although the average temperature during January was several de- 
gress below normal, the rate of loss did not increase, probably because of a 
small decrease in rate of release from the reservoir. In a period of 56 days, 
the indicated loss amounted to 2,177 acre-feet. 

The 1954-55 period began with a very low loss during December when the 
release rate was quite low and the average temperatures about 10°F. above 
normal. The rate of loss increased sharply during January and even more 
during February. The average temperature was a few degrees above normal 
during January and a few degrees below during February. The average dis- 
charge at the reservoir was increased during each of these months. The in- 
dicated loss during this period of 96 days totaled 3,649 acre-feet. 

The winter of 1955-56 began early with the reservoir release rate at 80 
cfs. The average temperatures ranged from about 10°F. below normal in 
November to about 4°F. below during February. Indicated water loss rose 
sharply during November and December when the discharge rate was main- 
tained quite steady. The rate of loss decreased moderately during January 
and considerably during February. The average rate of release from the 
reservoir was being decreased considerably during these months. A total 
loss of 5,178 acre-feet is indicated for this 110-day period. 

Conclusions drawn from Figure 10 are complicated by differences in tem- 
perature and other factors. However, we might say that generally the rate of 
loss is relatively low for the smaller discharges, particularly if the rate of 
release is decreased somewhat as the winter progresses. Higher losses are 
usually indicated for higher rates of release as is to be expected of a stream 
of this type which widens rapidly as the discharge increases in the range of 
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flow being considered. Decreases in the rate of flow during a given period 
tend to decrease the rate of loss in all ranges of flow experienced. Increases 


_ in the rate of release add substantially to the rate of loss probably owing to 


water breaking through the existing ice cover and flowing on top where it is 
readily frozen. 

The question arises whether all the indicated losses shown by this study 
are caused by the water going into storage in the form of ice. Possibly some 
flow is lost from the channel by seepage into the ground. There are no diver- 
sions during the winter months and the only consumptive use is a very minor 
amount for stock watering. 

September 25 to November 18, 1952 was selected for study of possible 
seepage losses. During this period there was no appreciable precipitation 
for several weeks, the reservoir release was quite steady for several days 
prior to the beginning and the ending of the period, base flow conditions pre- 
vailed on the tributaries, evaporation losses were low and there was no loss 
due to transpiration, there were no known diversions for irrigation or other 
purposes, and there were no complications due to freezing. 

The stream flow and weather records from August to December 1952 are 
shown on Figure 11. The reservoir release was steady at 76 cfs at the be- 
ginning of the study period and was steady at 110 cfs for several days prior 
to the end of the period. The release was reduced to about 25 to 30 cfs for 
several days on two occasions during October. Tributary inflow increased 
slightly during the period. 

The data for the study period are in Table 3. The method of computation 
used was the same as in the winter water loss studies (Table 2) except that 
an evaporation adjustment has been included. This adjustment was computed 
using an average river width of 75 feet and evaporation depths of 0.8 inch, 

2.0 inches and 0.5 inch respectively for the portions of September, October, 
and November included in the study period. Inflow from the ungaged inter- 
vening areas was computed on the same basis as it was for the winter periods. 

The tabulation for the total period shows a 99 acre-feet loss in the Heart 
Butte Dam to Lark reach but a 697 acre-feet gain in the Lark to Mandan 
reach. The net total for entire reach from Heart Butte Dam to Mandan indi- 
cates a gain of 598 acre-feet, approximately 6 percent of the Mandan flow. 
During the period, November 1-18, when the flow at all stations was almost 
constant it is noted that the gain at Mandan was 60 acre-feet, less than 2 per- 
cent of the Mandan flow. The reason for this indicated gain cannot be posi- 
tively stated; however, it is probable that it results from seepage inflow from 
the alluvium forming the valley fill in the lower portion of the basin. 

This study indicates that probably none of the indicated water losses dur- 
ing the winter months are due to seepage from the river. On the contrary, it 
appears that the losses ascribable to ice storage may actually be greater than 
indicated in Table 2 because seepage inflow is disregarded. 

It appears proper to conclude that the losses due to ice storage are at least 
as great as indicated by the computations. These ice losses are, of course, 
temporary as far as total runoff is concerned because they are recovered 
during the breakup period. 

On February 28, 1956, the writer examined the Heart River at a number of 
sites in an effort to obtain field data upon which to base estimates of the vol- 
ume of ice in the river channel. Most of the rural roads were blocked by 
drifted snow but the river was examined at four highway crossings in the 
reach being studied. The ice thickness was measured near each bridge and 


y 
q 
q 
a 
“a3 
a 

= 
4 


AUGUST TO DECEMBER 1952 


MAIN STEM GAGING STATIONS 
——HEART RIVER BELOW HEART BUTTE DAM, N. DAK. 
---HEART RIVER NEAR LARK,N. DAK. 

HEART RIVER NEAR MANDAN, N. DAK. +— 


! 


IN SECOND-FEET 


DISCHARGE 


+ + 


+ + + 


| 
| 


TRIBUTARY GAGING STATIONS 
| ANTELOPE CREEK NEAR CARSON, N. DAK. 
ae CREEK NEAR ALMONT, N. DAK. 
CREEK N.. DAK. 


IN SECOND-FEET 


+ 


Sots. 


DISCHARGE 


MAX. & MIN. 


Pi- 
TATION | TEMPERATURE 


+ - 
5 10 15 20 25 | 0 2025| 5 10 5 2025 10 1582025 | § 10 15 20 
AUGUST SEPTEMBER OCTOBER NOVEMBER DECEMBER 


FIGURE II 


PREC 


ASCE ERSKINE 1261-15 
a LOWER HEART RIVER BASIN DISCHARGE RECORDS | | at 
| | | ] | 
| | 4 4 
10 
| 
10 4° 
| 
| | | | 
= 
| 
4 
a 
all 


uepueW 03 
Wed SSOT HO NIVD 


June, 1957 


97 4787 SSOT HO NIVOD 
UPPUEW 
TeqOL 
uot 
eely 
uospnye year) 
Jeeu 
04 
wed SSOT YO NIVOD 
JeATY 
uot 


ZuTUeAIE 


quowTy seeu Appny 


uosse) yeerd 


wed MOTeq JeATY 


| 
1261-16 HY 3 
w N ao oO an w w a 
i & g 
£ 
¢ 2 ae 
w 
15 
2 
© N i a R i 3 
a 
a 
2 
° 
a 
= ; 
x an 
3 
° 
el § 
2 = 
oa 
»| 2 
fa) 
3 
' 
a F 
= 
e| @ 4 
7 A 
i 
& 


ASCE ERSKINE 1261-17 


the width of the ice cover was estimated by pacing at a number of sections 
within several hundred feet above and below each bridge. The ice thicknesses 
measured ranged from 2.4 to 2.9 feet and the stream widths ranged from 55 
to 110 feet with most width measurements in a range from 90 to 100 feet. On 
the basis of these samples an average ice thickness of 2.6 feet and average 
width of 90 feet was assumed for the total reach. The ice was found to begin 
about 2 miles below Heart Butte Dam. A computation of the ice volume in the 
reach between the dam and Mandan based on this information indicates an ice 
volume of 2,500 acre-feet. The computations based on the streamflow records 
indicate that the volume at this time was about 5,100 acre-feet, or practically 
twice as much. 

It is believed that the discrepancy between these two figures is largely due 
to the samples of ice thickness not being representative of the overall reach. 
A process known locally as glaciation is believed responsible for this situa- 
tion. During this process, water breaks through the ice from underneath or 
enters from the sides, flows over the ice and freezes upon it, thus causing 
abnormal thickness. Hydrographers working on this river as well as others 
in the general area, sometimes encounter ice thickness of from 4 to 6 feet 
when searching for sections for making current-meter measurements, It is 
not uncommon to find that water flows from holes cut in the ice in connection 
with current-meter measurements or gage readings, indicating that the water 
is flowing under appreciable pressure and would rise above the ice surface 
when cracks or other openings permit. Ranchers occasionally report cutting 
through 5 and 6 feet of ice to obtain water for stock. Where this glaciation 
process has been operative, it is probable that the width of the ice is sub- 
stantially greater than indicated by the samples used in the computation. The 
width sampling did not take into account the ice that forms in the mouths of 
the numerous small tributaries as the ice formation on the main stream 
forces backwater into these areas. 

These circumstances lead to the conclusion that a volume of ice storage 
computed from a few samples of ice thickness and width is likely to be un- 


reliable, and that the volume of ice storage should be determined by water 
loss studies, 
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ATTENUATION OF SOLITARY WAVES ON A SMOOTH BED 


Yoshiaki Iwasal 
(Proc. Paper 1262) 


SUMMARY 


The hydraulic characteristics of solitary waves have been considered in a 
preceding paper, (1) In this paper, further development of the study of solitary 
waves is treated. Especially considered is the mathematical analysis of the 
attenuaiion process of solitary waves. The theoretical results are compared 
to experimental data of Russell and Ippen. 

It should be also noticed that the total kinetic energy minus the potential 
energy is twice the amount of the kinetic energy of the vertical motion, as 
Starr(2) has already verified. 


I, INTRODUCTION 


Theoretical and experimental researches on the solitary wave, defined as 
a single elevation of water surface above the still water level moving with 
esha | celerity without change of form, may be found in many fruitful refer- 
ences.\", 

Although the purpose of the study on this problem was originally based on 
the mathematical interest, regardless of the phenomenon’s actual existence, 
it is well known that a close analogy exists between the hydraulic character- 
istics of a solitary wave and those of an incoming water wave in surfs near 
coasts. Therefore, there has recently been a renewed interest in obtaining 
the complete characteristics of solitary waves in order to apply them to the 
problems of coastal engineering.(4, 

The previously published papers(1, 6) were concerned with the hydraulic 
configuration of open channel flows with appreciable vertical acceleration. 

It was shown that the solution of a second order approximation of translation 


Note: Discussion open until November 1, 1957. Paper 1262 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
83, No. HY 3, June, 1957. 

1. Asst. Prof. of Hydraulics, Dept. of Civ. Eng., Faculty of Eng., Kyoto Univ., 

Kyoto, Japan, visiting fellow, Hydrodynamics Lab., Massachusetts Inst. of 

Technology, Cambridge, Mass. 
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wave flow was expressed in terms of Jacobi’s elliptic function including the 
solitary wave as one of the limiting cases, which also revealed some funda- 
mental hydraulic characteristics of the solitary wave. 

In this paper, further development of the study of solitary waves will be 
presented. Especially considered is the attenuation process of solitary waves 
by viscous shear within the laminar boundary layer between the main poten- 
tial flow under a passing solitary wave and a smooth horizontal channel bed. 
The mathematics of the attenuation problem was first treated by G. H. 
Keulegan(7) in 1948, based on the small amplitude theory of wave motion. In 
1953, J. W. Daily andS. C. Stephan(4) published an empirical formula for the 
attenuation, a result of much experimental data, and in 1955 and 1956 A. T. 
Ippen and G. Kulin(8,9) published results of further experiments and pre- 
sented a new approximate attenuation formula based on the Blasius’ formula 
of coefficient of friction experimentally verified by their experiment. 

The finite amplitude theory of solitary waves is here considered in the 
treatment of viscous damping in the laminar boundary layer. It should be 
noticed, therefore, that the total kinetic energy is larger than the potential 
energy by twice the amount of the kinetic energy of the vertical motion, as 
V. P. Starr(2) first verified in 1947. Assuming that the flow within the lami- 
nar boundary layer under a solitary wave is also translatory, as is the main 
flow outside of the boundary layer, and the velocity distribution is given by 
some suitable formulas, the viscous damping of a solitary wave may be 
evaluated by using the usual concept that the rate of energy loss of the wave 


is equivalent to the rate of dissipation of energy in the laminar boundary 
layer. 


Il. HYDRAULIC CHARACTERISTICS OF SOLITARY WAVES 


As described in the previous section, some hydraulic characteristics of 
solitary waves are now derived, considering translation wave flow with ap- 
preciable acceleration in the vertical motion. In this section hydraulic char- 
acteristics of solitary waves are presented for the convenience of discussion 
in a later section. But it should be prerequisitely noticed that in the theory 
discussed the velocity distribution in the horizontal direction is assumed to 
be truly translatory. 
The solitary wave profile is given as the Rayleigh’s type 


aysech” | (x-ct)| (1) 
where 7 = elevation of wave profile above still water level, a, = amplitude, 


Yo = Still water depth, x = horizontal distance along the wave measured from 


the wave crest, t = time, and C = wave celerity expressed by Vg(yg + ap) 
from consideration of the law of continuity. 


After transforming the coordinate system into the moving one, 


Zz 


Putting yo >> a), the expression of Boussinesq and of Keulegan for solitary 
waves is obtained. 
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The horizontal velocity in the main flow is given by considering the law of 
continuity, 


7 


and the vertical velocity is given by using the usual hydrodynamic concept of 
stream line, (1) 


where y = vertical distance from the channel bottom. 
The volume of the solitary wave per unit width is obtained by integrating 
” from negative infinity to positive infinity. 


= 


Now, consider the energy of solitary waves: in the small amplitude theory 
of wave motion, the energy is half potential and half kinetic. This forms the 
basic idea of the study on attenuation by Keulegan(7) as well as Ippen. (8) In 
1947, Starr(2) first integrated the energy integral for the periodic gravity 
wave anu applied it to the problem of solitary waves. His results explain that 
the kinetic energy is larger than the potential energy by twice the amount of 
kinetic energy of the vertical motion. But it is rather difficult to discuss 
quantitatively the hydraulic characteristics of solitary waves in Starr’s deri- 
vation. The same conclusions as Starr’s have been obtained by using the 
theoretical results of reference(§) and are as follows. 

The potential energy is given by 


where p = density of fluid, and g = acceleration of gravity. Substituting 77, 
expressed by Eq. (1), and integrating 


The kinetic energy of the horizontal motion is 
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The integration yields 


The kinetic energy of the vertical motion is given by 


Exus= P| 2) dx - 


Again integration yields 


The total energy of a solitary wave per unit width is expressed by the sum 
of Eqs. (5), (6), and (7), that is, 


+ I+ 


The above equation is the total energy of the finite amplitude wave. When the 
wave amplitude is very small compared with the original water depth, 
Yo >> ao, Eq. (8) becomes 


(9) 


Eq. (9) is the same as Keulegan’s expression, and it follows that the potential 
energy is equal to the kinetic energy in the small amplitude theory. Denoting 
E,, as the total kinetic energy, the following relation between the kinetic and 
potential energies is obtained: 
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This relation shows that the kinetic energy in the finite amplitude theory is 
larger than the potential energy, and the difference between them is twice the 
amount of the kinetic energy of the vertical motion, as Starr first derived in 
his gravitational wave thecry. Figure 1 shows the relation between the di- 
mensionless energy, expressed as E' = E/p Yo» and the relative amplitude, 

oor rate of energy dissipation when the solitary wave passes through the 
channel is given by differentiating Eq. (8) with respect to the time, that is, 


where s = travel distance of a wave. 


The above relation is one of the fundamental equations in evaluating the damp- 
ing process of solitary waves in the laminar as well as turbulent boundary 
layer. 


Ill, RATE OF ENERGY DISSIPATION IN THE LAMINAR 
BOUNDARY LAYER 


The rate of dissipation of energy per unit volume, @ , in the two-dimen- 
sional laminar boundary layer is given by 


2 


where u and v = horizontal and vertical velocities in the boundary layer, 
and # = dynamic viscosity of the fluid. It is sufficient to consider that the 

gradient of horizontal velocity in the y-direction is very large compared with 
other terms, thus the principal part of the above equation becomes 


(13) 


Damping of a solitary wave on a smooth channel bed is usually considered to 
be due to the viscous shear expressed by the above equation. If the velocity 
distribution in the laminar boundary layer is determined from the original 

Navier-Stokes equations of motion, the rate of energy dissipation can be ob- 
tained by integrating from the channel bottom to the thickness of the boundary 
layer, and from negative infinity to positive infinity. It is, however, very dif- 
ficult to integrate the original hydrodynamic equations. Although Keulegan 
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Relation between energy and 
relative amplitude 
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did, his theory implies the assumption that the thickness of the boundary 
layer extends to infinity, while the flow is actually confined. Under some ap- 
propriate conditions, the original dynamic equations become linear, anda 
rigorous solution is obtained, but the expression is so complicated, by being 
composed of a Fourier series, as to be impractical. This difficulty is mainly 
due to the finite thickness of the boundary layer. 

Ippen, Kulin, and Raza(8) first measured the velocity distribution in the 
laminar boundary layer by an electrical differential gauge, reported in 1955, 
and their results show that the velocity distribution in the boundary layer can 
be assumed to be linear or parabolic as in the steady regime. In this paper, 
the distribution of velocity is also assumed to be given by the Couette 
(Linear) and parabolic types of flow. Thus, by integrating the above equation, 
the rate of energy dissipation per unit width is given by, 


for Couette flow: 


for parabolic flow: 


(15) 


where § = thickness of boundary layer. It is, therefore, necessary to evalu- 
ate the process of development of a boundary layer under a passing solitary 
wave. 

As the thickness of the boundary layer may be considered to be rather thin 
compared with the main potential flow, except near and at negative infinity, it 
is assumed that the flow is unconfined. Hence, the development process of 
the boundary layer is simply given by the following equation (10): 


(20 +3") +c (16) 


Ue) ar Ox 


(ur) = 28 20, | 


where u* = frictional velocity, and @ and §* = momentum and displace- 
ment thickness defined by 


Transforming the coordinate system into the moving one with constant 
celerity 
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Assuming that the velocity distribution in the boundary layer is given by 
Couette and parabolic types of flow, the momentum and displacement thick- 
nesses become, respectively, 


for Couette flow: 


& 
» 


for parabolic flow: 


25 
= 75? 


Hence, the following boundary layer equations, substituting the above values 
into Eq. (17), are obtained, respectively, 


2 


SC- Sule 52 _ (18a) 
ox 
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Transforming the independent variable from x to 7 and introducing the di- 
mensionless parameters 2=Q0, § =Qi, and Eqs. (18a) and 
(18b) become 


2 + 12 C1 +A) (19a) 


where Rc = Reynolds number of celerity, determined by the given flow con- 
dition expressed by gi/2y,3 2/y , the upper positive sign means the equation 
for the positive direction from the origin to the infinity, and the lower nega- 
tive sign, the equation for the negative direction. As it is still difficult to 
solve the above equation, the following equation, as the first approximation 
for small values of Ag, can be assumed for parabolic flow: 
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Although this equation has a singular point at o = 0, the solution for the ap- 
proximate equation under the boundary conditions i, = 0 at o= 0 for the 
positive direction and iz, = i, at o = 1 for the negative one becomes 


2 


2.2. o 

tz Ao [/+ (/+ 9 

where subscripts 1 and 2 show the values in the positive and negative direc- 


tions, respectively. This is the approximation equation for the development 
of laminar boundary layer under a passing solitary wave. Figure 2 shows an 
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Fig. 2 Development of the 
boundary layer 
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example of the development process of the boundary layer. The above equa- 
tion implies that at negative infinity the thickness of the boundary layer be- 

comes infinite although the actual flow is confined. Hence, Eqs. (21) may be 
applied to the present problem from positive infinity to the point at where the 
thickness of the boundary layer is equal to the water depth. This point is cal- 
culated by the following relation, with sufficient accuracy, for parabolic flow: 


-4 
e~ 
‘A 


In the same way, the critical point for the case of Couette flow is given by: 


“az 


From this point to negative infinity, the thickness of the boundary layer will 
be considered to be equal to the water depth, that is, the fully developed 
laminar regime will be attained. Although the behavior of the development 
process in the boundary layer at the negative branch in the present analysis 
seems very doubtful, its actual behavior cannot be confirmed. 

Applying Eqs. (21) to this problem, the thickness of the boundary layer 
under the crest of a solitary wave is given by, 


for parabolic flow: 


2 
4? Qa. , 


for Couette flow: 


The theoretical thickness of the above expressions are very small, and about 
one-fifth of the experimental results presented by Ippen, Kulin and Raza. 
More detailed experiments will be needed to clarify the development process 
of the laminar boundary layer. 

Therefore, the expression of the rate of energy dissipation in the laminar 
boundary layer for parabolic flow from Eqs. (2), (15) and (21) is 
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For the case of Couette flow, the following relation of the rate of energy 
dissipation is derived in a similar manner: 


1 


(23) 


These above two relations are also fundamental equations in evaluating the 


process of viscous damping in the laminar boundary layer under a passing 
solitary wave. 
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IV. ATTENUATION FOR A SMOOTH BED 


As already shown in the previous section, the basic principle of attenuation 
of solitary waves is that the rate of energy dissipation in the laminar boundary 
layer under a passing solitary wave must be equal to the rate of energy loss 
manifested in the amplitude decrease of the wave. For the present case of 
parabolic flow, the following relation, from Eqs. (11) and (22) is obtained: 


(24) 
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Integrating the above equation under the beundary condition of (a/yo) = 
(a;/y,) at s = yields, 
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This is the two-dimensional attenuation formula based on distance traversed. 
If the channel is of limited width, however, some factor to include the fric- 
tional effect along the side walls must be added. If (1+ 2y,/B) is used as a 
first approximation, which is only correct for very low waves, the following 
relation for the rectangular channels is obtained: 


(26) 


20027 = 0.1499 Re (1+ 


In the same way, for the case of Couette flow, the attenuation equation is 


(27) 


(I+ BX) 


V. DISCUSSION ON THE ANALYTICAL RESULTS 


To verify the above theoretical approach, classical experimental data of 
Russell as well as recent data of Ippen will be considered. In Figure 3, these 
experimental data are plotted with the parametric notation A j: In Russell’s 
data, the temperature is assumed by Keulegan to be 15.8°C, oie in Ippen’s 
data the kinematic viscosity is assumed to be 1 x 107° ft”/sec. Although the 
experimental data are very scattered, the value of K' 


seems to increase as the ratio of Aj to A, increases. K' is not constant, 
but a variable of the amplitude as well as, rigorously speaking, the Reynolds 
number of celerity in this theory. It is constant for the small amplitude 
theory, that is, 1/12 for Keulegan’s expression, and 1/14 for Ippen’s expres- 
sion. Curves of K' versus Aji/Ao expressed by parametric values of Aj 
are also shown in the same figure. The dotted lines show the relation of the 
attenuation equation when Couette flow is assumed. Comparing these curves, 
it seems that Eq. (26), which represents the attenuation equation of a solitary 
wave in the laminar boundary layer by assuming parabolic distribution of 
velocity, is considered to be fair. More approximately, Ippen’s expression 
based on the Blasius’ formula of the coefficient of friction seems better than 
Keulegan’s equation. 

Although this approach to attenuation implies some physical assumptions 
which will influence the development process of the boundary layer to a 
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Fig. Relation between K’and Ai/Arc 


considerable extent, it cannot be ascertained whether these assumptions are 
correct or not, because the experimental technique is still too primitive to 
obtain detailed behavior in the boundary layer. After establishment of 


further progress in experimental technique, a discussion of these assump- 
tions will become helpful. 
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: solitary wave height above still water level 


: channel width 


: wave celerity 


--¢ : subscript “c”, referring to critical condition except Rc 
E : Wave energy per unit width 

Ek =: kinetic energy per unit width 

Eku : kinetic energy of horizontal motion per unit width 

Ekv : kinetic energy of vertical motion per unit width 


Ep =: potential energy per unit width 


E' : dimensionless energy of wave expressed E/ Bo. 

g : acceleration of gravity 

--; : subscript “i”, referring to initial wave in attenuation theory 
K' : damping coefficient in attenuation 

--g : subscript “o”, referring to main flow condition 

Re : Reynolds number of celerity expressed by 


: horizontal distance of wave travel 


: time 


: horizontal particle velocity 
: vertical particle velocity 

: wave volume per unit width 

: horizontal distance along wave measured from wave crest 
: vertical distance from channel bottom 


: subscript “1”, referring to values in positive direction of x 


--2 : subscript *2”, referring to values in negative direction of x 
: thickness of boundary layer 
: displacement thickness 
surface elevation measured from still water level 


: momentum thickness 


: dimensionless thickness of boundary layer expressed by 3 /ao 
: dimensionless wave height expressed by ad/Yo 

: dynamic viscosity 

: kinematic viscosity 

: density 
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: dimensionless surface elevation expressed by 7) /ao 


rate of dissipation of energy by viscous shear per unit volume 
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THE ROLE OF SEDIMENTATION IN WATERSHEDS 


Fred H. Larson,! M. ASCE, and G. Robert Hall2 
(Proc. Paper 1263) 


SYNOPSIS 


Federal participation in watershed programs depends upon the premise 
that benefits must exceed the costs of any proposal. A large part of the bene- 
fits are usually due to the reduction of flood water and sediment damages. 
Since program proposals may have widely divergent effects upon sediment 
and flood water damages, it is necessary that separate appraisal systems be 
set up for flood water damages and for sedimentation damages. 

This paper outlines the approach to sediment problems being made by the 
Department of Agriculture. A method is outlined to inventory the damages 
due to sedimentation and the benefits related to the reductions in these dam- 
ages attributable to watershed programs. 


Authority for the installation of watershed works of improvement by the 
Department of Agriculture has come in three stages. These are the Flood 
Control Act of 1936 and appropriations for “pilot watersheds” approved in 
1953. The third, and the latest, legislative authority for this type of work is 
Public Law 566, the Watershed Protection and Flood Prevention Act of 1954, 
as amended in 1956. 

There are 11 watersheds authorized for treatment under the Flood Control 
Act and 58 authorized as pilot watershed projects. 

The basic procedure used in developing these projects is also being used 
in developing plans for watersheds under Public Law 566. This law states 
that it is the “policy of the Congress that the Federal Government should co- 
operate with local public agencies to prevent erosion, flood water and sedi- 


ment damage and further the conservation, development, utilization and 
disposal of water.” 


Note: Discussion open until November 1, 1957. Paper 1263 ia part of the copyrighted 
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83, No. HY 3, June, 1957. 


1. Head, Eng. and Watershed Planning Unit, Soil Conservation Service, Upper 
Darby, Pa. 
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In the language of the Act, erosion, floodwater and sediment damage are on 
an equal footing. This is a significant reflection of the change in emphasis 
that has developed in watershed and flood control programs during the past 
generation. Prior to World War II, concept of flood control seldom went be- 
yond the protection of specific localities against inundation. 

In the watershed studies made by the U. S. Department of Agriculture, 
sediment has been considered an integral part of the problem on the grounds 
that it is deposited largely by flood flows or during flood-producing storms. 
Sediment damages are not integrated with water damages, sediment benefits 
are computed independently. 

The water conservation program of the Department of Agriculture empha- 
sizes the proper use of every acre. The widespread practice of strip- 
cropping, contour cultivation, terracing and other soil and water conserva- 
tion measures by landowners and operators reduce sedimentation. This 
effect is related to but not solely a function of the hydraulics of flow. 

Benefits from watershed treatment must be related to damages. This type 
and extent of the treatments proposed are predicated upon the benefits. 
Therefore, watershed projects usually start with a resource inventory, in- 
cluding a sediment damage survey. When damages are tabulated, the nature 
and scope of the proposed plans are indicated. The anticipated costs and 
benefits of alternate plans are then determined. 

Perhaps the clearest way to present the Department of Agriculture’s ap- 
proach to sediment problems in watershed planning is to outline the proce- 
dures used in developing this phase of the watershed plan. 

The investigation consists of four phases: 


1. Determination of damage from sediment. 

2. Determination of sediment sources and yields. 
3. Developing the plan of treatment. 

4. Benefit computations. 


In the sediment damage surveys it may or may not be important to go 
further into the type of damaging sediment than mere bulk properties would 
indicate. In many watersheds damages may be related to type of sediment, 


which in turn may be related to a source. This may isolate an area where 
control should be specified. 


Bottomland Damage Survey 


The bottomland damage survey is the starting place. Over the years it 
has developed into a rather standard procedure. A classification of the bot- 
tomland is made. This classification is usually by “damage reaches.” A 
damage reach may be an area of high sediment damage or it may be based on 
a hydraulic function, related cultural features or other considerations remote 
from sediment damages. Since the sediment damage study must be coordi- 
nated with other studies, it has been found more efficient to coordinate these 
studies from the beginning rather than interrelate them at a later date. 

Once the “reaches” have been selected, the geologist, by a random sam- 
pling method, selects the locations of the valley cross sections. These may 
not be the hydrologic cross sections. 

The bottomland damage survey is a field job. The geologist proceeds 
from valley wall to valley wall, recording the length along the section and 
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the percent of damage or benefit of each separate occurrence. The standard 
form (Fig. 1) provides for five types of damage to the land and an estimate of 
the percent of damage caused by each separate occurrence. Overwash, 
swamping, streambank erosion, floodplain scour and valley trenching are the 
types of damage enumerated. 


Overwash 


Infertile overwash is the deposition of relatively infertile sediment upon 
an old soil. It may take place on alluvial or colluvial soils and may take vari- 
ous forms, such as overbank splays, fans, or vertical deposits. Sterile sand 
or gravel commonly cause the damage, but silts and clays derived from sub- 
soil are harmful in some areas, especially when deposited on heavy soils. 

It is possible to have a soil improved by the deposition of more fertile 
overwash, or by changing the texture. For instance, a silty deposit may im- 
prove a very sandy soil. These areas are appraised and the benefit weighted 


against the damage on other areas to obtain an estimate of overwash damage 
in the watershed. 


Swamping 


Swamping is any impairment of drainage of bottomland or colluvial soils 
by sediment deposit. It may be caused by filling of stream channels with the 
products of erosion, which raises the watertable on the bottomlands, or 
formation of natural levees by sediment deposits which prevent proper sur- 
face drainage. Deposits of clay upon bottomland soils cause “puddling,” re- 


duce permeability and prevent internal drainage. 
In some cases a benefit, rather than a damage, may be found, such as 
when sand is deposited on very heavy textured soils. 


Streambank Erosion 


Streambank erosion is the recession of the bank of a stream. This reces- 
sion of the bank is usually accompanied by a replacement of material on the 
opposite side of the stream. 

In estimating damage, a comparison is made between the production of the 
replaced material and the uneroded area. Interviews with farmers, tree 
growth, and photogrammetric evidence help estimate recovery and the time 
involved. 

The “net” cross sectional area of material removed is recorded. Under 
fairly typical conditions the eroding bank may be six feet high, tapering to 
the six foot depth over a length of 30 feet. In that case, for example, the 
cross sectional area recorded would be 90 square feet, although the gross 
volume transported by the process would exceed that amount. 


Floodplain Scour 


Scour includes localized scour channels on the floodplain, sheet removal 
of large areas of soil from the bottomland, and erosion of terrace soils by 
flood flows. 

Since evidence of scour damage is often eliminated by farm operations, it 
is difficult to express this damage on a historic basis. Landowners observa- 
tions are an invaluable guide to this decision. 
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Often it is found that the amount of scour determined by the survey is 
typical of the year-to-year occurrence. In this case the damage is assumed 
to be the average annual damage. 


Valley Trenching 


Valley trenching is the term used to describe a gully on the bottomlands. 
It may be caused by a concentration of overbank flows. This in turn causes 
an eroding channel or a side tributary that partly or entirely abandons its 
channel and erodes a new passageway along or across bottomland to the main 
stream. 

Valley trenches are often transitory. Again, landowners observations are 
valuable as a guide to their history and rate of growth. 

The percent of damage caused by these occurrences is a determination the 
field man must make on each damaged area. His criteria for the damage 
estimates are developed with experienced soil and crop specialists in the 
area. These specialists help him determine, for instance, how much a six- 
inch deposit of sand will affect production from a particular bottomland soil. 
The combined judgments are based on pot studies (where yield comparisons 
are made between the sediment and unaffected bottomland soils), plot yield 
studies (where the yield from an affected plot is compared with yield from an 
unaffected plot often in the same field), and chemical analysis. 

Table 1 is an example of the guide developed to aid the field men in their 
determinations. 

Unfortunately, a great many of our watersheds offer few opportunities to 
compare affected areas with unaffected areas. Surprisingly few unaffected 
bottomland areas exist today. 

In summarizing data for the bottomlands, the lineal distance of damage is 
totaled and the “weighted average” of the damage is computed for each range. 
The several ranges for each reach are then summarized. In short, an 
“average” range for each reach is developed. 

The next step is to expand the damage to the reach of bottomland in ques- 
tion. If bottomland area is available from aerial photographs, the average 
range length is divided into the area to obtain a valley length factor. If acre- 
age is not available the mileage of stream valiey in each reach is measured 
from the best maps available. 


The acres affected by each occurrence are then computed by the formula: 


a=0.121 bc 

where a = acres affected 
b = av. lineal distance on each range, in feet 
c = mileage of valley in this reach 


When this process is completed we have developed an inventory of the 
bottomland. It shows the area affected to date by these different types of 
damage and the amount of damage, percentage-wise. 

This inventory summarizes conditions as they exist on the bottomlands in 
the watershed. It-is a total occurrence to date of the area and extent of the 
damage, or benefit, caused by sediment and erosion on the bottomlands. 

Due to the widely variable land use histories and physiographic conditions 
found in the watersheds, no single method has been found for converting the 
total damage to an average annual basis. 
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Table 1 


Soils Laboratory Analysis of Recent Sediment 


and Old Soils From Chigley Sandy Creek Valley 
Washita River Watershed, Oklahoma 1/ 


Available Available : 
Sample No. : pH > Phosphorous : Potash > Matte ‘Sand :Silt:Clay 


(1) Sediment 8.2 (basic) 0 (low) 


(2) Sediment 


8.5 (basic) 0 (low) 


(1) Old Soil 


6.6 (s. acid) 10 (med. ) 


(2) 01d Soil 


7-1 (neutral ) 20+ (high ) (high) 2.5 (high) 


Percentages o 


to Valley Lands as Related to Depth and Texture, 
Chigley San ] 


k Valley, Washita River Watershed, Oklahoma. 1/ 


Depth Increment Type of Sediment Percent Damage 


4-8 inches Fine and coarse sand and s 


4-8 inc 
or 
8-12 inches 


sand 


sand and gravel 


edimentation Damages and Sediment Contribution Rates 
eek Drainage Area, Washita River Basin, Oklahoma, 
-ontrol work plan (unpublished). USDA-SCS, Ft. Worth, 


in Chigley Sandy C 
a part of a flood 


Texas. 


1/ Marshall, R. L. - S 
r 


ub 
: : : > Percent : Percent Bes 
a 
162 (med.) 1.74 (med) 79 21 0 
e202 (high) 2.3 (high) 33.5 39.5 27 
Table 2 J 
a 
20 
Medium and coarse sand 40 
Fine and coarse sand Xe) 4 
12-24 inches Coarse 60 
12-24 inches Coarse 


\ 
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In the case of erosion damage we are fortunate in universally available 
aerial photographs. Many of the watersheds have been photographed two or 
more times, with a time lapse of 8 to 12 years between coverages. In areas 
with two sets of photos, it is possible to measure the changes in channel 
alignment that occurred between the flights and thus obtain reasonably accu- 
rate measurement of the amount of streambank erosion that occurred. This 
readily gives an average annual rate. Where only one flight is available the 
rate can be determined by ground measurements to the cutting edge of the 
bank from carefully identified points on the photos, and on the ground. It 
might be pertinent to mention that the line measured on the ground is a 
straight line between two identifiable points, one left and one right of the 
channel. This permits an exact rectification of the photographic scale on 
this line. 

Photographic measurements are used only to determine the average annual 
rate of occurrence. This is compared with the “gross” area affected as de- 
termined in the valley land damage study, primarily for the purpose of re- 
lating the “gross volume” of sediment produced, as determined by the bottom- 
land damage survey, to an average annual volume of sediment produced from 
this source. 

There is usually a degree of natural recovery on areas of bottomland 
damages. In monetary evaluation, the geologist provides the economist with 
data on the total area damaged and the loss in productivity. The estimate of 
the portion of the damage that could be expected to recover were it no longer 
subjected to flooding or deposition is also given. The economist can evaluate 
the annual loss in income from the damaged area based on the present dam- 
age. The loss in value of production from the non-recoverable portion of the 
sediment damage will continue after the installation of a program. The re- 
duction of deposition as a result of the program is applied to the recoverable 
portion to determine the benefits from the reduction of sediment damage ac- 
cruing to the program. 


Other Sediment Damages 


A number of sediment damages are appraised in addition to bottomland 
damages. Not all of these exist in any one watershed, however, their im- 
portance must be considered in any watershed survey. 


Reservoir Sedimentation 


This damage is the loss of storage capacity by sedimentation of water 
supply, power, recreational, desilting, flood control and other reservoirs. 
The first step in evaluating reservoir sedimentation damage is the com- 
pilation of a list of reservoirs in the watershed, by physiographic provinces 
or problem areas, including data on capacity, size of drainage area, use, 
date of construction, and cost of construction. In some cases, water supply 
reservoirs, alternate water sources and associated costs of development are 
determined. 

The second step is to determine the capacity loss and from this the dam- 
age caused by sediment. 

The damage computations may be based on either of the following: 


1. Damage to present structures. 
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2. Cost of restoration of capacity (dredging). 
3. Estimated cost of alternate source development. 


In the case of water supply reservoirs, the city or industrial engineers in 
charge of the facility can usually supply all of these answers except the rate 
of sedimentation. 

Sediment damage to power reservoirs is complicated. Not only the amount 
of sediment, but the location of this sediment in the reservoir area is impor- 
tant. Therefore, quite detailed studies, in cooperation with the power com- 
pany, are usually necessary. 

A great many power reservoirs have lost all reserve storage and are now 
operating on run-of-the river flow. Sand and other coarse material may then 
move across the silted reservoir, through the intake canals, and into the 
turbines and other equipment. 

This type of damage is evaluated as the cost of replacement parts, labor, 
loss of power and other costs directly connected with repairing the damage 
caused by the sediment. On the other hand, some power companies find it 
cheaper to prevent the damage rather than to repair it. They usually develop 
some method of dredging the sand from the intake canal and prevent its 


reaching the turbines. In this case the annual cost of dredging represents 
the damage. 


Damage to Transportation Facilities 


Sediment deposits damage highways and railways by collecting in ditches 
and culverts and upon the roadways, and by filling and constricting channels 
beneath bridges. These sediment deposits must be removed to keep the high- 
way and railways functioning properly. 

In surveying this type of damage the cost of sediment removal is usually 
obtained from the maintenance division of the highway department. In addi- 
tion to this cost of sediment removal, associated damages are quite high in 
some localities. By causing a malfunction in the drainage system of the 
roadway, damage from frost action is increased enormously. The stability 
of fills and shoulders is often destroyed and the diverted runoff may remove 
fill, road gravel and, in some cases, entire sections of highway surface. 

Highways are usually classified as to type of construction; sometimes the 
fact that they are state, county or township roads indicates the general type 
of construction and expedites the analysis of highway department records. 
In most areas the damage is obtained from the records of the highway de- 
partments. The data consist of the length of road and price level as well as: 


1. Annual cost of removing sediment from ditches, culverts, and road 
surfaces. 


. Annual cost of channel dredging to protect bridges and highways. 
. Annual cost of repairing washouts due to plugged culverts and ditches. 


Field and laboratory studies are sometimes necessary to identify the type 
and origin of this sediment. For instance, winter sanding is a particularly 
important distinction to make. As a general rule, visual inspection of a num- 
ber of damage points will enable the geologist to determine the source of 
damaging sediments. 

In many watersheds, for instance, road ditch sediments are fine materials 
originating from sheet erosion; while plugging the culverts will be bed 
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material of purely local origin from the streambed adjacent to the culvert. 
This type of breakdown is very important in developing a program to allevi- 
ate the damage or in predicting the benefits expected to accrue to the overall 
watershed plan. 

In some watersheds, harbor or canal dredging is necessary to provide or 
protect transportation facilities. In these cases, the annual cost of dredging 
is considered to be the annual damage, but those costs incurred due to water- 


shed sediments must be segregated from costs due to industrial waste, lit- 
toral drift, etc. 


Drainage Ditch Sedimentations 


Drainage ditches, irrigation ditches, and floodways are vulnerable to 
sedimentation. As they fill with sediment they lose their effectiveness. To 
determine the damage, usually it is only necessary to determine the costs of 
clean-out over the years and convert this cost to an average annual cost on 
the price level used for the report. 


Increased Flood Stages 


When channels become clogged by sediment, flood crests for the same 
discharges are constantly forced to higher levels and floodwater damage in- 
creases. 

Since river gages are located in stable reaches, this occurrence is not 
usually discernible from the gaging records of a stream. The reliable 
method of determining this aggradation is by comparing new surveys with 
previous cross section or channel profiles. Unfortunately, very few of 
these older surveys are available and in the case of old highway or railway 
surveys are open to question due to construction activities. Where such 


records are available, the damage due to the increased stage or frequency 
is computed. 


Sediment Damage to Property 


After most floods, deposits of sediment are found on streets and in homes, 
factories, sewers, wells, and other places where they cause damage. The 
cost of removing this sediment is considered the damage. 

It is necessary to obtain records for a period of years for cities, indus- 
tries and farms. These figures can then be reduced to an average annual 
damage during this period. In this regard, the problem of cleaning up after 
a flood, regardless of the amount of sediment deposited, is always present. 
The costs identified as sediment damages are the cost associated with the 
physical effort to remove a mass of material, as distinguished from washing 
down and scrubbing up always required after a flood. We call the “shoveling 
out” a sediment damage and the “scrubbing down” a floodwater damage. 


Determination of Sediment Sources 


Since the Department’s watershed projects encompass sediment damage as 
well as floodwater damage, it is encumbent upon the people developing the 
program to consider sediment sources, sediment entrainment and sediment 
transport. These parameters of the sediment problem are usually expressed 
as the sediment yield at the damage point in question. 
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Some or all of the following types of erosion may be important: sheet ero- 
sion, shallow gullies, deep gullies, valley trenches, streambank erosion, 
floodplain scour, roadside erosion, gravel pits, mine waste, industrial waste 
and construction activities. 

There are several different methods of computing the amount of sediment 
produced by sheet erosion. They were primarily developed to determine the 
amount of soil lost from the field (erosion) and consequently they do not ex- 
press sediment yield without due consideration being given to the processes 
of entrainment and transport. 

Measurements of soil loss from experimental plots on different soils, 
slopes, and with different climatic conditions and rainfall patterns under 
different kinds and combinations of land use were summarized by the Soil 
Conservation Service.3 

Measurements were available from nineteen research stations and 
covered periods of from 5 to 15 years. 

From the study and evaluation of these data, it was determined that the 
rate of sheet erosion is related to a numoer of major factors, and that cer- 
tain relationships exist between these factors. An empirical formula was 
developed: 


E =- Fx g§1.35, 10.36, pl.75, ¢ 
where E is the probable soil loss in tons per acre per year. 

F is a soil factor, based upon the erodibility of the soil and 
other physical factors. 

S is the steepness of slope in percent. 

P is the rainfall. The amount used is the maximum 30-minute 
rainfall expected in the locality from a 2-year frequency 
basis.4 

C is a cropping factor which may be the product of several 
factors related to the use of the land. 


These formulas have been found usable to estimate the amount of sheet 
erosion. 

The amount of sediment originating from gully erosion is determined by 
field surveys, the procedure being similar to the methods used for measur - 
ing streambank erosion. 

Sediment production rates for mine wastes, construction activities, roads, 
gravel pits and related sources are problems commonly encountered but 
with no common solution. Many roads and skid trails can be cross sectioned 
and the volume of erosion to date determined. If the history of that stretch 
of road is known and the extent of the cut involved in construction is deter- 
mined, the gross amount of erosion, to date, can be determined. Average 
annual amounts can then be estimated. 

Areas involved in housing developments, super highways, airports, strip 
mining and other activities are given independent consideration. The amount 
of sediment produced is usually estimated by the geologist, in consultation 
with the others on the planning party, as well as the engineers in charge of 
the development. Many of these sources are, here-today-gone-tomorrow 


3. G. W. Musgrave - The Quantitative Evaluation of Factors in Water Ero- 
sion, Journal of Soil and Water Conservation, July 1949. 
4. See U.S.D.A. Misc. Pub. No. 204. 
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2} | affairs, that have little influence upon a long range watershed program. On 
ae the other hand, many areas are being converted from agricultural to urban 
b use and the effect of this conversion upon sediment production, as well as 


4 hydraulic characteristics of the watershed, must be recognized. 
Sediment Source - Damage Relationship 


The total sediment yield from the various types of erosion, on an annual 
basis, are converted to percentages of total sediment. The relative impor- 
tance of each type of erosion can then be determined for specific damages. 

These figures apply only to the sediment as a whole within the watershed. 
The relative importance of the various sediment sources in the case of 
particular types of sediment damage will vary either way from these mean 
values. Some types of sediment damage are caused largely by one type of 
sediment. The location of the damage point, whether in the headwaters or 
at the mouth of the river, influence the relative importance of its sediment 
source. Some types of sediment, after erosion, are redeposited a short dis- 
tance from their original source. Some types of erosion, such as floodplain 
scour, may not occur in headwater areas. It is therefore necessary to de- 
velop for each type of sediment damage and for each damage area, the rela- 
tive importance of the various sediment sources. The measured relative 
source areas are used as a guide. These data are most usable when con- 
verted to percentages of damaging sediment and tabulated in the following 
form: 


SEDIMENT SOURCES IN PERCENT "X" WATERSHED 


- Flood- Road- Av. Annual 
sheet Gully plair side Dollar 
Erosion Erosic Scour E ion Total Damag 
65 10 LOO 
50 30 10 10¢ 31, LOC 
40 kc 429 SVU 
L U 4C i 100 21.70C 
- - 100 535 
- 100 1,600 
= 42 OC 
Developing the Program 
; a The tabulation of sediment sources by percent is a useful device, not only 
—_ to show the relation of the sources to the damage but also in computing the 


effects of the program. By referring to the above table, the effect of any 
erosion reduction in terms of benefit are evident. 

For example, 50 percent of the overwash damage, amounting to $15,550 
per year, is due to sheet erosion. Assuming that sheet erosion can be re- 
duced by 60 percent, the benefit from reduction of overwash damage would 
be $9,330 or $31,100 x .50 x .60 = $9,330. 

In developing a watershed program, a great many factors affect the final 
selection of the practices and measures to be recommended. In general, the 
items to be included in the plan are those items which offer the most benefit 
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with the least cost. It is quite obvious that the results of the damage surveys 
and the effects of the independent and interdependent measures that are pro- 
posed for inclusion in the plan must be presented in such a manner as to 
facilitate alternate appraisals. 

The effect of the land use changes can be predicted through the recomputa- 
tion of soil losses - with the land treatment measures installed, using the 
soil loss formula previously discussed. 

The effect of structural measures on sediment production can also be pre- 
dicted. 


The sum of these effects provides a basis for predicting the reduction of 
sediment damages. 


Sedimentation Considerations in Detail Designs 


Following the acceptance of the watershed plan by the local people and its 
approval by the Government, detailed designs for all of the works of improve- 
ment are made. 

The sediment storage requirements in floodwater detention reservoirs 
are computed for each individual reservoir and the data are usually summa- 
rized on a form similar to Figure 2. 

The estimated sediment yield used in the design of a particular floodwater 
retarding structure is based on watershed conditions above the structure. 
These computations, for a 50-year period, are based upon conditions at the 
time of the inspection, adjusted for the amount and effectiveness of land 
treatment and other measures to be applied in the watershed. 

The trap efficiency of the reservoir is estimated by the use of envelope 
curves based on capacity -inflow ratio.) The total capacity (sediment pool 
plus retarding pool) to emergency spillway level of the reservoir is used in 
determining the capacity-watershed or capacity-inflow ratio. 

The weight of the sediment in the reservoir is estimated from volume- 
weight determinations made on sediments in farm ponds or other reservoirs 
within the area. In computing the volume of sediment in dry pools and in the 
floodwater retarding area, ultimate compaction of the aeriated deposit is as- 
sumed. Compaction of volume of sediments below the normal pool level is 
based on ultimate compaction under submerged conditions. 

The distribution of sediment in the reservoir is dependent on the nature of 
the sediment, the topography and shape of the reservoir, nature of the sedi- 
ment, the topography and shape of the reservoir, nature of approach channel 
above it, and the time required to empty the retarding pool. Allocation of 
capacity for sediment accumulation in the reservoir is made in two general 
locations, the lower of which is termed the sediment pool (below the normal 
pool level) and the higher of which is termed the retarding pool. The volume 
of sediment estimated to occur in the sediment pool may fix the elevation of 
the crest of the principal spillway. 

The problem of bed load transport, both as related to floodwater detention 
structures and other installations deserve mention. Where total sediment 
transport rates are high, bed material accounts for a minor part of the total 


load and, while in no case is it ignored, it is realistically given minor con- 
sideration. 


5. Trap Efficiency of Reservoirs, G. M. Brune, Transactions A.G.U., Vol. 34, 
No. 3, 1953. 
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In the watersheds, such as the mountainous forest and pasture areas of 
the Northeast, where total transport is low, bed material transport may ac- 
count for most of the sediment damage and consequently warrants close con- 
sideration. A number of formulas of the Du Boys type have been, and are 
being used, as well as the procedure developed by Dr. Hans Einstein.6 

Unfortunately, these calculations, while requiring a rather high order of 
precise mathematics, give an index of magnitude rather than a precise esti- 
mate. This is not a deterrent in planning activities since most of the dam- 
ages caused by bed material (such as channel filling, bridge and culvert 
plugging, etc.) involve a physical effort for their relief. These damages 
are usually expressed in the cost of removing the material. Consequently, 
the benefit is expressed as a percent reduction in damages. 

Design problems are more difficult since they often involve a volumetric 
estimate. 

Since none of the existing bed load formulas are precise when applied to 
untested conditions, the procedure most generally followed is to compute the 
bed load transport by two or three formulas,? and from these computations 


select, by judgment, the most conservative result to use in design of the 
work of improvement. 


. The Bed Load Function for Sediment Transportation in Open Channel 


Flows, H. A. Einstein, U.S.D.A. Technical Bulletin No. 1026, September 
1950. 


. Engineering Hydraulics, page 795, edited by Hunter Rouse, 1949. 
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Proceedings of the American Society of Civil Engineers 


STILLING BASIN EXPERIENCES OF THE CORPS OF ENGINEERS# 


R. H. Berryhill,! A.M. ASCE 
(Proc. Paper 1264) 


SYNOPSIS 


The Corps of Engineers has designed and built more than 150 reservoir 4 
projects during the past twenty-five years. Most of these projects have been ? : 
built since the Congress passed the Flood Control Act of 1936. Some of the = 
experiences of the Corps of Engineers in the design, operation and mainte- - 


nance of the stilling basins below the spillways and outlet works of these pro- 
jects are described herein. Aliso described are several model and prototype 
correlations and some conclusions which may be drawn from the experiences 
described. This paper is one of a series of four papers assigned to a task 
force of the Design Committee of the Hydraulics Division on the general sub- 


ject “Energy Dissipators for Spillways and Outlet Works.” Soe 

a 

4 

INTRODUCTION 
A number of different types of stilling devices have been used for the safe Zl 
handling of flows released from the many spillways and outlet works built by -_ 


the Corps of Engineers. For spillways, the stilling devices used include 
conventional horizontal stilling basins with and without baffle piers, end sills 
and end sill controls, as well as sloping and stepped stilling basins, two-level 
basins, and bucket-type basins. Under the classification of bucket-type still- 
ing basins, the following types have been used: roller buckets, deflector 
buckets, and high deflector buckets with or without subdams downstream. 
Many navigation and power dams under the jurisdiction of the Corps of 
Engineers, with “run-of-the-river” flows, have experienced major spillway - 


Note: Discussion open until November 1, 1957. Paper 1264 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
83, No. HY 3, June, 1957. 

a. Presented at American Society of Civil Engineers Annual Meeting, 
Pittsburgh, Pa., October 18, 1956. a a 

1. Civ. Engr., Southwestern Div., Corps of Engrs., U. S. Dept. of the Army, wh 
Dallas, Texas. a 
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releases which have provided good tests of the adequacy of the stilling basin 
design. Notable among these is the Bonneville Dam spillway, which passed 
1,000,000 cu. ft. per sec. during the flood of 1948. Of the reservoir projects 
with flood control storage, however (and these comprise the majority) there 
are few spillways which have experienced large flow releases. This is of 
course brought about primarily by the fact that the flood control storage pro- 
vided will often contain the entire flood of record on this type of project, and 
only rare floods will produce major spillway discharges. Since the average 
Corps of Engineers reservoir project has been in operation less than 15 
years, a number of years may elapse before many of these projects experi- 
ence large spillway flows. 

Discharges from the outlet works at most Corps of Engineers reservoir 
projects are of course frequent. At several projects where the stilling basin 
is common to both spillway and outlet works, flows released from the outlet 
works have been more critical for design of the stilling basin than have the 
flows from the spillway. Energy dissipating devices that have been used for 
outlet works include most of the types used for spillways as enumerated 
above, and in addition include deflector or splitter blocks. 

Discussed in the following pages are some of the stilling basins and other 
stilling devices which may be of interest either because of certain flow ex- 
periences, model-prototype comparisons, or difficulties encountered during 
operation. For grouping purposes, the first section involves a discussion of 

' stilling devices below spillways where the discharges from the spillways are 
of primary interest or importance. The second section deals with stilling de- 
vices below spillways where the discharges from sluices located within the 
spillway weir section are of primary interest. Finally, stilling devices below 
certain tunnels and conduits are discussed, followed by conclusions which 
may be drawn from some of the experiences described. 


Spillway Stilling Devices 


McNary Dam (Washington and Oregon) 


At this navigation and power project, the floow flows of the Columbia River 
are passed through a 1,310-ft.-long concrete gravity ogee spillway containing 
twenty-two 50-ft. by 50-ft. gates separated by 10-ft.-wide spillway piers. 

The stilling basin is the conventional hydraulic jump type and is designed to 
pass 2,200,000 cu. ft. per sec. at pool El. 356.5. As shown in Fig. 1, the ap- 
ron is set at El. 228.0 and is 270 ft. long, with two rows of baffle piers and 
an end sill. The structure is unusual in that the maximum head on the spill- 
way crest is 65.5 ft., resulting in a discharge into the stilling basin of 1,680 
cu. ft. per sec. per ft. of basin width, or 100,000 cu. ft. per sec. for each gate 
bay. 

Extensive use was made of hydraulic models in the development of this 
project. The spillway models included a 1:100 scale general model, a 1:36 
scale sectional spillway model and a 1:10 scale single gate bay model. The 
apron length selected was equivalent to 3.2D9, where Dg is the maximum 
conjugate depth after the hydraulic jump. The apron elevation was set so that 
the maximum tailwater depth of 75.5 ft. was equal to 0.90D2 because of the 
anticipated effect of the baffle piers. It should be noted that the Froude Num- 
ber (V1/ gD1) at design discharge is only 3.3, where V1 and Dj are the 
velocity and depth, respectively, upstream from the hydraulic jump. 
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Fig. 1. Stilling Basin Dimensions, McNary Dam. 
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The reservoir pool was raised in December 1953, upon completion of most 
of the project features. Prototype operation tests were made in June and July 
1955 with normal pool elevations of 335 and 340 and with spillway Bays 15, 16 
and 17, numbering from the north bank, discharging together and with Bay 16 
discharging alone. The purposes of the tests were to determine the magni- 
tude of the pressures on the spillway crest, gate piers and baffle piers by 
means of piezometers installed in Bay 16, and to determine the conformity 
between performance of the spillway models and the prototype. The pres- 
sures obtained for the crest, gate slots and piers will not be discussed here, 
except to state that very good correlation with the models was obtained in 
nearly all instances. 

The baffle pier pressures obtained from the prototype tests are shown in 
Fig. 2. It is noted that piezometer B-1 is located in the middle of the up- 
stream face and B-2 is located at the downstream end of the one-ft. radius 
along the south side of the baffle. The opposite side of the baffle coincided 
with the centerline of Bay 16. A third piezometer, located in the top of the 
baffle, was apparently clogged during the test period. The tailwater elevation 
was approximately 269.5 for all tests made with reservoir pool at elevation 
335. With pool elevation 340, tailwater elevation was approximately 268.5 for 
all conditions except for bays 15, 16 and 17 operating together with gate open- 
ings of 1, 3, 5 and 10 feet, when tailwater elevation was 264. 

The low pressures shown for piezometer B-2 indicate that some pitting and 
erosion of the sides of the baffles may be expected in the future. However, 
the flow conditions during the test period are not considered normal. Since 
normal spillway operation at McNary Dam is accomplished by nearly uniform 
gate settings for all the bays, the larger gate settings would only be used dur- 
ing periods of large river flows and higher tailwater elevations. For example, 
full gate openings would normally be accompanied by a tailwater level about 
20 to 21 ft. higher than that used in the tests. The operation of a single gate 
at a large opening will occasionally be required for short periods of time 
when a load rejection at the powerhouse necessitates an immediate release 
from the spillway. Based on the above considerations, it is apparent that the 
pressure conditions shown in Fig. 2 are not as critical as the data seem to 
indicate. 

No comparable information on baffle pier pressures in the models was 
available. However, pressures near the cavitation range were noted in the 
model for large flows of short duration. 

An indication of the hydraulic losses incurred from reservoir pool to the 
upstream row of baffles is given by the graph for piezometer B-1 in Fig. 2. 
Total losses with full gate openings varied from 6.9 ft. for pool El. 340 to 
10.9 ft. for pool El. 335, while as percentages of the total available head the 
losses were 6.5 and 10.7, respectively. It is hoped that additional test read- 
ings in the future will either verify these values or narrow the range of 
results. 

The observed stilling basin action was satisfactory, although as might be 
expected, the large concentrations of flow in the test bay created extremely 
rough and turbulent flow conditions. With pool elevation 340 and free flow in 
Bay 16, the hydraulic jump remained in the stilling basin with tailwater 16 ft. 
(28%) lower than the minimum tailwater required for a jump as determined in 

the 1:36-scale sectional model. Apparently the side inflow which could not be 


duplicated in the sectional model was sufficient to partially stabilize the jump 
in the prototype test bay. 
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Fig. 2. Baffle Pier Pressures, McNary Dam. 
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Fig. 3 shows the turbulent jump action which occurred during passage of 
the June 1956 flood with flows of 735,000 to 787,000 cu. ft. per sec. As nearly 
as can be determined by the author from comparison of photographs, the pro- 
totype stilling action was as predicted by the models. 

It should be noted that the basin length at McNary Dam was made some- 
what shorter and the maximum end sill velocities allowed were somewhat 
higher than would be considered desirable at some locations, because of the 
very hard durable basaltic rock which lies below the McNary end sill. 

It is interesting to note that the downstream portion of the stilling basin 
training wall is only 7 ft. high beyond a point 105 ft. downstream from the P. L. 
of the spillway bucket curve. The wall was made low in order to provide bet- 
ter escape conditions at certain flows for the Columbia River Salmon migrat- 
ing upstream. No harmful effects were noted as a result of this feature in the 
general model or in the prototype with the 1956 flood flows. Possibly the low 
Froude Number accounts for this performance. 


Bonneville Dam (Washington and Oregon) 


This navigation and power project was completed in 1937. During the ensu- 
ing twenty years, the large daily flows and annual flood flows of the Columbia 
River have subjected the stilling basin baffle piers to extremely heavy pound- 
ing for prolonged periods of time, resulting in gradual erosion of the concrete 
over large areas of the baffle piers and floor. It is estimated that about 
1,000,000,000 acre-ft. of water have passed over the spillway to date (1956). 
Details of the erosion and the repairs made have been described else- 
where.l, 2. However, several features are of particular interest and are 
mentioned briefly here. 

Nine baffle piers were repaired in 1947 with 9-in. radius on each leading 
edge and a 2.5-ft. radius connecting top and back slope. This design was de- 
veloped as a result of model tests conducted in 1941 in a vacuum tank at 
Carnegie Institute of Technology. During the period 1947 to 1954, these 
“Carnegie” baffle piers sustained considerably less general wear than the 
original type of baffle pier without rounded edges, as shown in Fig. 4. In the 
first of these baffle piers to be reshaped in 1947, deep erosion was later 
found in areas where unsatisfactory bond had been obtained between the new 
and old concrete. The change in design had been effected by removing the 
outer 6 in. of concrete on all exposed faces. During the recent (1954) repairs, 
a minimum of 12 in. thick concrete was provided and anchored with No. 8 
hooked reinforcing bars grouted 2 ft. into the original concrete. 

One baffle pier in the upstream row was repaired in 1947 with a 3/8-in.- 
thick steel armor plate all around, anchored with bolts spaced on 12 in. cen- 
ters. All of the armor plate had been ripped off and the sides were deeply 
pitted when inspected in 1954. 

During the great Columbia River flood of 1948, the stilling basin was sub- 
jected to a severe test. For a one-day duration, the peak river flow was 
1,000,000 cu. ft. per sec., whereas the design discharge is 1,600,000 cu. ft. 
per sec. Stilling action was very turbulent, as is shown in Fig. 5. Eye 


1. Symposium on Cavitation, “Experiences of the Corps of Engineers,” by 
John C. Harrold, Transactions, A.S.C.E., Vol. 112, 1947, p. 35. 

2. “Bonneville Dam Spillway Erosion Repaired,” by Donald H. Basgen, 
Engineering News Record, Apr. 21, 1955. 
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Fig. 3. McNary Dam Stilling Basin During 1956 Flood Flows, 
(a) 735,000 cu. ft. per sec., (b) 787,000 cu. ft. per sec. 
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witnesses reported large, almost instantaneous surges near the spillway 
abutments in the order of 25 to 30 feet vertically. The total flow remained 
over 900,000 cu. ft. per sec. for 18 days duration. 

Recent observations of flows in the south half of the stilling basin, where 
the downstream row of baffle piers has been filled in solid to form a sill, in- 
dicate more turbulent surface conditions than the original basin. Model 
studies had indicated that the new design basin with solid end sill on the baf- 
fle deck would deflect the flow upward, producing more turbulence on the sur- 
face but providing a large, slow-moving eddy in an upstream direction on the 
bottom of the river bed downstream from the lower apron. 


Center Hill Dam (Tennessee) 


The spillway at this project is a high concrete gravity ogee weir with 
roller bucket type stilling basin. The basin is 470 ft. wide and is divided into 
three levels of 7-ft. increments, vertically. Radius of the bucket is 50 ft. and 
maximum head available from pool to bucket invert is 233 ft. At the lip of the 
bucket, the angle of tangent with the horizontal is 45 degrees. Design dis- 
charge is 458,000 cu. ft. per sec. and the unit discharge is approximately 973 
cu. ft. per sec. per ft. width of basin. 

Construction work was suspended from March 1943 to January 1946 be- 
cause of the war effort, and during this time the flow of the Caney Fork River 
was diverted through the low spillway monoliths and the sluices. A large 
quantity of rock was spilled in the bucket during construction of a rock fill 
cofferdam adjacent to the spillway. This rock was not removed and the vio- 
lent eddies set up by sluice flows during the diversion period resulted in the 
formation of several large oval-shaped cavities in the bucket. The abrasive 
action apparently continued over a period of eight years until 1951, when the 
damage was discovered. The maximum spillway discharge during this period 
was 41,000 cu. ft. per sec. At this flow, waves about 3 ft. high were caused 
by the jets from the sluices breaking through the tailwater. A release of 
40,000 cu. ft. per sec. was maintained for 8 days. Normal spillway releases 
since completion of the dam occur about 10 days a year at a rate of about 
5,000 cu. ft. per sec. 

No difficulty has been experienced with spray action from the bucket dur- 
ing operation of the powerhouse located to the left of the spillway or the 
switchyard located just downstream, although the stilling action of the bucket 
is somewhat more violent than a conventional hydraulic jump type basin. The 
greatest amount of spray occurred with sluice discharges during the diversion 
period. Since the powerhouse was placed in operation, there have been no re- 
leases through the sluices. 

It was found that flow releases should be evenly distributed across the 
spillway, since the back roller that is formed immediately below the bucket 
will tend to deposit fill against the downstream lip of the bucket. With one 
section of the lip left low, unbalanced releases may easily cause gravel to be 
deposited in the bucket. A level bucket would be preferable in this regard. 

The Wolf Creek and Dale Hollow projects, located nearby, have similar 
spillways and roller bucket basins. Since Wolf Creek Dam was completed in 
1950, the maximum flow release has been about 40,000 cu. ft. per sec. Con- 
tinuous releases have been made for 60 days averaging about 22,000 cu. ft. 
per sec. Since closure of the Dale Hollow Dam in 1943, the largest release 
has been only 10,200 cu. ft. per sec. No erosion has been apparent at either 
of these projects. 
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Baldhill Dam (North Dakota) 


At this project, the 140-ft. long ogee spillway section is joined to a chute 
of 15% slope, slightly flaring in plan, and a level apron 70 ft. long and with 
flare in plan of 1 in 6, averaging about 155 ft. wide, complete with two rows 
of baffle piers and an end sill. Maximum head available from pool to apron 
is 51 ft., and the design discharge is 43,100 cu. ft. per sec. 

A flood of record proportions occurred in April 1950, before completion of 
construction, which required placing the spillway in emergency operation. 
Since much more than normal storage space was available, the maximum re- 
lease was maintained at 3,150 cu. ft. per sec. Comparisons between model 
and prototype at this flow indicated good agreement. In both, the hydraulic 
jump began and ended on the chute section, and the flow was distributed fairly 
uniformly across the width of the basin. 

For discharges above 10,000 cu. ft. per sec., the model indicated that the 
jump would begin on the chute and end on the level stilling basin floor. Down- 
stream depths for this type of flow in the model were fairly consistent with 
those obtained from empirical relationships developed by Carl E. Kindsvater, 
M. ASCES3 for flow in a straight chute and basin. Larger prototype flows than 
the 1950 releases would be desirable for comparison, but cannot readily be 
obtained because such flows would cause flooding conditions in communities 
downstream from the dam. 


Conchas Dam (New Mexico) 


This project was completed in 1939. The spillway consists of a concrete 
gravity ogee weir 340 ft. wide and 178 ft. high from crest to toe. The basin 
consists of a horizontal apron 134 ft. long with one row of 8-ft. high baffle 
piers and a 12-ft. high end sill. The basin length is unusually short for a 
project of this size, and is equivalent to about 2.3D2, where Dg is the maxi- 
mum conjugate depth after the hydraulic jump. In addition, the ends of the 
training walls are tapered back 1.1:1 from the end sill. Spillway design dis- 
charge is 182,000 cu. ft. per sec. Model tests were made of the hydraulic 
jump action prior to construction. 

In 1944, the basin was dewatered for inspection. Total releases up to that 
time had amounted to two million acre-ft., of which 700,000 acre-ft. were 
from spillway releases during 1941 and 1942, and the remainder from sluice 
releases. Inspection revealed that most of the basin was in good condition. 
Behind the baffle pier attached to the training wall on either side, however, 
the floor was eroded to a depth of 9 to 12 inches in the area bounded by the 
baffle pier, training wall and end sill. The wear also extended into the train- 
ing wall and downstream face of baffle pier as much as 6 in. This wear was 
caused by eddies which formed in the “dead” spots behind the two end baffle 
piers. In later basins, the end baffle piers have been placed a short distance 
away from the training walls in order to prevent such eddy action. In addi- 
tion to the above, erosion on the level apron downstream from two of the 
regulating sluices varied from a depth of 3 in. near the baffle piers to virtual- 
ly none at the spillway toe. Repairs consisted of replacing the concrete be- 
hind each of the end baffles. 

Downstream from the end sill, the exit channel side slopes are protected 


3. “The Hydraulic Jump in Sloping Channels,” by Carl E. Kindsvater, 
Transactions, A.S.C.E., Vol. 109, 1944, pp. 1107-1154. 
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by a 10-ft. blanket of sand and gravel covered with a 3-ft. layer of derrick 
stone on slopes of 1 1/2:1. Substantial displacement of the derrick stone oc- 
curred after the 1941 and 1942 flood releases and subsequently. Causes of 
the failure were thought to be slippage of the over-burden along the surface 

of the shale and sandstone, influenced by saturation of the overburden, under- 
cutting of the revetment toe by the violent eddies from the stilling basin, and 
vibration of the revetment which had been noticed earlier during spillway 
releases. 

Observation of the prototype spillway flows up to 38,000 cu. ft. per sec. 
indicated that the hydraulic jump was completely submerged, whereas model 
tests had indicated that a free jump should occur for flows less than 75,000 
cu. ft. per sec. for all anticipated tailwater conditions. Since the prototype 
tailwater elevations were within the expected limits, although near the maxi- 
mum anticipated, the reasons for this discrepancy between prototype and 
model flow action are not known. 

As noted elsewhere ,4 degradation of the streambed downstream from this 
project is recognizable for a distance of more than 40 miles. The streambed 
near the dam has lowered a maximum of 10 ft. to rock outcrops. 


Folsom Dam (California) 


The spillway for this project is a high concrete gravity ogee weir topped 
with five 42-ft. wide by 50-ft. high “operating” tainter gates and three similar 
“emergency” tainter gates, 42-ft. wide by 55-ft. high. The service stilling 
basin below the five operating gates consists of an 8:1 ramp, 175 ft. long 
joined to a level apron, 147 ft. long. A 15-ft.-high vertical end sill is also 
provided. The height of the weir crest above the level apron is 303 ft. Below 
the emergency gates, a high level bucket is provided. Design discharge for 
the service stilling basin is 200,000 cu. ft. per sec., although all structures 
were designed to be safe against failure during the spillway design flood of 
567,000 cu. ft. per sec. Since the service basin is 242 ft. wide, the design 
unit discharge is 826 cu. ft. per sec. per ft. of width. The stilling basin de- 
signs were evolved in conjunction with a series of model tests. A level still- 
ing basin was chosen for the service basin instead of a flip-bucket because of 
increased safety to the adjoining powerhouse and because maintenance costs 
would be reduced, with less eroded material entering the tailrace area. 

During the flood of December 1955, a flow of 64,000 cu. ft. per sec. was 
released from the five service gates, providing a unit discharge of 264 cu. ft. 
per sec. per ft. of width. The pool level was 36.7 ft. above the spillway crest 
and 340 ft. above the level apron. Tailwater was 63 ft. above the level apron. 
As predicted by the model, stilling action was less than complete, and moder- 
ately high velocity flow continued downstream from the end sill for a consider- 
able distance, as shown in Fig. 6. In the model, a simulated discharge of 
110,000 cu. ft. per sec. resulted in bottom velocities of 16 ft. per sec. imme- 
diately downstream from the basin end sill. However, with a scour hole simu- 
lated in the exit channel in this area, the model indicated bottom velocities of 
only 12 ft. per sec. and upstream in direction. 


4. “Sediment Studies at Conchas Reservoir in New Mexico.” by D. C. 
Bondurant, Transactions, A.S.C.E., Vol. 116, 1951, p. 1288. 
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Arkabutla Dam (Mississippi) 


This project has a 300-ft. long uncontrolled concrete chute spillway with 
ogee weir and 3:1 chute slope. Design discharge is 89,000 cu. ft. per sec. 
The stilling basin as developed from model tests in 1940 is 95 ft. long, with 
level apron, two rows of baffle piers, and an end sill. In 1953, the spillway 
went into operation for the first time, with a maximum discharge of 7,100 cu. 
ft. per sec. Measurements were made of the spillway discharge and water 
surface elevations along the chute. Stilling action was observed to be satis- 
factory and well-distributed across the basin. 

It was noted that tailwater for the 1953 flow was 16 ft. higher than that 
used for the model tests. After inquiring, it was found that the outlet channel 
was not constructed to the bottom grade and width as used in the model study. 
The channel had been excavated to a bottom grade varying from 5 ft. higher 
at the lower end to 13 ft. higher at the upper end than the bottom grade used 
for the determination of the tailwater data for the model tests. Also, the 
spillway exit channel was densely grown up with vegetation, resulting ina 
roughness coefficient “n” in Manning’s formula of 0.070 during the 1953 flood. 


Lock and Dam No. 1, Mississippi River (Minnesota) 


The ogee, Ambursen type fixed crest spillway at this navigation project is 
574 ft. long and has a level concrete apron approximately 80 ft. long at an 
elevation 30 ft. below the spillway crest. The concrete apron is too high for a 
hydraulic jump except for higher flows. In 1952, the maximum flood of record 
occurred and a discharge of 69,200 cu. ft. per sec. passed over the spillway. 
The hydraulic jump apparently formed near the downstream edge of the apron. 
Tailwater depth at this time was 21 ft. Excessive scour occurred below the 
apron, and a hydraulic model study was made to determine the most feasible 
method of modifying the apron. As a result, a baffle wall was constructed 
with 2-ft. by 2-ft. openings which would move the hydraulic jump well up on 
the apron. The baffle wall was placed about 30 ft. downstream from the toe 
of the weir. Action of the basin with a discharge of 21,500 cu. ft. per sec. is 
shown in Fig. 7. 

There are 26 navigation dams on the Upper Mississippi River. At many of 
these projects, the spillway structures are founded on pile structures in the 
sand bed of the river. In the usual case, the low concrete sill is followed by 
a concrete stilling basin about 50 ft. long containing two rows of baffles, an 
end sill and approximately 50 ft. of derrick stone below the end sill. In a few 
instances, this stone has been undermined by erosion of sand and gravel lo- 
cated further downstream. The usual cause of this loss has been concentrated 
flow when some of the spillway gates were bulkheaded for extended period of 
time. Remedial works usually require that a timber mattress be constructed, 
floated into position, sunk, and the entire mattress covered with a one-ft. 
layer of one-man stone. This work has been relatively permanent. The gate 
operation schedule on a number of these projects founded on piling is so ar- 
ranged that the maximum velocity over the end sill does not exceed about 5.5 
ft. per sec., as determined from model tests with erodible beds. 


Demopolis Lock and Dam (Alabama) 


This navigation project was completed in November 1955. A flood occurred 
in April 1955, when the dam was essentially complete except for several 
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cofferdam cells which were still in place. The maximum discharge was 
85,000 cu. ft. per sec., with pool at elevation 79.6 and tailwater at 74.0, or a 
net head of 5.6 ft. In the model, the most critical stilling basin conditions 
were obtained with a flow of 80,000 cu. ft. per sec. and a net head of 5.9 ft. 
The observed performance of the prototype appeared the same as the model. 
No erosion was reported. 

The stilling basin is a roller type, with two sections. One section is 585 ft. 
long and has a 27.5 ft. level apron set at El. 19.0, surmounted by a 4.4 ft.- 
high sloping end sill. The other section is 895 ft. long and has a 29-ft. level 
apron set at El. 45.0, with a 3-ft.-high sloping end sill. Spillway crest is at 
El. 73.0. 


Old Hickory Dam (Tennessee) 


This is a low-head power and navigation project with normal lift of 60 ft. 
A bucket-type energy dissipator, with bucket radius of 34 ft., was provided 
in order to deflect low to medium spillway flows upward and away from the 
streambed. At high flows, the spillway will be submerged and the overflow 
will ride the upper layers of the tailwater. This project is still under con- 
struction. However, diversion flows up to 34,000 cu. ft. per sec. have been 
passed through 4 of the 45 ft. wide gate bays, creating more turbulent stilling 
action than will be experienced after the project is completed, since tailwater 
was appreciably lower than normal during this period. Stilling action ap- 
peared turbulent but satisfactory. 


Spillway Stilling Devices with Sluice Flows 


Bull Shoals Dam (Arkansas) 


This dam is a straight, concrete gravity structure 284 ft. high, with an 
overflow spillway located in the central portion of the valley. Flow over the 
spillway is controlled by 17 tainter gates, 40-ft. wide by 28-ft. high. Release 
of floodwaters when the reservoir pool is below spillway crest is accom- 
plished by 16 sluices, 4-ft. wide by 9-ft. high, controlled by hydraulic slide 
gates. As shown in Fig. 8(a), the stilling basin was constructed with a 
stepped apron 808-ft. wide and 203.3-ft. long, with elevations varying from 
437.2 to 446.5, and a 4-ft.-high end sill. The apron had a 1.0 ft.-high step 
located 35.3 ft. downstream from thé toe of the dam, a 1.5 ft.-high step lo- 
cated 47.3 ft. downstream from the toe and 2.5 ft.-high steps located 71.8 and 
98.8 ft. downstream from the toe. Each step was shaped as a quadrant of an 
ellipse, with semi-minor axis equal to the height of the step and semi-major 
axis of twice this amount. This design is similar to the basin at Cherokee 
Dam (TVA) and was adopted after model tests were made in 1944-45. The 
alternative design was a more costly basin with horizontal apron at El. 435 
equipped with two rows of baffle piers and end sill. Stilling action for sluice 
flows was a controlling factor in the basin design. The dam was completed 
in 1951. 

After approximately 2,600,000 acre-feet of water had been released 
through the sluices, extensive damage was noticed near several of the sluice 
outlet portals. The basin was unwatered in October 1952 and a careful inves- 
tigation made. Erosion, varying from a depth of less than 2 feet to a maxi- 
mum of 3.5 feet had occurred in the 5-foot thick concrete slab below the outlet 
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portals of twelve of the sixteen sluices. Typical of the four or five most 
severely eroded areas was monolith 22, shown in Fig. 9(a). The severe 
erosion on top of the first step appeared to be the result of cavitation. Once 
started, it seemed probable that the erosion process may have been acceler- 
ated by dislodged concrete and broken reinforcing steel. However, neither 
the surface character of the eroded holes nor the stone and steel indicated 
that abrasion was the primary cause. There was clear evidence of impact 
abrasion or sandblast effect on the faces of the second and third steps down- 
stream. 

An excellent example of the erosion of the first step in its early stages 
was afforded by monolith 35, shown in Fig. 9(b). The sluice in this monolith 
had passed only about 5,000 acre-feet of flow, as compared to 420,000 acre- 
feet for the sluice in monolith 22 described above. 

Temporary repairs were made in the damaged areas while awaiting the 
results of prototype pressure cell tests, new model tests and new analyses 
leading to a permanent solution. The damaged steps were cut back and re- 
placed with smooth, upward sloping slabs. The slabs replaced the two up- 
stream steps at all conduit outlets except two. At the outlet in Monolith 21, 
adjacent to the left training wall, all four steps were replaced. At Monolith 
23, where the sluice had not been in operation, the undamaged steps were re- 
tained for testing purposes. 

Electrical pressure cells were installed in the stilling basin in conjunction 
with the repair work described above. At the top of the first step in Monolith 
23, three cells were installed along the sluice centerline at distances of 
2.67, 3.67 and 4.67 ft., respectively, from the upstream end of the step. 
Reservoir pool El. during the prototype tests was 638.4 and tailwater just 
downstream from the end sill was at El. 452. Conduit discharge was approxi- 
mately 3,300 cu. ft. per sec. at full gate opening and 1,300 cu. ft. per sec. at 
half gate opening, with corresponding velocities of about 92 and 82 ft. per sec., 
respectively. At the half gate opening, instantaneous pressure readings in the 
cavitation range (-35.8 ft. of water) were obtained at the pressure cell fur- 
thest upstream. The average pressure reading for this cell at half gate was 
+3.0 and the maximum was +24.3 ft. of water. At full gate opening, the three 
cells were rendered inoperative. The conclusion derived from the results of 
these tests was that the damage had resulted from cavitation. 

It is interesting to note that the 1944-45 model tests, using piezometers 
connected to a manometer bank, gave a fair prediction of the average proto- 
type pressure (+6.5 ft. versus the +3.0 ft. reading quoted above). Of course 
the piezometers could not reflect the transient pressures encountered. 

Observation of the prototype flow from the conduit in Monolith 23 with 
tailwater completely repelled because of a partial cofferdam indicated that the 
spread of the jet was complete over the width between the flared guide walls 
before it reached the end of the guide walls and was almost complete at the 
location of the first step. At the latter location, the jet was slightly thicker in 
the center than on either side, and compared closely with the cross sections 
and longitudinal profile of the unsubmerged jet as measured in the 1944-45 
model tests. 

Observations of the prototype stilling basin performance with sluice dis- 
charges and with ramps installed in lieu of the first two steps indicated that 
tailwater elevations below approximately elevation 459 were insufficient to 
create a hydraulic jump and consequently the conduit jet extended beyond the 
basin end sill with considerable energy for some distance downstream. The 
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Fig. 9. Erosion at Bull Shoals Dam, (a) Eroded Step of Monolith 22, 
(b) Pitting at Top of First Step of Monolith 35. 
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deficiency in tailwater was verified in the model, and an analytical study also 
indicated that a tailwater elevation of 460 to 464 would be required for an 
adequate hydraulic jump. In order to increase the effective tailwater depth 
for sluice flows, the final design for the stilling basin additions incorporated 
a 12-ft. high end sill, as shown in Fig. 8(b). The upstream face of the end sill 
was provided with a 1:2 slope in order to increase the lateral spread of indi- 
vidual conduit flows over the sill. Since a relatively large drop-down occurs 
over the high end sill, a 24-ft.-long secondary basin with 2-ft.-high sill was 
added. One row of 8-ft. high baffles, 7.5 ft. wide spaced 8.5 ft. apart, was 
added to the main basin to improve stilling action for large spillway flows. 
The elliptical step furthest downstream was left intact, as shown in the figure. 
Repairs were completed early in 1955. 

The concrete mix for the original basin contained manufactured sand and 
about 4 1/2 bags of cement per cu. yd. For the basin revisions, the mix 
utilized natural sand, 5 1/2 bags of cement per cu. yd., 1 1/2 in. maximum 
size aggregate, and a water-cement ratio of 0.41 to 0.45. 

From soundings made of the exit channel in the first 100 feet downstream 
from end sill in 1948 and again in 1953, erosion in the order of 2 to 5 feet was 
found in a number of locations in the moderately hard limestone. 


Nimrod Dam (Arkansas) 


This project was completed in 1942. The spillway is composed of a con- 
crete gravity ogee weir with a conventional hydraulic jump type stilling basin, 
consisting of a level apron 88 ft. wide and 198 ft. long, surmounted by a single 
row of 6-ft. high stepped baffle piers and a 4-ft. high stepped end sill. Seven 
Sluices, 6-ft. wide by 7.5-ft. high are provided in the base of the spillway and 
discharge into the basin. The sluice inverts intersect the level apron at a 
slope of 4:1. In plan, the flaring sidewalls of the sluice outlets intersect the 
spillway bucket. The sluices are separated into groups of 2, 3 and 2 by two 
longitudinal dividing walls, 14-ft. high, extending through the basin to the end 
sill. 

The action of the dividing walls is particularly striking under some condi- 
tions of operation. As shown in Fig. 10, three distinct water levels are evi- 
dent within the basin. In this photograph, the three center sluices were oper- 
ating with a net head of 51 ft. The walls function to prevent the formation of 
large eddies within the basin, thereby allowing maximum spreading of the 
Sluice flows. In addition, the amount of gravel carried into the basin is held 
toa minimum. Stilling basin action is satisfactory, although with two end 
sluices operating alone, there is a noticeable concentration of high-velocity 
flow for a considerable distance downstream from the basin. 

The basin was dewatered and inspected in 1952 after almost 9,000,000 
acre-ft. of water had been released through the sluices. Rough concrete was 
found at the outlet of the third sluice, with maximum erosion depth of 5 in. 
Reinforcing steel was not exposed, however. Erosion depths of 4 in. were en- 
countered at the outlets of the two adjacent sluices. There was practically no 
erosion at the end sluices. Since there was no relationship between the quan- 
tities of flow released and the erosion, and in fact the maximum erosion had 
occurred at the sluice which had discharged the least, it was concluded that 
the characteristics of the various concrete pours in the basin were the pri- 
mary factors contributing to the amount of erosion encountered. The baffle 
piers were in good condition, with only slight erosion on the upstream faces of 
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the baffle piers located in the center of the basin. No repairs were made to 


the stilling basin, but another inspection will be made in the next few years 
as a precaution. 


Norfork Dam (Arkansas) 


This project has a concrete gravity ogee spillway with conventional hori- 
zontal-type stilling basin 568 ft. wide by 185 ft. long. Two rows of 8-ft.-high 
baffles and a 6-ft.-high end sill are provided. Eleven sluices, each 4-ft. wide 
by 6-ft. high, are evenly spaced across the base of the spillway. The average 
head for sluice releases is about 170 ft. 

An inspection of the Norfork stilling basin in 1945 has been described 
elsewhere.” It was noted that approximately 2,000,000 acre-ft. of water had 
passed through the stilling basin at that time, mostly from sluice releases, 
and that some of the baffle piers had been damaged quite severely. The dam- 
age was attributed to erosion by recirculated debris which had not been 
cleaned out of the basin when the project was completed in 1944. That this 
was the case was confirmed in 1953, when the stilling basin was again in- 
spected. An additional 2,000,000 acre-ft. of water had passed through the 
Sluices since the previous inspection. Comparison of conditions in the basin 
in 1953 with photographs taken in 1945 indicated a very small amount of addi- 
tional wear had occurred in the 8-year period. The additional wear was 
noticed in slight rounding of the upstream edges of the baffle piers. Several 
rough, shallow cup-shaped eroded spots were found on the apron upstream 
from the baffles near the right side of the basin. No repairs to the baffles or 
the apron were considered necessary. 


Fort Gibson Dam (Oklahoma) 


This project was completed in 1950, except for the powerhouse. The spill- 
way consists of a concrete gravity ogee weir controlled by 30 tainter gates, 
40-ft. wide by 35-ft. high, separated by 10-ft.-wide piers. Total width of 
spillway is 1,490 ft. Tensluices, each 5.67-ft. wide by 7-ft. high, are located 
in the base of the spillway, along the center line of each of ten successive 
spillway monoliths. The inverts of the sluices are horizontal at their outlets 
through the spillway weir, and are located about 16 ft. above the stilling basin 
apron. The sluice outlets are the “jet deflector” type, with a tetrahedral 
concrete block 3-ft. high located within the exit portal. Maximum reservoir 
pool El. is 582.1, spillway crest El. is 547, basin floor El. is 477 to 485, and 
maximum tailwater El. is 551. The stilling basin as model tested and con- 
structed for the spillway design flood of 919,000 cu. ft. per sec. is a roller- 
type basin composed of a 60-ft. long apron and a 6-ft. high stepped end sill. 
Below the sluices, however, the model indicated that a somewhat longer apron 
would be desirable in order to obtain adequate stilling action and uniform 
spreading of the sluice flow before passing over the end sill. Accordingly, 
the basin was made 20 ft. longer in the section downstream from the sluices, 
with dividing walls separating this portion of the basin from the remainder. 

The stilling basin was unwatered for inspection in 1953. The basin was 
found to be in good condition except in the area occupied by the ten sluices. 
The apron below several of the sluices was eroded to a depth of 3 to 4 in. in 
spots. Tune end sill steps were eroded a like amount in several locations. 


5. Transactions, ASCE, Vol. 112, 1947, p. 101. 
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Below the end sluice on the left, eroded depths up to 8 in. were found in the 
apron, with reinforcing steel exposed in some spots. Most of the latter ero- 
sion was in the center third of the apron and from about 15 to 35 ft. from the 
basin dividing wall. In addition, the steps on the end sill directly in line with 
this sluice were eroded to a maximum depth of 6 in. 

The basin in the section of greatest wear had been subjected to a discharge 
of more than 100,000 cu. ft. per sec. in 1949 while being used for diversion 
during the construction of the dam. At that time, a construction bridge with 
steel I-beam legs was fastened to anchor bolts embedded in the stilling basin 
apron. A barge belonging to the contractor broke loose during the peak of the 
1949 flood and swept out one span of the bridge. At the time of unwatering for 
the next step of construction, it was noted that large boulders had lodged 
around the legs of the construction bridge and by continual movement had 
scoured a large, deep area near one of the legs. After the construction period 
described, each of the sluices had released about 2,000,000 acre-ft. of water 
by the time of the inspection in 1953. 

Most of the erosion is believed to have been started by the action of eddies 
transporting gravel and boulders during the diversion stage. Subsequently, 
Sluice flows deepened the eroded areas. The fan-shaped pattern of discharge 
from several of the sluices was visible in the eroded basin floor below, but 
erosion in such instances was only to a depth of about 1/4-in. Since the usual 
head on the sluices was about 55 ft., and the outlet deflectors spread the flow 
in fan-shaped patterns into tailwater depths of about 5 to 10 ft. in the stilling 
basin, the concrete mortar was evidently not resistant to water action. 

It was noted that gruut placed three inches deep in anchor bolt recesses 
for the construction bridge had washed out. 


Pine Flat Dam (California) 


Two levels of sluices are provided through the high gravity spillway weir 
at this project. The upper level of sluices discharges onto the downstream 
weir face and flip bucket basin, whereas the lower tier of sluices discharges 
through “deflector type” sluice outlets located in the downstream face, under- 
neath the lip of the spillway flip bucket. Photographs of flow releases through 
the sluices indicate satisfactory spreading of the flows. It is understood that 
some repairs were necessary in eroded portions of several of the tetrahedral 
jet deflectors in the lower tier of sluices. It is also understood that the 
Sluices have operated with maximum heads in the order of 280 ft. 


Outlet Works Stilling Devices 


Lucky Peak Dam (Idaho) 


This project was completed in 1955. The outlet works consists of an in- 
take structure, a 23-ft. diameter steel-lined tunnel, and an unusual gate- 
controlled manifold at the downstream end of the tunnel with a design capacity 
of 30,000 cu. ft. per sec. As shown in Fig. 11, outflow is regulated by six 
5.25-ft. by 10-ft. hydraulic slide gates discharging freely into the atmosphere. 
One 30-in. diameter hollow jet valve is provided for close adjustment of the 
outflow. Tunnel flows are directed by six individual flip buckets into an 
unlined stilling pool, excavated in the rock and overburden of the river chan- 
nel. Since the top of flood control pool El. is 3060 and the floors of the gate 
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passages are at El. 2827, heads up to 233 ft. and velocities up to 122 ft. per 
sec. are attained with small gate openings. Model tests were made of the 
manifold and stilling area because of the high heads and unusual nature of the 
facilities provided. 

In May and June of 1955, the reservoir was filled to a maximum El. of 
3059 by lowering Arrowrock Reservoir upstream, resulting in heads up to 
232 feet at the slide gates. Irrigation releases up to 5,000 cu. ft. per sec. 
were made through one or two gates on a rotational basis at openings ranging 
from 5 to 8 ft. during the reservoir filling operation and until July 14, except 
for the period June 28 to July 4, when three gates were operated at 3 feet. 
Releases up to 5,000 cu. ft. per sec. were made through five gates at openings 
of 2 ft. until July 26. It was then apparent that this 13-day period of operation 
had resulted in erosion of the concrete invert and walls downstream from the 
gates, so releases for the following five days or until July 31 were made from 
gate No. 1 (upstream) at 8 ft. opening. Thereafter, all releases were made 
from gate No. 3 which had been out of service previously. 

The observed flow conditions during this period indicated satisfactory dis- 
persal and stilling action of the jets. As shown in Fig. 12, the action of the 
flip buckets provided a spectacular view. 

The outlets were inspected on July 7, 1955, at which time some minor ero- 
sion, not over 1 in. deep, was noted on the walls downstream from the stoplog 
slots. This was believed to have been caused by cavitation induced by flow 
past small vertical “ridges” of concrete that protruded only 1/8 to 3/16 in. 
from the walls at the trailing edge of the 12-in. taper at the downstream end 
of the stoplog slots. Additional, more extensive erosion was noted on the left 
wall near the floor downstream from gate No. 1. The wall was eroded to a 
depth of 6 to 8 in., and extended to a height of 2 ft. above the floor for a dis- 
tance of 3 ft. downstream from the gate frame. Some pitting was also noted 
in the gate frames and seals nearby. 

The outlets were inspected again on August 4, when it was found that opera- 
tion with 2-ft. gate openings had caused extensive erosion of the level con- 
crete floor and flip buckets downstream from the gates, as shown for gate 
No. 1 in Fig. 13. Fig. 13(a) shows that the eroded areas extended in sym- 
metrical patterns from the downstream edge of the gate frames. Maximum 
depth of erosion was 40 in. 

As a result of these inspections, it was concluded that erosion was not 
caused at gate openings of about 3.5 ft. and larger. At small gate openings, 
cavitation pressures are apparently produced by penetration of the downward 
jetting flow from each gate slot into the normal horizontal sluice flow, with 
rebound on the invert and resultant low-pressure areas, similar to the jet re- 
bound effects noted elsewhere.4 For larger gate openings, the jets apparently 
cannot penetrate the greater depths of horizontal sluice flow with sufficient 
force to produce pressures in the cavitation range. It should be noted that the 
gate lip was of standard design with 45 degree bevel to the lowest tip of the 
gate, located at the downstream edge of the gate slot. 

Apparently, the jet action from the gate slots was a result of a 4 in. offset 
at the downstream edge of each slot, since similar erosion has not occurred 
at other projects with comparably high heads where there are no such offsets. 
The eroded concrete was a 4-bag mix with 3-in. maximum size aggregate. 


4. Boulder Canyon Project Final Reports, Bureau of Reclamation, Part VI, 
Bulletin 1, “Model Studies of Spillways,” 1938, pp. 163-182. 
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A 12-in. strip of concrete located adjacent to the downstream edge of each of 
the gate frames, composed of a 6.25-bag mix with 3/4-in. maximum aggre- 
gate had practically no erosion. 

Repairs to the outlets consisted of installing 3/4-in.-thick steel plate 
liners on the floors and walls of the outlet channels. At gate No. 1, the lining 
extends to the end of the flip bucket whereas on the other gates it extends on- 
ly to the ends of the piers. Piezometer taps were placed in the floor and left 
wall at gate No. 1, so that prototype pressure data could be obtained. 

Two hydraulic models of the outlets were tested during the design stage, 
one a comprehensive model of the manifold, flip buckets and a portion of the 
river channel downstream, to a scale of 1:28.75, and the other a 1:16 scale 
model of the main conduit and one branch lateral, including gate slots and 
flared walls downstream. Although cavitation-producing pressures were 
noted with the 1:16 scale model in the vicinity of the gate slots with larger 
gate openings (but not at full opening) positive pressures were recorded for 
gate openings below 7 ft. The model did indicate the tendency of the gate at 
part opening to deflect the flow downward and outward into the gate slots. 

Prototype pressure tests were made in July 1956, utilizing the piezome- 
ters placed in the floor and left wall below gate No. 1 during the repair period. 
Two longitudinal rows of 5 piezometers each had been placed in the floor at 
distances downstream from the face of the gate seal of 5, 6.5, 9, 11, and 14.5 
ft., respectively. Row A was placed 6 in. to the left of the gate centerline. 
Row B was placed from 24 to 30 inches to the right of the centerline, in the 


area where the erosion pattern had originated. Some of the test results are 
shown below: 


AVERAGE PIEZOMETRIC PRESSURES IN FLOOR BELOW GATE NO. 1 


Ave e Piezometric Pressure in ° ater 
I-Ft. Gate Gate 3-Ft. Gate Gate 5-Ft. Gate 
ning ni Openi ni | Openi 
Row A Row B Row A Row B Row A Row B Row A Row B Row A Row B 


4.3 49. 


As expected, pressures were lowest at 1- and 2-ft. gate openings, and at 
these critical gate openings the pressures for the upstream piezometers in 
Row B were substantially lower than the corresponding piezometers in Row A. 

One piezometer was located in the left sidewall 1.5 ft. above the floor and 
3 ft.-8 in. downstream from the face of gate seal. The pressures obtained at 
this piezometer corresponded closely to those obtained with the first two 
piezometers of Row A. 

In the author’s opinion, these data serve to verify the earlier supposition 
that the concrete erosion was caused by fluctuating pressures in the cavitation 
range induced by high-velocity jets originating in the gate slots. 


Denison Dam (Texas and Oklahoma) 


The outlet works at this project consists of three 20-ft.-diameter circular 
conduits which discharge into an upper stilling basin terminated by an ogee 


oy 
‘+ 
fr Gate 
Seal, ft 
5 2 #21.7 417.7 
6.5 1 #10.1 
9 3.1 -5.3 -0.1 -0.3 -0.5 /3.2 
11 46.1 42.6 -1.5 £3.6 #1.5 = 
14.5 71.8 40.70 43.3 f1.6 42.1 -1.0 #0.4 
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13b. Erosion at Lucky Peak Dam. Eroded Area in Flip Bucket. 
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weir, and thence into a lower stilling basin. Spreading of the flow from the 
conduits into the upper basin is aided by divisional piers approximately 80 ft. 
long. The two-level stilling basin was adopted in order that the basin and 
outlet channel might be founded on a thin layer of Goodland limestone, rather 
than on the Walnut marly clay and Trinity sand below. Model tests were 
made of the adopted design in 1939 and 1940. 

As shown in Fig. 14, the upper basin is 335 ft. long and a maximum of 172 
ft. wide, and has two rows of 10-ft.-high baffle piers. The up-stream row of 
baffle piers lies 209 ft. from the outlet portal. At the end of this basin the 
weir crest rises 19.5 ft. above the floor. Design flow is 66,750 cu. ft. per 
sec. 

The lower basin is approximately 20 ft. below the upper, and extends 140 
ft. downstream from the intermediate weir. Width of the basin is constant at 
172 ft. Two rows of 8-ft. high baffle piers and a 7-ft. end sill are provided 
in this basin. 

In 1947, prototype pressure and air-demand tests of the conduits were 
made. The maximum release was 60,000 cu. ft. per sec., or 90% of design 
flow. As shown in Fig. 15, the hydraulic jump action was satisfactory but 
very turbulent. Observation of the high-velocity flow (about 65 ft. per sec.) 
emerging from the outlet portals indicated that the upstream row of baffle 
piers was absorbing large impact forces, since the hydraulic jump began 
forming quite near the baffle piers. Good mixing action was observed within 
the hydraulic jump. Occasional surges overtopped the 44ft. high walls of the 
upper basin, but small quantities of water were involved and no damage was 
done. Flows over the intermediate weir, although rough, were surprisingly 
uniform considering the turbulence of the jump in the upper basin. Stilling 
action in the lower basin may also be described as quite turbulent. Down- 
stream from the end sill, strong currents continued for some distance. Pro- 
jecting river banks about 2,000 feet downstream suffered some damage 
through undercutting and sloughing. The model had indicated bottom veloci- 
ties of 16 to 18 ft. per sec. below the end sill for this discharge. 

For the condition of unbalanced flow into the upper basin from one or two 
conduits, the baffle piers and stilling basin weir are subjected to severe 
impact. 

As closely as can be determined, the model tests had predicted the proto- 
type performance of this stilling basin in every respect. The 1940 model test 
report states that the hydraulic jump was held in the upper basin by the baffle 
piers, and that to overcome this condition would have required excessively 
high sidewalls. Also, the report states that flow in the exit channel beyond 
the end sill had relatively high bottom velocities of 13-16 ft. per sec. for both 
small and large discharges. Since the ability of the limestone to withstand 
the sustained velocities was questioned in the report, about 50 ft. of concrete 
paving had been provided below the prototype end sill to improve this condi- 
tion. 

The upper and lower basins were inspected in 1948, after 4 years of opera- 
tion and 8,500,000 acre-ft. of water releases. The upper basin was in excel- 
lent condition and no evidences of pitting or cavitation was evident. There 
were a few chipped spots at random locations on the baffle piers, which were 
believed to have been caused by large steel rollers from a broken roller chain 
on one of the Broome gates. Although the rollers were made of hard steel 
alloy, they were rounded into near-spherical shapes when recovered. 
The lower basin was not dewatered but was inspected by a professional 
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Fig. 15. Stilling Action at Denison Dam with 90% of Design Discharge, 
(a) Upper Basin, (b) Lower Basin. 
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diver and by engineers using the diving gear. The basin was found to be in 
good condition and the only notable evidence of wear was found on the floor 
at the downstream edges of the end baffle piers, where a small amount of pit- 
ting was noted. A considerable number of large rocks were found between 

the baffle piers and the end sill. The diver reported one rock to be in the 
order of 8 to 10 feet long and 3 to 4 feet thick. 

The limestone exit channel was carefully inspected but no signs of wear 
were noted from the sustained high velocities referred to above. 

In the Red River channel downstream from the limestone ledge at the out- 
let works, degradation of the alluvial bed has been progressive over the years 
since construction of the dam was initiated in 1939. The effect on tailwater 
elevations at the outlet works in the 15-year period 1939 to 1954 has been a 
lowering of 1.5 ft. at 40,000 cu. ft. per sec., 4.5 ft. at 20,000 cu. ft. per sec., 
and 6.5 ft. at 5,000 cu. ft. per sec. Little additional degradation is expected 
in the future, since minimum flow releases from the power plant have been 
observed to have only about one ft. of slope in the five miles downstream to 
the river control section. 


Kanopolis Dam (Missouri) 


This project was completed in 1948. The outlet works consists of a tunnel 
14 ft. in diameter and a horizontal jump type stilling basin with two rows of 
5-ft.-high baffle piers and a 4-ft.-high end sill. From the tunnel portal, the 
sidewalls flare in plan at 3:1 until a basin width of 52 ft. is reached at the lo- 
cation of the downstream row of baffle piers, and then extend parallel for the 
remaining distance to the end of the basin. At the location of the upstream 
row of baffle piers, the stilling basin is 40 ft. wide. The length of sloping 
apron downstream from the portal is 42 ft., whereas the horizontal portion of 
the apron extends for a distance of 78 feet. Capacity of the outlet works at 
full reservoir pool is 5,700 cu. ft. per sec., but normal maximum releases 
are limited to about 3,000 cu. ft. per sec. 

A progressive scour problem existed in the outlet channel below the basin. 
A small channel was eroded through the 18-inch rock paving, and some ero- 
sion of the banks had occurred by 1952. In order to stabilize the outlet chan- 
nel, the 18 in. depth of grouted rock paving (on 6 in. of gravel) was repaired 
for a distance of 200 ft. downstream from the basin, and a concrete cut-off 
wall, 8 ft. deep by 3 ft. wide, was placed 100 ft. downstream from the basin. 

Field observations indicate that a hydraulic jump will occur for discharges 
up to the normal maximum of 3,000 cu. ft. per sec., but for a flow of 5,500 cu. 
ft. per sec. photographs indicated considerable spray action from the impact 
of tunnel discharges on the upstream row of baffle piers. The result was in- 
complete stilling action and considerable turbulence in the outlet channel. 
Later inspection revealed that the sides of the upstream row of baffle piers, 
as well as the floor immediately adjacent, were pitted to a depth of 1 to 2 in. 
The evidence as shown in photographs clearly indicates that the pitting was 
the result of cavitation. A rough check indicates that tailwater depth was 
about 0.7 of the hydraulic jump depth, Dg, at a discharge of 5,500 cu. ft. per 
sec. 


Muskingum Reservoirs (Ohio) 


During the 1940 flood in the Muskingum River Basin, observations and 
tests were made at the outlet tunnels and conduits for Atwood, Bolivar, 
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Leesville, Mohawk and Wills Creek Dams. These projects, with the exception 
of the latter, have conventional hydraulic jump stilling basins with level apron, 
two rows of baffle piers and end sill. The Wills Creek basin has one row of 
baffle piers and no end sill. Because of limitations on the channel capacities 
downstream, it was not feasible to open all gates at these projects so as to 
observe hydraulic action under conditions approaching maximum capacity of 
the outlet works. However, flows of 16 to 71 percent of maximum capacity 
were observed. In each case, the hydraulic jump was well within the basin, 
with good energy dissipation and well-distributed, acceptably low velocities 
over the end sill. Photographs taken of the stilling action at Bolivar and 
Mohawk Dams were compared with similar photographs of hydraulic model 
tests conducted in 1935. In each case, the prototype performance was equal 
or better than that anticipated as a result of the model tests. 


DISCUSSION AND CONCLUSIONS 


The performance record of the stilling devices at Corps of Engineers 
reservoir projects has been good. Stilling action has been equal or better 
than anticipated during the design and model testing stage, except in the few 
instances described above. Some erosion of the energy dissipating structures 
has been experienced at several projects, and these instances have been dis- 
cussed in this paper because of their interest and the possibility that others 
may apply the results to the prevention of similar erosion on future projects. 
However, most of the structures have withstood the erosive forces of high- 
velocity flow releases with little or no damage. 

Of the stilling basins below spillways at Corps of Engineers projects which 
have experienced large discharges, probably the most severely tested of the 
group has been the Bonneville Dam basin, with an experienced discharge with- 
in 62 percent of design discharge and with about one billion acre-ft. of total 
flow during past years. Stilling action has been adequate to date, but recent 
modifications to the stilling basin will improve the basin performance. 

Stilling basins below outlet works have been tested near design discharge 
more frequently. Model-prototype conformity has been good. 

It is interesting to note that the causes of the erosion experienced include, 
in the author’s opinion, the following: Apron not long enough for complete 
stilling action, apron not deep enough for the formation of a fully effective 
hydraulic jump, structural shapes within the stilling basin which created 
areas of cavitation action, abnormal conditions brought about by diversion of 
river flows during construction, inadequate stilling action brought about by 
other, unforeseeable factors and misconceptions which were prevalent in hy- 
draulic design procedures at the time the project was designed, including in 
some cases model testing procedures, inability of the concrete to resist the 
erosive forces of the water action in cases where “adequate” stilling action 
was provided by the basin design, and combinations of the above. 

No damages were reported from inadequate structural design. Many of 
these types of structures are probably somewhat overdesigned as compensa- 
tion for some of the unknown forces which may be inherent or developed in the 
structures. Exceptions to this statement may be instances such as the deter- 
mination of the centrifugal forces developed in bucket curves where the forces 


i= 
| 
a 
a 
= 
a 
# 
2 


ASCE BERRYHILL 1264-35 


involved are larger than has been generally realized, as noted by D. B. 
Gumenski, M. ASCE.§ 

In the author’s opinion the experiences presented in this paper lead to the 
conclusion as suggested by Gumenski that the hydraulic design of these struc- 
tures is the most important single phase in their development. Although the 
engineer responsible for the hydraulic design of a project receives from 
model tests important and at times indispensable help, particularly when no 
adequate means of analysis exists (and in any case the model aids in verify- 
ing the adequacy of the design), the primary responsibility for acceptance of 
a particular design must rest with the hydraulic designer. He must determine 
what actual conditions may develop in the prototype which cannot be foreseen 
by model tests. 

Baffle piers at Corps of Engineers projects have withstood erosion quite 
well in most instances. Whereas the erosion of the baffle piers at Bonneville 
Dam was understandably severe, as described above, and erosion to a lesser 
extent occurred in the first few years at Norfork Dam and several other pro- 
jects where gravel deposits were left adjacent to or in the basin upon com- 
pletion of construction, baffle piers of many projects have had little or no 
erosion. An example of the latter is Denison Dam, where the baffle piers 
are subjected to heavy impact forces and yet have had very little erosion to 
date after passing 15 million acre-ft. of flow. 

For standard hydraulic jump type basins composed of a level apron, two 
rows of baffle piers and end sill, an apron length measured from the P.L. of 
the spillway bucket curve to the lower end of the end sill equal to 3Dg appears 
to provide adequate stilling action for most of the projects studied, insofar as 
the limited prototype information at or near design flows indicates. Excep- 
tions to this criterion were the McNary Dam basin with low Froude Number, 
and basins on especially erodible foundations, where a somewhat longer length 
was considered desirable. 

It has been found desirable to unwater many of the stilling basins periodi- 
cally for inspection. Substantial quantities of rock and gravel have been re- 
moved during some of these inspections. 

In some recent instances, conventional hydraulic jump type stilling basins 
with baffle piers and end sill have been designed for full energy dissipation 
for the greatest flood of record or a somewhat larger flood, while not allowing 
the hydraulic jump to sweep out of the basin for the spillway design flood. A 
high end sill aids in confining the jump to the basin for this condition. In 
terms of the spillway design flood, the resulting basin dimensions can often 
be reduced to provide a basin length of 2 1/2 Dg and apron floor set high 
enough that the tailwater depth equals about 0.8D92. 

With regard to the height of basin training walls for hydraulic jump type 
basins, experience with hydraulic models (and a few prototype experiences 
such as at McNary Dam) indicates that some overtopping of the walls, say 
about 5 feet for the spillway design flood, is not harmful to the stilling action. 
If the basin were closely confined by an earthen dam, this reduction would not 
be made. 

Present practice of the Corps of Engineers provides for superior quality 
concrete in critical areas of hydraulic structures. Where velocities exceed 
40 ft. per sec. frequently or for prolonged periods of time, the water-cement 


6. “Design of Side Walls in Chutes and Spillways,” by D. B. Gumenski, 
Transactions, A.S.C.E., Vol. 119, 1944, p. 355-372. 
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ratio is held to a maximum of 5 gallons per bag, slump at about 2 inches, and 
maximum size of aggregate at 1 1/2 inches. 
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SYNOPSIS 


The recirculation of water after it has been used as a coolant in thermo- 
electric plants and returned to the river or canal from which it was original- 
ly drawn is of great importance to the economy of operation of such plants. 
Among the many factors that affect the amount of recirculation, the most im- 
portant ones are the channel depth and discharge, the distance between intake 
and outlet of the flow diverted for cooling, the degree of heating of the divert- 
ed water, and the relative amount of diversion. The effects of these factors 
have been singled out for an extensive analytical and experimental study, the 
results of which, together with design criteria derived therefrom, are 
presented in this paper. 


I. INTRODUCTION 


A part of the water that has served as a coolant in a steam plant and been 
returned to the river or canal of origin often finds its way back to the intake 
in the form of an overlying wedge-like layer, from which it is once again 
taken (together with new water) through the cooling system. The amount of 
recirculation of heated water, which directly affects the efficiency of opera- 
tion of a plant and often lowers its maximum capacity, depends on the geome- 
try of the intake, the design of the outlet, the distance between the intake and 
the outlet, the discharge and mean depth of the original flow in the river or 
canal, the relative amount of diversion, and the degree of heating of the di- 
verted water. The complexity of the recirculation problem necessitates the 
separation of these factors into groups, the effects of which can be studied 
separately and then compounded subsequently. 


Note: Discussion open until November 1, 1957. Paper 1265 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
83, No. HY 3, June, 1957. 
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In the present investigation, the last four factors were chosen for careful 
analytical and experimental study, because they form an organic unit of pre- 
dominant importance in the determination of recirculation. To eliminate the 
effects of any particularity of the design of the intake and outlet, the simplest 
forms were used for these structures: for two-dimensional flow the intake 
was a bottom slot extending all the way across the concrete channel used to 
Simulate a canal in the laboratory; for three-dimensional flow the intake was 
a rectangular side opening with its bottom flush with that of the channel; in 
most cases the outlet extended all the way across the channel and discharged 
the heated water in such a way that downstream from the outlet the tempera- 
ture was essentially uniform. In the few available field tests, the temperature 
distribution changed continuously.(10) 

As will be seen from the body of the paper, the amount of recirculation is 
directly and simply related to the areas occupied by the warm water and the 
cool water, respectively, at the intake—i.e., to the position of the common in- 
terface of the warm and cool layers. If the intake and outlet conditions are 
essentially the same as those prevailing in the experiments, the problem of 
recirculation reduces to the problem of determination of the form of the warm 
wedge, and the results presented herein can be directly used either for pre- 
diction of recirculation or for selection of the distance between the intake and 
the outlet. However, the question then arises: how useful would the results 
be for other designs at the outlet and the inlet? 

The amount of recirculation for any particular geometrical design of the 
intake depends chiefly on the amount of diversion, the degree of heating, and 
the location of the interface at the intake—the last of which, for a first ap- 
proximation, does not depend on the particular intake form used. Thus the 
method for determination of the form of the warm wedge presented herein is 
directly useful for the investigation of the effect of the intake design. As to 
the outlet, its design can be expected to affect the form of the warm wedge 
only through its effect on the local position of the interface, the equations 
governing the form of the warm wedge being the same for all outlet designs. 

Thus the determination of the form of the warm wedge not only provides a 
means of predicting the amount of recirculation for intake and the outlet condi- 
tions similar to those prevailing in the experiments, but also plays a definite 
and specific role in the determination of the effects of intake and outlet de- 
signs on the amount of recirculation. In short, the method for the determina- 
tion of the warm wedge presented herein occupies a central position in the 
solution of the recirculation problem. 

Before proceeding, the notion of river or canal, to which the present paper 
applies, must be understood in the sense that the Froude numbers are not too 
small—in particular not less than those of the numerous experiments cited 
herein and of the other field tests to which reference is made. Otherwise the 
waterway behaves like a lake and the intake conditions change completely, as 
analyzed theoretically by Craya(1) and experimentally by Gariel. 


Il. ANALYSIS 


The flow in the canal can be divided (see Fig. 1) into three zones: that up- 
stream from the intake, that between the intake and the outlet, and that down- 
stream from the outlet. In the first zone the warm wedge is stagnant and the 
fluid in it has no mean velocity. For the particular outlet used in the 
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experiments, the third zone is unimportant, because thorough mixing at the 
outlet destroys downstream stratification. However, the case of stratifica- 
tion in zone three is included in the analysis. The determination of the form 
of the warm wedge—which will be the chief concern of the analysis, because 
this form is intimately related to the amount of recirculation—is complicated 
by the relative uncertainty of the interfacial shear and the apparent lack of 
starting points for the determination of the interface. It will be shown in the 
following sections that the critical depth (yet to be defined) of the lower layer 
occurs at the intake and the outlet, thus providing the starting points needed 
for the upstream and the middle zones, respectively. With a somewhat ration- 
al evaluation of the interfacial shear, the one-dimensional equations of motion 
for the two layers then permit the determination of the interface. Thus the 
following seemingly unrelated sections are individually necessary and all to- 
gether sufficient for this determination, and must not be mistaken for discon- 
nected considerations presented without a central theme. 


A. Equations of Motion 


The flow in both the upper (warm) and lower (cool) layers will be in gener- 
al non-uniform. If the vertical component of the velocity is neglected and on- 


ly the mean velocities in the respective layers are considered, the equations 
of motion are 


for the upper layer and 


for the lower layer. In these equations, which were given by Schijf and 
Schénfeld,(5) x is the distance measured from a fixed point to any variable 
point on the interface in the downstream direction, hy and hg are the depths 
of the warm and cool layers, uy and ug are the corresponding mean velocities 
(positive in the downstream direction), p and p +4 represent the densities 
of the upper and lower fluid, respectively, and g is the gravitational accelera- 


tion. Sg is the slope of the canal bottom; the energy slopes, S;,_ and S2e¢, are 
defined by 


Sie= Pgh, Soe “egha (3) 


in which the shear stresses, T 4 at the bottom and Tj at the interface, are 
given by 


=f; Ug | U2) Ce | (4) 


It should be understood throughout this analysis that the density differences 
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considered will be so low (Ap/p< 0.003) that the variations in the free- 
surface elevation (i.e., in the total depth h) will be negligible in comparison 
with the variations of the position of the interface. In other words, one may 
put h = hy + hg = const. 


By subtracting Eqs. (1) and (2), and utilizing the equation of continuity and 


Eqs. (3) and (4),(9) one can reduce the differential equation for the position of ~ 
the interface to the dimensionless form rT: 
- (1-7) - (-1) +Ca” 
in which 7) = ho/h, a= {i/f, qi/q2, and 


the latter being the densimetric Froude number (after Keulegan(2)), based on 
the approaching canal flow (ug denotes the velocity of the undisturbed flow). 
Integration of Eq. (5) can be performed after the proper end conditions for 
the different zones have been determined. It will be more convenient to cal- 
culate the distances x from the downstream end section upwards. Then the 
sign of x has to be changed in Eq. (5). 


B. The Critical Depths 


In connection with the determination of the form of the warm wedge, the 
critical depth is defined to be that for which the tangent to the interface is 
normal to the bottom of the channel—i.e., 


A(x/h) (7) 


From Eq. (5) it follows directly (the subscript c indicating the critical condi- 
tion) that 


2 

which was obtained already by Stommel and Farmer(3) in a different way. In 
the physically real range of 0<"<1, and for Ca # 0, represented in Fig. 2, 
there are two critical depths for a given Froude number Fo: a lower anda 
higher one, between which, in accordance with Eq. (5), the interface rises in 
the upstream direction. If the value of (FoQ)max is exceeded, no increase in 
the elevation of the interface in the upstream direction can take place. For 


Cg = 0, which corresponds to a stagnant upper layer, the critical depth is 
given by its usual definition: 
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In computing the critical depths, the total depth h was assumed to be con- 
stant. In practice h is always known, because the canal is designed for a cer- 
tain depth of homogeneous fluid. Otherwise (if the value of h were not given), 
there would be an infinite set of paris of solutions for Eq. (8) as proved by 
Benton. (4) 

As will be shown in the next article, of the two critical depths the lower 
one will be established immediately upstream from the outlet as that which 
corresponds to the end of the stratified flow. Upstream from the outlet the 
depth of the lower layer is increasing, but it never can exceed the upper 
critical value. 


C. Conditions at the Outlet 


In the zone downstream from the outlet the stratification ceases to exist, 
due to the thorough mixing produced by the form of outlet used in the experi- 
ments and a homogeneous fluid of density p * flows farther downstream. 
Tentatively one may assume that the middle zone has to end with the critical 
depth or with supercritical flow—in other words, that the conditions down- 
stream from the outlet do not have any influence on the middle zone. This 
can be justified by finding the downstream depth at the outlet, which would 
correspond to the critical depth upstream from the outlet. 

As indicated in Fig. 3, single primes and double primes denote the quanti- 
ties upstream and downstream from the outlet, respectively. This notation 
will be used henceforth unless otherwise stated. By applying the momentum 
principle to the lower layer, as done by Yih,(6) one has 


2 


Here qo‘ = qg", because the discharge in the cool layer does not change 
across the outlet. For the upper layer one has to assume that the returned 
warm water while flowing through the lateral canal has no momentum in the 
x-direction, and that both discharges q,‘ and q;"’ are directed outward from 
the “control area.” Therefore, the impulse-momentum equation for the upper 
fluid will be 


hy 


Downstream from the outlet one may assume without appreciable error a 
uniform velocity distribution, and by considering that hy' + hg = hi" + hg" = 
h = constant, one obtains 
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(2) 


Equations (10) and (11) specify hg '/h as a function of qg"/qy"" and qy'/q2'. 

In other words, for every value of the upstream parameter qyj'/q2', which 
represents the ratio of the discharges in the middle zone, there exists a cer- 
tain ratio of the downstream discharges q2''/qj" (or a definite value for 
hg'/h) for which the depths immediately upstream and downstream from the 
outlet are conjugate. If the actual ratio of qg''/qj"' would be greater (or the 
actual value of hg'/h would be less) than indicated by the conjugate state, the 
conditions downstream from the outlet would influence the middle zone by 
submerging it. 

The upstream parameter qj'/qg' can vary between the limits 0<q,'/qo9'< 1, 
where the value of 1 corresponds to 100% diversion and is more conducive to 
submergence than the value of zero. The simultaneous solution of Eqs. (10) 
and (11) then gives(9) for the conjugate state 


= 0001 (13) 


ha 204; 
= 0.4; 


These conditions practically never appear in practice, and as a consequence 
at the downstream end of the middle zone the critical depth will be 
established. 

The case in which the two layers of densities p ; and 99 are still present 
downstream from the outlet—as is often true in field conditions—can be 
treated in a completely similar manner, by solving simultaneously the mo- 
mentum equations of the upper and lower layers. Due to the fact that the 
three densities p;, Py, and o9 are very close to unity, the result obtained 
will be practically identical with Eq. (13). In other words, at the downstream 
end of the middle zone the critical depth is also established for the case that 
the stratification downstream from the outlet is not destroyed, thus again giv- 
ing the necessary starting point for the calculation, at the outlet. 


D. Conditions at the Intake 


In the neighborhood of the side intake used, the flow in both layers can be 
considered as a potential flow in a manifold efflux,(12) at least at elevations 
far from the interface. In the lower layer the flow is approaching from the 
upstream direction, whereas in the upper layer the direction of the approach 
is opposite to it (Fig. 4). At the mouth of the intake there are corresponding - 
ly two stagnation points, and the separation from the walls of the lateral canal 
is on opposite sides for the two layers. Near the stagnation points the veloci- 
ties are rather low and the pressure in the stagnant layer is increased by the 
velocity head, whereas in the moving layer the pressure does not change very 
much from the initial value. The created pressure difference will result ina 
tendency for the lower fluid to rise near the downstream corner of the intake, 
and an opposite tendency for the upper fluid to sink near the upstream corner. 
This will cause a secondary helicoidal motion near the intake entrance, with 
the consequence of mixing in the lateral. The “penetration zone” will be 
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greater at the downstream corner, because the velocity of approach of the 
cold layer is always greater. 

If the transition at the intake is treated as a sudden enlargement (as indi- 
cated by the shaded area in Fig. 4), an analysis(9) similar to that of Jaeger 
for a homogeneous fluid(11) gives 


-1) (14) 
2(x-A) 


In this equation 1) = hg '/ho denotes the ratio of the two conjugate depths, h9' 
being the depth downstream from the intake, the quantity 


2 (15) 


is the local densimetric Froude number, and the parameter # denotes the 
ratio of the discharges in the lower layer: 8 = q2'/q2. 

Another approach for the transition at the intake would be that applied for 
the neighborhood of the outlet, where the whole width of the canal was includ- 
ed in the calculation and the four discharges (upstream and downstream, for 
each layer) were different. The application of the momentum principle(9) 
would then result again in Eq. (14) with the only difference that now the 
parameter @ is replaced by B2. 

It seems that the sudden-enlargement method is more convenient for ap- 
plication to the side-intake conditions, whereas the discharge-variation 
method (second approach) applies more conveniently to a bottom-slot intake 
which extends across the whole width of the main canal. At the outlet the 
warm water is returned to the canal in a jet-like form, which assures that 
the discharges are distributed more or less along the whole width of the canal, 
so that the second approach is more appropriate. 

The analysis and graphical representation of Eq. (14) were given by 
Jaeger(11) and the graph is reproduced in Fig. 5. If the upstream conditions 
are given, together with the parameter 8, which actually represents the di- 
version ratio, the other conjugate depth can be determined from this figure. 
Assuming critical conditions upstream from the intake (F9 = 1) and taking an 
average 8 = 0.5, the downstream conjugate depth would then be around 50% 
higher than the critical depth upstream from the intake. This oftentimes will 
be impossible, because the depths in the middle zone are limited—i.e., they 
have to be less than the “upper critical depth,” given by Eq. (8) and in Fig. 2. 
The transition then occurs in such a way that a supercritical flow of the lower 
layer is established downstream from the intake, followed by a jump (Fig. 1). 
Figure 2 indicates that the jump will be in close proximity to the intake. 

In brief, the depth of the lower layer immediately upstream from the intake 
has to be critical, and this gives the starting point for the calculation of the 
upstream zone. 


E. Form of the Warm Wedge 


From the earlier analysis it follows that the flow conditions in all three 
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zones are independent, and hence the form of the warm wedge can be deter- 
mined in each of them separately. Since it is known that the critical depth is 
established at the ends of the upstream and middle zones, from these points 
on the calculation can be performed by using Eq. (5). 

For the upstream zone the upper layer is stagnant (Cg = 0). Substituting 
the value of Eq. (9) into Eq. (5) for x = 0 (x is measured upstream from the 
intake) and integrating, one has 


4X (14%) 


It was assumed that the resistance coefficients and the Froude number do not 
change along the warm wedge. By putting 7 = h2/h = 1, the total length L of 
the warm wedge in the upstream Zone is obtained (Fig. 6). (The experimental 
points in this and the next figures are the results of the experiments dis- 
cussed in greater detail in Section II.) 

The other extreme case is Cg = 1, which corresponds to 100% diversion, 
the discharges in the two layers then being equal. Integration of Eq. (5) for 
this case gives 


2 
~ + (1-9) -1 2(4-n)- 


where the letters A and B denote symbolically the lower and upper limits. 
Equation (17) is represented graphically in Fig. 7 for a common experimental 
value of 0.150 for Fo'—the Froude number of the middle zone. For 0<CQ<1 
the general solution is more complicated, and the best way of carrying out 
the calculation is by numerical integration of Eq. (5). 

The length of the warm wedge in the middle zone is given in Fig. 7. There 
is a certain adjustment which has to be made in connection with the length of 
the warm wedge in the middle zone. In the downstream neighborhood of the 
intake where the expansion takes place (Fig. 4), the discharge qg' per unit 
width in the lower layer is larger than farther downstream. Therefore, the 
ratio Cg = qi'/qa2' is smaller. With decreasing Cg the slope of the interface 
according to Eq. (5) tends to flatten, so that the upward curvature of the inter- 
face (Fig. 7) disappears. The expansion zone is approximately 1 to 2 times B, 
where B is the width of the canal. The total length of the warm wedge is 
therefore increased approximately by this amount. If the initial ratio of Cg 


is close to zero, the expansion is negligible and the flattening effect can be 
neglected. 


wi 
ia 
4 
2 (n*- , { ( 3_p%), 
h F2 7 3 n ° 
co] 
a 


1265-11 


Bottom-slot intake 

o Side intake 

@ Wedge in middie zone 
+ Side intake and outlet 


00143 © 0.156 
0145 0.156 
® 0.160 


& 
| ASCE BATA — 
0.8 
aa i 

j 
0.0 i 

a | 0 q 
a 
Fig.6 Length of the Worm Wedge 
0.8 3 

8 

h Eq.(I7) for F, = 
0.4|__ 
O 
| Fo = 
© 0. 
02 4 

h 
| 0.0 10 20 30 40 
0.5 x 1.0 LS 
f+ q 
mal Fig.7 Shape of the Warm Wedge in the Middle Zone for F,=OIS a 


ay 


1265-12 


HY 3 June, 1957 


One of the principal quantities required in the numerical evaluations is the 
resistance coefficient for turbulent flow at the interface, fj, or its ratio to the 
bottom resistance @ = f/f, which will be treated in more detail in the next 
article. In accordance with the experimental evidence(8) fj varies along the 
warm wedge from fj = 0 near the critical end section to fj/f = 1 at its up- 
stream leading edge. This is especially true for the upstream zone with stag- 
nant upper layer, whereas in the middle zone the opposite flow in the two 
layers soon imposes the final value of the resistance at the interface. 

To account for the aforementioned change of f;/f along the interface one 
may assume that the actual shape of the warm wedge can be obtained by tak- 
ing different values of @ = fj/f along the wedge and combining the different 
segments into one curve. For Fo = 0.447 (Fo? = 0.20), which is the mean 
value for the experiments, the relative shape of the warm wedge for the up- 
stream zone is obtained this way(8) and it is represented in Fig. 8. Although 
the procedure is somewhat arbitrary, the curve obtained for the shape of the 
wedge is closer to actuality than the curves for constant interfacial shear. 

As the calculations show, the effect of the Froude number on the relative 
shape of the wedge is small and can be neglected in a first approximation. 

It may be mentioned that for f = 0 and Cg = 0 Eq. (5) gives the form of the 
intrusive saline wedge in an estuary, for which, by a procedure similar to 
that of the warm wedge,(8) the relative shape is given by 


5 


3 6.4 


(18) 


x_ 


Equation (18) is represented graphically (Fig. 9) for two different values of 
Fo, together with Keulegan’s experimental data.(2) One sees that the agree- 
ment is very good. 


Resistance in Established Stratified Turbulent Flow 


For the cases of flow between fixed plates (Fig. 10, No. 1), stratified flow 
(No. 2) and free-surface flow (No. 3) the hydraulic radius, R, has the values 
of hg/2, H/2, and hg, respectively, where H is the height of the extrapolated 
velocity profile (with hg<H<2 hg). The resistance coefficient, f, does not 
vary much for the three cases. 

It is convenient to assume (arbitrarily) for all stratified-flow cases (No. 2) 
the same hydraulic radius: R=hg. However, to take care of the variations 
of the hydraulic radius, then the resistance coefficient (denoted by f*) must 
be assumed to be variable. The energy gradient for turbulent flow can thus 
be expressed by 


4H/2 29 «4h, 2g 


One can further conceive that the variable resistance coefficient, f*, is com- 
posed of a constant resistance, f, on the bottom, and a variable resistance, fj, 
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at the interface—i.e., f* = f + fj—and then from Eq. (19) it follows that 


= - (20) 
Oo = H 

In order to obtain a solution for H/hg one can express the shear stresses 
for the upper and lower fluids and equate them. In the upper layer the inter- 
facial shear will remain the same if one substracts or adds a constant veloci- 
ty to the velocity profile. The interfacial shear is, therefore, the same in 
cases (a) and (b) in Fig. 11 and is given by 


f, 2 (21) 


where f, is the resistance coefficient at the interface for the upper fluid and 
uj is the velocity at the interface. The former will be discussed in more de- 
tail at the end of this section and in the section on design criteria. 

For the lower fluid one has by geometric considerations (Fig. 10) 


Equation (21) divided by 


and equated to Eq. (22) gives 


ott 


2 


(24) 


2 


where Umax is the maximum velocity in the lower layer. 

From the logarithmic velocity distribution of the extrapolated velocity pro- 
file the values of uj/Upax and u2/umax can be obtained, (8) which after sub- 
stitution in Eq. (24) give 


2- 


(25) 


a. £8 1- 254 tn 
ha” 1-2.5A-2.54 (7 - 4) - 4) 


where 


=. 
VE 


j 
j 
| 
(22) 
To ha il 
£72. (23) 
8 
4 


BATA 


Extrapolated velocity profile 


Free surfoce 


Interface 


Velocity Distribution in the Upper Fluid 


Fig.12 Resistance Coefficient at the Interface 
in Turbulent Flow 


ASCE 1265-15 Al 
| 
To To 2To 
Fig.10 Shear and Velocity Distribution in the Lower Fluid 3 
0. 
Mie 4 Umean * Uj Uj q 
Fig | 
1.0 
3 
le fi LA | 
f 
0.0 
0.0 OS 1.0 IS 


1265-16 HY 3 June, 1957 


Equation (25) together with Eq. (20) gives the solution for the interfacial re- 
sistance which is plotted in a graph in Fig. 12. A slight variation can be ob- 
served with the absolute value of f. Equation (25) is valid for stagnant upper 
layer. In case the two layers are moving in the opposite directions, the 
velocity distribution is approximately the same as if the interface were re- 
placed by a fixed wall, which is represented in Fig. 12 by a direct-correlation 
straight line between the interfacial resistance coefficient fj and that of the 
upper fluid, fj. 

The problem now lies in the determination of the coefficient fj—i.e., what 
resistance is acting on the upper fluid, which is moving on the lower one. 
Without precise measurements of the velocity distribution in the upper fluid 
this question can not be answered. Because further evidence is not available, 
it is here assumed that for the experimental conditions in the laboratory 
flume the interface was smooth. Falling droplets indicated that there is much 
greater disturbance near the zero-velocity line than at the interface itself. 
One could therefore assume the formula for smooth surfaces for calculating 
f{, at the interface, whereas for the bottom resistance, f, both smooth and 
rough surfaces can come into consideration. Knowing the ratio f;/f one can 
determine with the aid of Fig. 12 the interfacial resistance coefficient sought— 
L2., {,/f. It is characteristic that even for f}/f = 1 one obtains fj/f<1, which 
indicates that there is a finite velocity at the interface. Theoretically the in- 
terface will be brought to rest only if {j/f>«. 

It is here briefly mentioned that by an analysis similar to that of turbulent 
flow, the variable resistance coefficient in laminar flow can be obtained,(8) 
and it is given by 


3+N 4 (26) 
F = 384 3+4N R 


in which R represents the Reynolds number of the lower fluid (calculated 
with 4 ho for the characteristic hydraulic radius), and the quantity 


hy Ue (27) 
N= 
denotes a parameter (“ with the corresponding indices represents the dy- 


namic viscosity). If the discharges are low (especially in the laboratory), 
laminar flow may occur in the middle zone. 


F. Recirculation 


As defined in previous works, (7) the recirculation is that percentage of the 
returned warm water which enters the intake and by which the average tem- 
perature of the intake water is raised. The same amount of cool water 
necessarily by-passes the intake and continues to flow farther downstream. 
The recirculation can be expressed as 


Ti -Ta (28) 
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where Tj is the temperature of the intake, Tj is that of the upper fluid, and 
T2 is that of the lower fluid. 


Limiting Condition for No Recirculation 


The intake entrance will not act as a constriction, because the flow through 
the intake is regulated by the gates at the end of the lateral canal and not by 
the intake entrance. The drop in elevation of both the free surface and the 
interface at the entrance of the lateral will, therefore, be negligible, and one 
can suppose that the amount of recirculation will be proportional to the area 
of the intake occupied by the warm layer at that place. This assumption was 
already made in(7), where a different approach was followed. 

For the expulsion of the warm wedge from the upstream zone, according 
to Eq. (16), it is required that Fo 21. However, for no recirculation it is 
necessary that the depth of the lower fluid downstream from the intake be al- 
so equal to the total depth. This is the case when the distance between the in- 
take and outlet is just equal to the length of the warm wedge in the middle 
zone. The limiting condition for no recirculation, therefore, is given by Eq. 
(16) and represented in Fig. 6 for Co = 0 with {L/h corresponding to the dis- 
tance between the intake and outlet. From the related Froude number, Fo’, 
the necessary discharge in the middle zone can be calculated; in other words, 
the maximum diversion for no recirculation can be determined. 

For the expulsion of the warm wedge even from the middle zone it is 
necessary that Fo’ ei. 

It is interesting to notice that for 100% diversion the condition of no re- 
circulation would be practically impossible, because in that situation there 
would be no flow in the middle zone, the interface would be horizontal, and 
the warm water would reach the intake no matter how far it might be from 
the outlet. Of course, for great distances the heat loss through evaporation 
and radiation would change the analysis. 


The Amount of Recirculation 


If the upstream and downstream depths at the intake are not equal, the 
most reasonable assumption is that the amount of recirculation will be pro- 
portional to their mean value, because gravitational effects are small for the 
range of temperature differences usually encountered except for very small 
diversions. The upstream depth is the critical depth for a stagnant upper 
layer, defined by Eq. (9). The downstream depth, t, is governed by the jump 
(Fig. 13), and it is less than the elevation of the calculated interface would 
indicate at the same place. The amount of recirculation can therefore be cal- 
culated by 


Lat 


It is difficult to estimate how much t is less than hg’ in Fig. 13. A straight- 
line variation of the t-line between the upper and lower critical depths, sug- 
gested by the fact that the middle portion of the interface and also its flat- 
tened end portion are fairly straight, gives for Ce =1 
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Outlet 


Fig.1I3 Assumption of a Straight Line for the Variation of t-Line 


t-he OL (30) 


h-2he 4-2n br 


in which AL denotes the distance between the intake and outlet and Ly =A L is 
the length of the t-line, where A is an empirical coefficient, having a value of 
1.28 for Cg = 1 on the basis of experimental data.(9) 

Introducing into Eqs. (29) and (30) the value of the critical depth, given by 
Eq. (9)—which holds approximately for the middle zone, too, for low values of 
Fo'—the amount of recirculation is obtained in the form 


From Eq. (17) the parameter f L/h is a function of Fo’ and Ca, so that the 
recirculation is 


where the functional form is given by the previous equations. Equation (32) 
states that the recirculation is dependent upon the Froude numbers of the up- 
stream and middle zones (corresponding to the two depths at the intake), the 
relative distance between the intake and outlet, f AL/h, and the ratio Cg of 
the discharges in the two layers of the middle zone. 

The validity of Eq. (31) or Eq. (32) is subject to two limitations. First, if 
the upstream Froude number Fo is greater than 1, the first term on the right 
side of Eq. (31) would be negative, and from Eq. (9) the upstream depth at the 
intake would turn out to be higher than the total depth, which has no physical 
significance. The validity of Eq. (31) is therefore limited to FoS$ 1. For Fo 
increasing above 1 but Fo' = const, the downstream depth would remain con- 
stant, and the upstream zone would tend to push the wedge away from the up- 
stream edge of the intake, by keeping its own depth also constant and equal to 
h. As a result the recirculation decreases only slowly with increasing Fo. 
Second, there is a lower limit of Fo, below which Eq. (31) does not hold. 
Equation (31) was derived from Eq. (29), so the second and third terms on its 
right side represent together the influence of the downstream depth at the 
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intake. When the sum of these two terms becomes equal to the first one, the 
recirculation becomes 


(33) 


and the upstream and downstream depths at the intake are equal. The eleva- 
tion of the interface at the intake is then mostly governed by the upstream 
depth with some influence of the downstream jump, and the recirculation 
would follow Eq. (33), instead of Eq. (31), with values somewhat below it. 

Equations (31) and (33) are represented in Fig. 16, together with the limit- 
ing Froude number, (Fo')jjm, for which the length of the warm wedge is just 
equal to the distance between the intake and outlet. Because there is no re- 
circulation for this case, (F¢Q') jim = constant coincides with the abscissa— 
i.e., with r = 0. 

It is interesting to note that the diversion ratio, defined as D = Qj/Qc, 
where the indices i and c refer to the discharges at the intake and in the canal, 
does not explicitly influence the recirculation. Generally speaking, for a cer- 
tain value of diversion, and for the other conditions remaining constant, one 
can have different values of recirculation. For example, for the limiting con- 
dition of no recirculation, by varying the discharges in the canal and intake in 
the manner that their difference remains always the same, all four parame- 
ters in Eq. (32) remain constant, but the diversion ratio changes. Only if the 
discharge in the canal or of the plant is assumed to be constant will the diver- 
sion ratio influence the Froude number Fo' in the middle zone directly. The 
role of the diversion ratio as a parameter is then justified. 


lil. EXPERIMENTAL VERIFICATION 


Experiments, for purpose of verifying and supplementing the foregoing 
analysis, were conducted in a concrete flume 132 ft. long, 5 ft. wide, and 15 in. 
deep. The lighter fluid was obtained by heating the water with steam. In 
about the half of the 72 runs(8) the intake was a regulated bottom slot across 
the whole width of the flume, whereas in the other half(9) it was represented 
by a rectangular 15-inch-wide side canal; it remained at a fixed place for all 
experiments. In one series the outlet was represented also by a 15-inch rec- 
tangular side canal, with its axis at 90° to the axis of the main canal. Paral- 
lel experiments were carried out with a discharge distributor as an outlet by 
which the warm water was distributed uniformly across the whole width. Be- 
cause of its mobility and because no principal difference was found between 
the two arrangements for the shortest distance (between the intake and outlet) 
investigated, the discharge distributor was used exclusively thereafter. From 
the distributor the water fell into a system of stilling screens at the free 
surface. 

The velocities in the canal were too low to be measured by the instruments 
available. Only the points of zero velocity were obtained by injecting dye, the 
change of the velocity direction being rather sharply marked. The tempera- 
ture measurements, on the other hand, were very thorough and precise. 
Thermocouples on a movable carriage were arranged on two verticals at 
every 0.5-in. height. Moreover, with movable thermocouples every point 
could be reached. A 48-channel indicating potentiometer enabled both great 
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accuracy (0.1°F) and rapidity in reading to be attained. A representative ex- 


periment is given in Fig. 14 and the results will be discussed in the following 
articles. 


A. Form of the Warm Wedge 


The temperature measurements indicate that there is a diffusion zone near 
the interface or dividing streamline. Because the two fluids have quite simi- 
lar characteristics, the heat is probably diffused at the same rate in both 
fluids so that the middle of the diffusion zone can be considered as a good 
first approximation for the interface. In the upstream zone the measured 
points of zero velocity coincide with the upper temperature limit, as suggest- 
ed by Fig. 11. In the middle zone, on the other hand, the zero-velocity line 
has been found to coincide with the interface, according to Section II E. 

The measurements indicate that there is a slight decrease in the tempera- 
ture of the upper layer in the upstream direction, with a sharp change at the 
leading edge (Fig. 1). The densimetric Froude number accordingly varies 
slightly along the wedge, which had little effect on the form of the wedge, as 
found in the analysis. The recirculation, however, has to be calculated with 
T , at the intake in Eq. (28), in order to obtain the accurate discharges in the 
middle zone. This temperature decrease, which is related to the mechanism 
of the diffusion process between the two layers, and which can be observed as 
a salinity decrease in the saline-wedge intrusion, has to be further investi- 
gated. 

In order to compare the analysis with the experiments one has to determine 
the interfacial shear. The ratio fj/f is a function of the Reynolds number 
based on ug and 4h9, which in the upstream zone varied from 15,000 to 70,000 
and in the middle zone from 4,300 to 35,000. For the longer wedges the 
Reynolds numbers are closer to the lower values, and for that range it was 
found by separate measurements in the flume that the resistance coefficient 
on the bottom is f => 0.03. According to Fig. 12 the ratio fj/f depends on f;/f, 
where for fj it was proposed to take the value for smooth surfaces which in 
this case is fj + 0.025. With these values fj;/f = 0.66 for Cg = 0 and fj/f = 0.83 
for Cg = 1. For the available field measurements on the Chicago Sanitary 
and Shipping Canal(10,8) with Reynolds numbers from 3 to 5 million one has 
{1/f = 0.010/0.014 = 0.70 and from Fig. 12 f;/f = 0.60. 

With the foregoing data the total length for Cg = 0 can be calculated, be- 
cause it will be well represented by an average or mean value of the inter- 
facial shear. In the calculations of the relative shape of the warm wedge, on 
the other hand, it is necessary to take care of the variations of the interfacial 
shear along the wedge, as explained in Section I E. Figures 6 and 8 show that 
the agreement between the analysis and the experiments for CQ = 0 is very 
good. 

The form of the warm wedge in the middle zone (for Cog = 1) can be verified 
only for portion of the length of the interface, because the total length is al- 
most never achieved in the experiments. In order to group the experimental 
points by the application of Eq. (17), data with approximately the same Froude 
number, Fo’, have to be chosen. The data in Fig. 7, for different discharge 
and temperature values, give good agreement with the theoretical curve. 

The data for the transition sections at the intake and outlet verify that the 
critical depth is established at both places. The critical depths at the intake— 
taken at a distance 3/2 h from the intake—are represented on Fig. 15, and they 
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show satisfactory agreement with Eq. (9) for both the bottom slot and the side 
intake. The critical depths at the outlet section can be obtained only by ex- 
trapolation, because the measurements at the outlet itself are impossible due 
to the great fluctuations in the mixing zone. These depths, although following 
Eq. (9), are not represented on Fig. 15. 

Finally the qualitative picture of the potential flow was verified in the ex- 
periments by dyeing the streamlines. The effect of the spiral motion in form- 
ing two pools of cool and warm water at the two corners of the intake can be 
very definitely seen in the field measurements. (8) 


B. Recirculation 


The experiments were performed for three different distances between the 
intake and outlet, from which the data for fA L/h = 0.93 are presented in Fig. 
16. Diagrams of recirculation for values 0.372 and 1.36 of fA L/h are avail- 
able in (9), The experimental data follow the theoretical curves fairly well. 
The limiting value of Fo' for no recirculation agrees especially well with the 
theoretical value. If the length of the warm wedge is less than the distance 
between the intake and outlet, a slowly moving eddy is formed in the upper 
layer, from the leading edge of the warm wedge to the intake. This eddy, 
however, does not affect the recirculation, which remains zero. 

In order to check if the experimental data for the three distances between 
the intake and outlet are consistent with Eq. (31), from the four parameters 
in Eq. (32) Fo, Fo; and Co will be chosen constant and only f A L/h will be 
left variable. With Fo = 1%tor which the first term in Eq. (31) vanishes), 

Cg = 1 and Fo’ = 0.150 (which is the best defined by the experimental points), 
Eq. (31) is plotted in Fig. 17, together with the experimental data. The agree- 
ment is very good, so this graph enables one to plot the recirculation curve 
for Fo’ = 0.150 according to Eq. (31), and for any distance between the intake 
and outlet on the basis of both theoretical and experimental support. The 
lines Fo' = const in Fig. 16, between Eq. (31) for Cg = 1 and Fo' = 0.150 on 
the one hand, and the abscissa r = 0 on the other hand can be obtained by in- 


terpolation, a greater accuracy being unwarranted due to the changes in CQ 
from 1 to 0. 
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The amount of recirculation in case of the two-dimensional bottom slot is 
much less, compared to the three-dimensional side intake(8) and even for 
100% diversion (which in this case seems to be the dominant parameter) it is 
only about 10%, regardless of the distance between the intake and outlet. This 
is because the slot is covered by a layer of cool water, and the warm water 
can be drawn down only through the action of occasional vortices. 

The effect of the intake and outlet forms on the recirculation will now be 
discussed. Although this problem was not thoroughly investigated, some use- 
ful suggestions, nevertheless, can be made. In one experiment the upper and 
lower halves of the intake entrance were blocked, successively, and in an- 
other one the width of the intake canal was about half of that normally used. 
The recirculation did not change more than 2 or 3% in all these cases. Thus, 
for the range of the temperature difference and the diversion ratio prevailing 
in these tests, the intake form for rivers seems to have little effect on the 
recirculation. In this connection, the reader may be referred to work by 
Craya.(1) 

At the outlet there is an intensive mixing between the warm water and the 
by-passing cool water, which results in a temperature drop between the re- 
turned hot water and the upper layer of the middle zone. This temperature 
drop—or, better, the percentage of the plant discharge which mixes with the 
cool water passing by the outlet, as proposed by Laursen and Hubbard(7)—is 
a measure of the intensity of the mixing, and it is different for various forms 
of the outlet. From Fig. 1, in which a typical temperature variation of the 
free surface—i.e., of the upper layer—is given, one sees that the mixing zone 
extends only to a very short distance from the outlet. If the calculation of the 
Froude number, Fo’, is based on the temperature outside of the mixing zone, 
the effect of the outlet form is excluded. This was done in all experiments, 
because the study of the influence of the outlet form requires a separate 
analysis. 

From a series of experiments,(9) in which the cross section and the form 
of the outlet are varied, one can conclude that if the warm water is returned 
to the river in a more jet-like form (i.e., with smaller cross-sectional area 
of the outlet) and directed toward the bottom, the mixing process is more in- 
tensive and the recirculation is less. The amount of recirculation—the other 
conditions being unchanged—varied in a range of more than 100% for the dif- 
ferent forms. Screens are very undesirable near the outlet section, and they 
should be moved upstream in the outlet canal. As a general rule the condi- 
tions for decrease in recirculation require a violent introduction of the warm 
water in the main canal, which is opposite to conditions required for naviga- 
tion. The best outlet form is best determined by specific model studies. 


IV. DESIGN CRITERIA 


The river section with greatest velocity will generally be chosen if several 
places for the plant site are available, even with the regulation of the section 
considered coming into consideration. With fixed plant discharge and with the 
best outlet form known, the problem is then usually to select the most eco- 
nomical distance between the intake and outlet. 

The safest design would require such a distance that no recirculation 
would take place even for the lowest discharges in the river. If the minimum 
discharge in the river is greater than the plant discharge, the Froude number 
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F¢' for the middle zone can be calculated and from Fig. 6 (for Cg = 0) the 
ratio f A L/h, and so the distance AL between the intake and outlet, can be 
directly determined. If the Froude number Fo' turns out to be greater than 
1, the warm wedge is expelled even from the middle zone and all warm water 
flows downstream from the outlet. On the other hand, if the minimum river 
discharge is equal to or even less than the plant discharge, recirculation can- 
not be avoided. The question is then how far upstream from the outlet should 
the intake be located from the standpoint of over-all economy. 

Several values for AL should be chosen and for each a diagram of recircu- 
lation as a function of the river discharge determined with the aid of Fig. 17, 
as explained in Sections II F and III B (see Fig. 16). From the hydrological 
data one can determine the number of days on which recirculation would 
occur. The amount of recirculation per degree and per day is related to the 
accompanying loss in power production, and an economic study will then show 
what distance between the intake and outlet represents the optimum. 

If recirculation is allowed to occur, the calculation has to proceed by trial 
and error, because in the field the temperature difference between the intake 
and outlet is constant rather than that between the canal and outlet. Also it 
was assumed in the analysis that the shear at the interface will be the same 
as for a smooth surface, which gave satisfactory results in the laboratory ex- 
periments. Under field conditions this assumption may not be valid and the 
resistance at the interface may be greater, with the final effect that the slope 
of the interface will be steeper. This, would, however, decrease the recircu- 
lation, so that the actual conditions in the field can only be safer than indicat- 
ed by the calculation and the laboratory experiments. 


V. CONCLUSIONS 


On the basis of the good agreement between the analytical and experimental 
results presented in the foregoing pages, the following conclusions concerning 
the recirculation of cooling water at a steam plant may be drawn: 

The three zones of flow under consideration—upstream from the intake, 
between the intake and outlet, and downstream from the outlet—are not in- 
fluenced by one another, and the critical depths are established at the end 
sections of the upstream and middle zones. The length and the shape of the 
warm wedge in each of the first two zones can be determined by applying the 
equations of non-uniform flow to the two layers of warm and cool water, and 
by considering that in the middle zone the elevation of the interface cannot 
exceed the upper critical depth. 

In the three-dimensional case with side intake and outlet, the amount of 
recirculation is directly and simply governed in magnitude by the depths im- 
mediately upstream and downstream from the side intake; these in turn are 
functions of the Froude numbers of the upstream and middle zones, the rela- 
tive distance between the intake and outlet, and the discharge ratio of the two 
layers in the middle zone. The functional dependence of the amount of recir- 
culation on the four parameters has been determined analytically. The condi- 
tion of no recirculation requires that the distance between the intake and outlet 
be equal to the length of the warm wedge. Oftentimes, however, this distance 
would be very great, and it is more economical to allow a certain amount of 
recirculation during the low-water stage in the river. 

The amount of recirculation for the two-dimensional case of a bottom-slot 
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intake, determined only experimentally, is much less than that for a side in- 
take, because the slot is covered by a layer of cool water which is only occa- 
sionally penetrated by warm water through vortex action. 
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NOTATION 
Except as otherwise noted, the subscript 0 refers to conditions of undis- 
turbed flow; subscript 1, to overlying warm layer; subscript 2, to underlying 


cool layer; subscript c, to critical state; subscript i, to interface. A single 


prime (') indicates zone between intake and outlet; a double prime ("'), zone 
beyond outlet. 


CQ =41/a2 

resistance coefficient 
Froude number 
acceleration of gravity 
depth of flow 

length of penetrating wedge 
rate of flow per unit width 


recirculation 


Reynolds number 


energy slope 


bed slope 


temperature 

mean velocity 

distance in downstream direction 

q2'/a2 

ho/h 

hg"/hg 


mass density 
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intensity of shear 


intensity of bed shear 
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A STUDY OF BUCKET-TYPE ENERGY DISSIPATOR CHARACTERISTICS 


M. B. McPherson,! A.M. ASCE, and M. H. Karr,2 J.M. ASCE 
(Proc. Paper 1266) 


SUMMARY 


Throughout this study the discharge was distributed uniformly over the 
channel width, obviating inclusion of the effects of spillway piers, unequal 
crest gate settings, slotted buckets, or laterally flared spillways or chutes. 
A bucket exit angle of 45° was used throughout, having been found by others 
to be most effective, under a variety of service conditions. A 45° exit angle 
has been recommended for high head dams.(1) 

Both a (1 on 1) and a (1 on 2) entrance slope were intensively investigated. 
In both instances the discharge, entrance head, tailwater depth and bucket 
radii were varied over a wide range of magnitude and of performance. 

General performance characteristics have been delineated. Required 
bucket radii for the various service conditions are specified. Energy dissi- 
pation effectiveness is demonstrated in a comparison with tailwater require- 
ments for a hydraulic jump in a horizontal channel. Information directly use- 
ful in approximating spray wall heights is given. A sample preliminary 
design calculation directly utilizes the material presented. 


Definition of Terms 


Referring to Figure 1, hj is the total available head relative to the bucket 
invert, measured as close as practicable to the bucket roller; hp is the depth 
in the bucket above the bucket invert, hg is the average surge peak depth 
relative to the elevation of the bucket invert, and hg is the tailwater depth 
relative to the bucket invert. The bucket radius is designated as R, and q is 
the discharge in cfs./foot width. 


Note: Discussion open until November 1, 1957. Paper 1266 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
83, No. HY 3, June, 1957. 
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Basis for Using Bucket Invert as Datum for hg and hg 
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Regions of high velocity downstream from the bucket tend to be concen- 
trated near the surface in this type of stilling basin. The energy dissipation 
downstream from the bucket, resulting primarily from the expansion of the 
live stream in the surge, can be regarded as a surface phenomenon for the 
most part. 

Included in Figure 2 are some of the findings from a recent model study(2) 
drawn to scale. With three original floor configurations, one solid and two 
movable, and with two stabilized scour patterns, the surface profile was un- 
changed for the same pool level at (a) and tailwater control level at (e). 

Other tests on the same model with a long horizontal movable floor at the 
elevation of the lip, and recent tests with a long solid horizontal floor at the 
same level as the bucket invert, produced the same surface profiles using the 
same hy and hg, shown at (d), as before. Although the tests were limited to a 
specific q, hj and hg, they offer some proof that the phenomenon is indeed 
concentrated at the surface. Nevertheless, the expanding surge scoured the 
movable floor. 

In the aforementioned model tests, it was found that with a movable floor 
at the level of the bucket lip, some bed material was thrown into the bucket 
roller (making several circulatory trips before being expelled). However, 
when the movable floor was arranged at or near the level of the bucket invert 
(for a horizontal distance at least as far as the surge peak), there was neg- 
ligible bed movement near the bucket. Until conclusive evidence to the con- 
trary is obtained, it is recommended that the channel floor immediately down- 
stream from a bucket be set at the elevation of the bucket invert, in general. 

Assuming that the above results were reasonably typical, all tests in the 
current study were conducted with a horizontal floor at the same elevation as 
the bucket invert. 

In a model study recently reported,(3) in which a movable bed was em- 
ployed, it was found that the point of maximum scour moved slightly down- 
stream with increasing discharge, but the location of the scour hole was un- 
affected by the elevation of the channel floor. More important, the water 
surface profile was unaffected by the elevation of the channel floor. The mov- 
able bed was set at the elevation of the lip, at the bucket invert, and below the 
bucket invert. However, the range of tailwater for these particular tests was 
apparently somewhat limited. 


Basic Parameters 


The surge height, hg, is principally dependent upon hg and hp. If hg is very 
large the surge will be “drowned,” as will hp. If h2 is quite small, the surge 
will form a free trajectory (herein called a “free-jet”) and the bucket roller 
will be swept out leaving only a live stream in the bucket—typical perform- 
ance for a “flip bucket.” 

The flow momentum as it leaves the bucket (or enters the surge) includes 
a velocity which is a function of (hy - hp). 


For the overall phenomenon in question, with effects of boundary friction 
and surface tension considered negligible, 


hp (q, hy, hg, R, g). 
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One of the several combinations of dimensionless parameters which satisfy 
the above equation is, 


h 
h 
l fe h, 1 


In a like manner, hg can be included, rather than hp, with 
= (q, hy, hg, R, g), 


for which one of the several combinations of dimensionless parameters 
satisfying this equation is 


As will be shown later, specific service conditions are independent of the 


magnitude of R/h, (used later as hj/R for convenience), and for these 
conditions 


and 


The parameter which includes the discharge can be equated to a Froude num- 
ber but in design computations it is more directly usable as it stands. 


Characteristics of A (1 on 1) Entrance Slope (45° Exit) 


Pertinent results obtained with a (1 on 1) entrance slope are given in 
Figure 3. As h, is progressively reduced, the flow leaving the bucket even- 
tually springs free in the manner of a “flip bucket.” The value of h,/hy 
where this occurred was not affected by the magnitude of hj/R, in any of the 
tests performed throughout the study. This service condition is designated 
“free-jet.” An upper limit for a “free-jet” of hp = 0.2 hg was selected, as 
shown. This limit is not necessarily the exact point at which a ‘free-jet” 
occurred, but for the data obtainable below this limit the jet was always free, 
and above this limit a roller was always obtained. For the lowest points 
plotted below the hp = 0.2 hg line, hp equals the depth of the live stream in 
the bucket. Below approximately hp = 0.2 hg, a very small decrease in hg 
resulted in an abrupt change to a “free-jet.” 

Above the “free-jet” range, with large values of hj/R, a “pulsating surge” 


was obtained for some or all runs of a constant ers series. A “pulsat- 

ing surge” is a pronounced vertical, unsteady motion of the surge, and pro- 

duces an accentuated fluctuating condition where the surge strikes the floor. 
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Inasmuch as a ‘pulsating surge” would undoubtedly be much more aggressive 
from the standpoint of scour, operation under this condition would seldom be 
desirable. 

The only obviously erratic data obtained were those for a “pulsating 
surge.” Ignoring all data for which a “pulsating surge” was obtained, the 
remaining data (for h, greater than 0.2 hg), representing good roller action 
and a normal surge, were plotted in Figure 3. The trends indicated by these 
remaining data were well defined, and the scatter of points could in no way be 
attributed to either a variation in R or hj/R. Referring to Table 1, note that 
most of the series represent a large range of all variables. 

Any scatter of the data points can be attributed to three factors: 


1. The parameter ers was not of exactly the same magnitude for all 
ch 
o l 
series (Table 1). (However, in Figure 3, points lying to the left of a curve 


are generally for a reer value less than that indicated for the curve, 
— 


and points lying to the right of the curve are for a greater parameter 
value). 


2. hp fluctuates somewhat making precise measurements of hp difficult (par- 
ticularly in low range of hp/h,). 


3. Evaluation of dj was not perfect because of the supercritical flow in the up- 


stream channel, and the velocity Vj, was obtained indirectly from a mea- 
surement of dj. 


In Table 1 it may be noted that in some instances a complete range of all 
variables was not accomplished. The lowest value investigated of hj/R = 2 is 
not necessarily a minimum. Lesser values, however, are of little practical 
utility. 

In the last column of Table 1 are shown the ranges of hj/R for which good 
roller action can be obtained. The upper limits should be regarded as mar- 
ginal, inasmuch as a ‘pulsating surge” was obtained for an occasional run in 
a test series for the higher h;/R limits. 

The ranges of h,/R given in the last column of Table 1 are recommended 
for design. The curves appearing in Figure 3 are closely representative of 
the data for good roller action and are recommended for design. (Only about 


half of the data taken appears in Figure 3, the remainder being for test runs 
with a “pulsating surge”). 


Characteristics of A (1 on 2) Entrance Slope (45° Exit) 


In Figure 6 are plotted all data points for which good roller action was ob- 
tained. The curves of Figure 3 for a (1 on 1) entrance are reproduced in 
Figure 6, and appear to be reasonable approximations of the (1 on 2) data. As 
with the (1 on 1) entrance, the scatter of points could in no way be correlated 
with either the variation in R or hj/R and can probably be attributed to the 
same factors as given for the (1 on 1) entrance, above. 

Ranges of test variables, given in Table 2, are adequate but not as broad 
as those investigated for the (1 on 1) entrance. 

A ‘pulsating surge” did not occur in any of the tests with a (1 on 2) en- 
trance. Instead, for higher values of hy/R, a long, high surface jump was 
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formed originating upstream from the P.C. of the bucket and extending well 
downstream from the bucket lip. The bucket did not appear to contribute to 
the jump roller action, since the latter was merely a standing jump originat- 
ing far up in a sloping channel. The maximum hy/R given in the last column 
of Table 2 is for good roller action, but with a sloping channel type of jump 
impending. The sloping channel jump data are not shown in Figure 6; if the 
points were plotted, they would lie to the left of the curves. 

Note that the upper limits for hy/R cited in the right column of Table 2 are 
for good roller action. Higher values of hj/R will not produce either the ad- 
verse operation of a “pulsating surge,” or bucket roller action, but merely a 
sloping channel type of jump described in the preceding paragraph. 


The “free-jet” upper limit of hp = 0.2 hg is the same as for a (1 on 1) en- 
trance slope. 


Effect of Entrance Slope 


Tests were performed on a model of the spillway shown in Figure 2, having 
a (10 on 7) approach and a 45° exit.(2) With rE = 0.055 (design head 
hy 

and q), and hj/R = 3.4, a good fit with the curves of Figure 3 was obtained. 
Equally consistent were the data from the original model tests. For these 
data hj was taken as the difference between headwater pool waterlevel and 

the bucket invert, inasmuch as the original model tests had proven that the 
loss of head due to friction over the spillway was relatively small. 

Special tests were performed on a modified model of a chute spillway(4) 
having a chute slope of 0.104, terminating in a free-trajectory vertical curve 
which becomes tangent to the bucket at an angle of 33° with the horizontal. 
(Laterally flared side walls in the vicinity of the vertical curve of the model 
were replaced by parallel walls for the special tests). These tests were for 
the following parameter values: 


/ 
3/2 
ve hy 
0.038 
0.058 


0.070 


Roller action was fairly good up to an hy,/hy of about 0.4, above which a jump 
was formed upstream from the bucket, originating in the trajectory curve or 
in the chute and extending downstream from the bucket lip. The data (for 
hp/hg > 0.4) was reasonably consistent with the curves of Figure 3 despite 
the complex approach conditions. Since the head loss in the chute was ap- 
preciable, hj was measured just upstream of the vertical trajectory curve. 
Broadly interpreting the somewhat limited investigation of entrance slopes 
described, the curves of Figure 3 should provide a reasonable anticipation of 
performance for approach slopes from (1 on 2) to (10 on 7), provided that 
hj/R is not too large. 
With an approach slope of (1 on 1), the bucket roller becomes “crowded” at 
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h;/R values higher than those recommended and the surge becomes unsteady 
and pulsates. With small slopes, such as (1 on 2), the relatively longer 
bucket roller is drowned at higher hj/R values and the bucket roller and 
surge merge into a single jump apparently indifferent to the presence of the 
bucket. In both cases the bucket is too small, relative to the depth of the 
entering live stream, to govern or control the action of the surge. 


Energy Dissipation Characteristics 


For a non-submerged hydraulic jump in a horizontal channel, a convenient 
appraisal of entrance conditions is obtained using a Froude number, such as 


1 , where V, and D, are the velocity and depth at the entrance to the 
= 

jump, with hy = Dj + v2 /2g. Since hy is measured from the same datum for 
both the or dissipator (bucket invert) and a hydraulic jump in a horizontal 
channel, it may be shown that 


from which it may be seen that a given value of is equivalent to a 


h, 2/2 


1" Equivalent values of Fi for this study are tabulated in 


specific value of F 


Figure 5. 
For a standing hydraulic jump in a rectangular, horizontal channel 


Equation (2) is graphed in Figure 5. The test points plotted in Figure 5 
were obtained by taking values of hg/h, from Figure 3 at the intersections of 
the curve for a given equivalent Froude number with the lines for hp = 0.2 hg, 
0.4 ho, 0.6 hy and 0.8 hog, and by combining these curve values of hg/h, with 
values of hy/D, calculated from Equation (1) for the same Fy: h2/hy . hy/Dj = 
h2/D}. 

For hp = 0.2 hg (the approximate borderline for “free-jet” conditions), the 
bucket dissipator is in close competition with a jump on a horizontal apron. 
For hp = 0.4 hg and Fy as high as about 16, the bucket dissipator has tailwater 


requirements which are fairly close to those for a jump on a horizontal floor. 
Referring again to Figure 3, note that the tailwater requirements for the given 
upper limit for a “free-jet” (hp = 0.2 hg) are not much less than for hp = 0.4 
hg. Therefore, when a “free-jet” condition must be avoided (for example 
proximity of electrical equipment), it is recommended that hp = 0.4 hg be used 
as a limit in design, assuring roller action and efficient energy dissipation. 
Inasmuch as the data for a (1 on 2) entrance slope closely approximated the 
curves of Figure 3, for a (1 on 1) entrance, and since the data from the special 
test with a (10 on 7) entrance were also in agreement, the hg/D, curves of 
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Figure 5 should be nominally representative of this wide range of entrance 
slopes. This is correct, of course, only for true roller action, without either 


a “pulsating surge” or a submerged jump (i.e., provided hy/R is not beyond 
the recommended maximums). 


Surge Heights 


Approximate surge-height characteristics are given in Figure 4 for a (1 on 
1) entrance slope and in Figure 7 for a (1 on 2) entrance slope. This informa- 
tion should be helpful in setting top elevations of spray or training walls. An 
attempt was made to measure an average of typical hg for all runs, but this 
was difficult because of the spray and unsteadiness inherent in the surge. The 
curves are for the same runs as those plotted in Figures 3 and 6, and are an 
average fit to the data. Hence they must be regarded as approximate and not 
necessarily conservative. 

Note that the curves of Figures 4 and 7 differ. Corresponding values of 
hg/hy for the (1 on 2) approach slope are less than those for the (1 on 1) slope. 
(For the (10 on 7) slope they were greater than for the (1 on 1) slope). 

The only variable which is definitely influenced by the entrance slope is 
the peak height of the surge, hs. Comparing Figures 3 and 6, no obvious in- 
fluence on hp or hg is apparent, and therefore the tailwater requirements for 
energy dissipation are virtually identical. Apparently the surge trajectory is 
flattened as the approach slope is lowered, and vice versa, without materially 
affecting any of the other variables. The upper limit of a “free-jet” is desig- 
nated in Figures 4 and 7 as approximately hp = 0.15 hg. 

The horizontal distance downstream from the bucket lip to the peak of the 
surge was measured for all runs, but analysis of the data is inconclusive at 
this time. 

This study did not include the determination of the distance downstream 
from the bucket necessary for expansion of the surge. It is no doubt obvious 
that the surge travels a much greater distance as a “free-jet” than with 
bucket roller action. Maximum scour occurs in the vicinity of the terminus 
of the surge trajectory, and the channel can be raised beyond this area without 


crowding the surge (see Figure 2). Appraisal of the horizontal distance to the 
scour hole is needed to complete the study. 


Illustration of the Use of Study Data for Design 


In Table 3 are presented a hypothetical spillway profile and a tailwater 
rating curve. The accompanying calculations are preliminary since hj, is the 
gross head without any correction for the spillway head loss. For the first 
trial the tailwater is found to be too low to provide roller action, with a “free- 
jet” indicated for all rates of flow: 


Values of 


Available from "(1) First From Figure 3, approximate 


fe i 3/2 Trial", via "(a) Tailwater minimum to avoid a "free-jet": 
1 Curve", of Table 3: 


0.267 0. 
0.042 0.263 0.32 
0.036 0.259 0.30 
0.025 0.237 0.25 
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(Also shown is the fact that the tailwater would likewise be insufficient to 
hold a jump on a horizontal apron. In either case, an end sill could be em- 
ployed to raise hg). 

In the second trial, in which the bucket invert is lowered 7-feet, roller ac- 
tion is assured at all rates of flow. The lowest value of hp/hg happens to be 
0.4, the recommended safe minimum for good roller action. The calculated 
points are plotted in (c), a reproduction of part of Figure 3. 

The approach slope happens to be (10 on 7), and for ERE = 0.04 a 

gh” 
value of hj/R = 5 is selected since in Table 1, for the lesser approach slope 
of (1 on 1), a value of 6 would provide only marginally satisfactory perform- 
ance. The approximate appraisal of hg was obtained using Figure 4. 

Various design combinations of q, hy and hg will yield unique trends when 
plotted on Figure 3; the given sample is but one such combination. A sample 
computation for a high head has not been included because effects resulting 
from air-entrainment would be unknown factors. 


RECOMMENDATIONS AND CONCLUSIONS 


Results of a study limited to a 45° exit angle and a uniform distribution of 
discharge over the channel width have been presented. 


Recommendations for design: 


a) A thickness of bucket lip less than R/10 was used in this study. For 
most applications this thickness is adequate to accommodate 
reinforcing. 


b) It is recommended (on the basis of limited information) that the down- 
stream channel be excavated to the elevation of the bucket invert, at 
least in the region adjacent to the bucket lip. 


c) The curves of Figure 3 (reproduced in Figure 6) should be usable for 
approach slopes from (1 on 2) to (10 on 7). The curves are representa- 
tive of good roller action up to specified maximum values of hj/R. 
“Free-jet” performance can be expected whenever hp is less than 0.2 
hg, in general. The maximum values of hj/R recommended in Table 1 
indicate probable commencement of a “pulsating surge,” whereas in 
Table 2 the maximums delineate the termination of bucket roller action 
and commencement of a sloping channel type jump apparently uninflu- 
enced by the presence of the bucket. The magnitude of these upper 
limits of hj/R were almost identical, but there is no way to tell whether 
a ‘pulsating surge” would or would not occur at these same limits for 
an intermediate approach slope. 


d) The absence of any influence on good roller action due to the size of the 
bucket radius has been demonstrated. An examination of Figures 3 and 
5 will show that a wide range of service conditions has been studied. 
The wide coverage of test variables listed in Tables 1 and 2 sustain the 
direct value of the results in design. 


e) Within the restrictions of a 45° exit angle and a uniformly distributed 
discharge only one major item remains to be investigated: the horizon- 
tal distance from the lip of the bucket to the scour hole (or the 
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termination of the surge expansion). This item is of importance only 
in guaranteeing an adequately large basin. The horizontal natural chan- 
nel requirement in Figure 2 (which is to scale) is only about 2.4 hg; a 


sing apron in place of the bucket in this instance would be much 
longer. 


The required natural channel length could best be investigated using a 
movable bed model. With the presentation of this study, if sufficient interest 
is exhibited by the profession in this particular item, the sponsors of the cur- 
rent study are agreeable to a future cooperative extension of the program. 
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TABLE 


June, 1957 


RANGES OF VARIABLES INVESTIGATED - (1 ON 1) ENTRANCE 


AND 45° EXIT 


Range of Range of Range of Range of 
q hy, h2, 
Ve feet efs./ft feet 


0.089 2.3 
-0.091 


Test Range Range of hj/R 
of for Good 
h,/R Roller Action 


| 
Radii, 
R used, 7 
feet 
0.35 
0.60 
1.20 
0.058 32:0 09+ 2.32 2.8 2-6 
-0.061 0.35 ay 
0.60 
1.20 
0.039 2.4-3.5 0.8-1.5 0.7-2.1 0.30 2.0 11.2 2-6 
-0.043 0.35 a 
0.60 
1.20 
0.025 2.3- 7.2 0.5 2.9 0.6-3.8 0.30 1.9 - 11.4 3-6 
-0.029 0.35 
0.60 
1.20 
0.012 3.4-6.2 0.5-1.1 0.6-2.2 0.30 2.9 - 11.6 3-8 | 
-0.013 0.60 4 
1.20 ‘ihe 
0.010 4.1- 7.7 0.5-1.3 0.6-3.7 6.4-7.2 %-6-7-2 
-0.011 1.20 ul 
‘ 0.005 7.2 0.5 0.7 - 3.0 1.20 6.0 2-6-2 iu 
a 
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TABLE 2 


RANGES OF VARIABLES INVESTIGATED - (1 ON 2) ENTRANCE 


AND 45° EXIT 


Range of Radii, Test Range Range of hj/R 
q> R used, of for Good 
ote./ft. feet h,/R Roller Action 


2.3 - 2.4 . ° : 2. 9.7 


Range of Range of 
q hj, 
0.085 2.8 - 2.9 
-0.092 0.4 
1.2 
0.059 2.6-2.7 1.4-1.5 0.9- 2.0 0.3 2.2 - 8.9 2 - 6 
0.6 
1.2 
0.040 2.4 - 2.5 0.9 0.8- 2.0 0.3 2.0 - 8.4 
-0.042 0.4 
1,2 
pe 0.025 2.4 0.5 0.6- 2.0 0.3 2.0 - 8.0 3 - 6 hee 
-0,027 0.4 
0.6 
0.013 2.3 0.3 O.4-1.9 1.2 1.9 2-2 
i 
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TABLE 3 
SAMPLE COMPUTATIONS—PRELIMINARY DESIGN 


(a) Tailwater Curve 
(1) FIRST TRIAL: 


he 
ft 


5910.025 13.3 
21/0.010] 8.8 


Referring to Fig. 3, the 
above values of h,/h, are all 


to the left of the q 


curves. Therefore, free-jet 


50 100 127 action for all rates of flow. 
q, cfs/ft (At maximum q, TW is insuffic- 
ient to obtain a hydraulic 
jump ina 9-7 channel; 
(b) Spillway Profile h./D, = 7-8 but h,/D, = 10 is 
needed - Fig. 5) 


(2) SECOND TRIAL: 
(Try lowering bucket 7') 


The results are plotted 
in (c). Maximum q near hy, = 
O.4hz and equivalent F, = 8.2. 
Referring to Fig. 5, h2/D, 
for bucket same as for hydraulic 
jump with horiz. floor. 
With h,/R e 5, R = 13.\! 
From Fig. maximum h, /h, 
is approx. = 0.17. Therefore, 
hs. 31.5". 


02 
ho/h, 


0.3 0.4 


223 
H,|~h h 
TW 10] 60 127 0.048] 16.0] 0.267 
EL 9} 59 | 108] 0.042] 15.5] 0.263 
8 58 fa’ 2 0.259 
6| 56 0.237 | 
210 3} 53 0.166 a 
| 
: 
~h, Hs |~h, he |he/h 
rt| ft n™| ct | 
(I 
10/67 | 127/0.041] 23.0/0.343 
ho 9/66 | 108/0.036| 22.5/0.341 q 
207 8/65 | 91]0.031|22.0]0.339 4 
(2) | 21/0.008/15.8]0.263 
(c) From Figure 3 ad 
A. 
0.3 
© 
| 
‘ 0.1 0.2 
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FIGURE 2 
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FIGURE 7 


LL 
portals 
| 
q 

| 

| 

| | A 

| ay. if 

a 
| 
0” 10° 20° 306° 40° 50° 60° 70° 80°96 4 


Journal of the 
HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 


MEASUREMENT OF SEDIMENTATION IN TVA RESERVOIRS*# 


E. H. McCain,! A.M. ASCE 
(Proc. Paper 1277) 


INTRODUCTION 


Investigations of reservoir sedimentation are of significant importance 
wherever storage reservoirs are an integral part of water conservation and 
utilization projects. The broad objective of any reservoir sediment investi- 
gation is to provide information upon which to base estimates of the number 
of years that will be required, first, for sufficient sedimentation to occur to 
interfere with overations, second, for sedimentation to fill the reservoir to 
an elevation w ere its useful life will be ended and, third, to determine cur- 
rently the actu 1l volume of the reservoir for water control operations. 

The TVA is presently investigating sedimentation in 30 reservoirs in the 
Valley, comprising a total of 1150 active sediment ranges. The magnitude of 
these investigations is so great as to warrant special consideration of the 
equipment used in making the investigations. TVA has developed or assem- 
bled commercial equipment which has materially increased the accuracy of 
the hydrographic investigations as well as the efficiency of the fieldwork, 
thereby, substantially reducing the cost of the investigations. 


Equipment 


The equipment used consists essentially of the following: 


1. Truck and trailer 

2. Boat and motors 

3. Two-way radios 

4. Distance measuring equipment 
5. Echo sounding equipment 


The equipment is designed to provide complete mobility and permit the 
entire outfit to move readily from a reservoir in one part of an area to that 


Note: Discussion open until November 1, 1957. Paper 1277 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
83, No. HY 3, June, 1957. 


a. Presented at Knoxville Convention, American Society of Civil Engineers, 
Knoxville, Tenn., June 8, 1956. 

1. Sedimentation Specialist, Hydr. Data Branch, Tennessee Valley Authority, 

Knoxville, Tenn. 
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in another, perhaps several hundred miles distant, in a day’s time. Figure 1 
shows the truck, trailer, and boat ready for movement. The truck is a 3/4- 
ton power wagon equipped with four-wheel drive. The four-wheel drive is an 
important feature in launching the boat into a reservoir and removing it since 
many reservoirs do not have improved ramps. The truck is also equipped 
with power brakes for safety and a power winch on the front of the truck that 
is used to assist with removing the boat from difficult places. The trailer is 
designed especially for moving boats. It has four wheels, two in tandem, 
which are mounted with knee-action devices. This permits the bed to be 
lowered and facilitates the launching and loading of the boat in relatively 
shallow water. All four wheels are equipped with electric brakes which are 
operated from the truck cab. 

The boat is of wood construction, 20 feet long and approximately 5-1/2 feet 
wide. It is a flat-bottom boat with a scow bow which provides for stable 
operation. Its operational weight is approximately 2 tons and it draws about 
one foot of water. It is powered by 2 Scott-Atwater 33 horsepower outboard 
motors. The motors are equipped with electric starters and have a forward, 
neutral, and reverse gear. A unique feature of these motors is the Bail-A- 
Matic device which provides for automatic bailing of the boat when the motors 
are operating. When the motors are operating under full power, they will 
move the boat at a speed of from 12 to 14 miles per hour. 

Good communication is very essential for the efficient operation of a sedi- 
ment investigation party. To provide communication, the TVA party uses 4 
two-way FM radio telephones operating on a frequency of 166.825 MC. the 
truck and the boat have permanently installed General Electric mobile com- 
binations with a carrier power of approximately 8 watts. This provides a 
range between these two radio telephones of from 10 to 15 miles. A portable 
radio telephone is used by the instrument man on the bank and one is main- 
tained as a standby. These radio telephones are the Hallicrafters “Littlefone” 
with an output power of one watt. They have a line-of-sight range of about 5 
miles and are powered by self-contained rechargeable wet cell batteries. 

The total weight of the portable unit is 14 pounds. 

Figure 2 is a view of the interior of the boat showing the equipment. The 
view is from the right stern of the boat looking toward the left bow. The long 
metal box in the right rear of the boat houses most of the power supply and 
controls. There are four 6-volt batteries connected to supply 12-volt power 
for operating the echo recorder, radio, and starting the motor generator. 
They also supply 6-volt power for starting the outboard motors. In the same 
box directly to the rear of the batteries is the battery charger. The battery 
charger operates off of 110-volt AC current from the motor generator and 
charges the battery with 12-volt DC current. To the rear of the battery 
charger, are the controls with a control panel at the rear end of the box. The 
controls are for the radio including a speaker, microphone outlet, and an out- 
let for a headset, starter controls for the outboard motors, and controls for 
the battery charger. Directly opposite the battery box on the left side of the 
boat is the motor generator. The motor generator supplies 750 watts of 110- 
volt AC power which is used for charging the batteries and operating the 1/2 
horsepower 110-volt AC electric motor on the distance measuring device. 
Between the motor generator and the battery box is the transducer well. The 
well houses the echo sounding recorder transducer which transmits the sound 
pressure wave and receives the returning echo. Forward of the transducer 

‘well and underneath the operating platform is the fuel supply. The fuel supply 
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for the outboard motors consists of six 6-gallon cans of gasoline which is 
sufficient for one-day operation. The gasoline cans are connected by short 
rubber hoses to copper tubes which run along each side of the boat back to 
the motors. The fuel supply for the motor generator is a two-gallon can con- 
nected directly by copper tubing to the motor generator. 

Figure 3 shows the distance measuring device with the motor generator. 
The measuring device consists of a steel frame on which is mounted a stand- 
ard automobile wheel, rewind motor, tensiometer, fixed circumference mea- 
suring wheel, and the necessary controls. The automobile wheel is used as a 
drum for 6000 feet of .039-inch diameter music wire. The 1/2 horsepower 
110-volt AC rewind motor is mounted on a movable frame and is belt- 
connected through a jack shaft with the drum. By lifting the motor frame by 
a lever fastened to the top of the main frame, the drive belt will loosen, 
thereby giving a clutch between the motor and the drum. The tensiometer 
consists of a sheave mounted in a sliding bracket with a short rod fastened on 
one end of the bracket which actuates the plunger in a small pressure cylin- 
der taken from an automobile brake. The pressure in the cylinder is trans- 
mitted with fluid through a tube to a pressure gage. 

Figure 4 shows the fixed 2-foot circumference measuring wheel which is 
geared to a Veeder counter which indicates the distance. The opposite end of 
the Beeder counter shaft is geared to a cam which actuates a microswitch 
which automatically causes a fix mark on the echo sounder chart at 50-foot 
intervals of distance along the range. There are two controls for the brake 
in the automobile wheel drum. The lever near the counter operates the me- 
chanical brake and the wheel near the end of the frame operates the hydraulic 
brake. In operation the wire on the drum is paid out through the sheave on 
the tensionmeter and over and around the fixed circumference measuring 
wheel to the bank. The tension on the music wire, while the boat is moving 
along the range, is controlled by the brake on the automobile wheel and the 
distance is recorded automatically on the echo sounder chart or can be read 
directly on the Veeder counter. The operator has a table fastened to the 
counter from which he determines the proper tension to maintain on the wire 
at various distances by observing the pressure gage which is mounted near 
the counter. When sounding of the range is completed the wire is released 
from the bank and rewound on the drum at the rate of approximately 600 feet 
per minute. A foot-operated guide is provided to distribute the wire on the 
drum as it is being rewound. 

The distance device was rated in the field between known distances at 
1000-foot intervals up to 6000 feet. It was rated to read correct distance up 
to 2000 feet by balancing the error due to sag and the stretch due to tension. 
Because of excessive sag it was necessary to increase the tension for dis- 
tances between 2000 and 6000 feet, therefore, distances in this range must be 
corrected from a curve. The accuracy of the equipment has consistently 
proven to be within one foot up to 2000 feet and within 2 feet from 2000 to 
6000 feet. This is well within the limits of accuracy that we must consider. 

The echo sounding equipment used is the Bludworth model ES-123 Depth 
Recorder which is designed especially to record relatively shallow depths of 
water automatically and with a high degree of accuracy. It is rugged, light- 


weight and compact, simple in operation, requiring a minimum of maintenance 


and will stand considerable exposure to the elements. Depths of water from 3 
to 180 feet may be permanently recorded in three ranges (0 to 60, 60 to 120, 
and 120 to 180 feet) at the rate of 300 soundings per minute so that abrupt 
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changes in depths of water are readily observed with the usual boat speed of 
about 3 miles per hour. The equipment operates automatically when the 
power switch is turned on, the only attention required being selection of the 
depth range, an occasional simple check of the battery voltage, and a check 
of the motor speed. Total power required is 12 volts DC at 14 amperes, or 
168 watts. 

Figure 5 is a view of the recorder case showing the operating controls. 
The controls consist of a volt meter, motor speed indicator, power switch, 
fix button, range selector, and a gain control. 

Figure 6 is an interior view of the echo recorder showing the electronic 
circuits and chart propulsion mechanism which generate the electrical im- 
pulses and receive the reflected pulse, accurately time the interval between 
the generation of the pulse and the reception of the reflected pulse, translate 
the interval into units of depth in feet and record the depth permanently on the 
recorder chart. While the paper moves normally at the rate of 1-1/2 inches 
per minute, the drive motor speed may be adjusted t5 percent from normal, 
to compensate for differences in the temperature and salinity of the water. 
An adjustable draft lever, which may be locked into position, is provided to 
compensate for elevation changes in the water surface and also that record- 
ings of depths below the surface or below the bottom of the boat may be made. 
The outboard transducer housed in the well in the center of the boat contains 
a transmitting and receiving magnetostriction type unit in a single cast alumi- 


num housing. It is equipped with cords and plugs for connection to the re- 
corder. 


Field Investigations 


Plate 1 shows the layout of sediment ranges in Douglas Reservoir. This 
layout is typical for the storage reservoirs in the Valley. The ranges should 
be located giving consideration to important local tributaries and the drainage 
area of these tributaries and their probable sediment characteristics. The 
width of the reservoir should be given consideration. In storage reservoirs 
subject to considerable drawdown, closer spacing of the ranges within the 
drawdown reach is desirable as the deposition and movement of sediment 
throughout this range is greater than in other parts of the reservoir. For 
example, in Douglas Reservoir which is subject to drawdown, ranges are 
spaced about two miles apart in the lower and extreme upper ends of the 
reservoir and 1 to 1-1/2 miles in between. 

Prior to the filling of the reservoir the selected ranges are permanently 
marked at each end and a cross section taken along the range by conventional 
methods. On some of the main-stream reservoirs which are relatively shal- 
low and have comparatively smooth bottoms, the ranges were marked and the 
original cross sections were taken by hydrographic methods soon after the 
reservoir was filled. Utilizing the hydrographic equipment available, this 
method was more economical in establishing the original cross section where 
the terrain was not so rugged as on the storage reservoirs. Subsequent in- 
vestigations on all reservoirs are carried out by hydrographic methods 
utilizing the equipment described. 

The field party required to carry out the sediment investigations consists 
of 4 men—the party chief who supervises the work and operates the distance 
measuring device, the echo sounding recorder operator, a transitman who 
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keeps the boat on a range line, and the motorboat operator. In sounding a 
range a light aluminum boat powered by a 16-horsepower motor is used in 
conjunction with the sounding boat. This boat is for the use of the transitman 
and saves numerous trips across the reservoir with the heavier boat. The 
transitman goes ahead in the light boat to the range to be sounded and clears 
and flags one of the markers. He then meets the sounding boat at the oppo- 
site marker. The transitman sets the instrument up over the marker and 
sights to the marker across the reservoir. An iron pin is then set on the 
range line near the water edge and a tape distance taken to it. A distance is 
taken to the water edge and the wire from the distance measuring device is 
fastened to the iron pin. The echo sounding recorder is started and the boat 
proceeds along the range line at a speed of approximately 3 miles per hour. 
The transitman, located on the bank where the wire is fastened, keeps the 
boat on the range line by means of radio communication to the boat operator. 
During normal weather conditions the boat stays within 2 feet of the range 
line. When the boat reaches the opposite bank, the transitman is informed by 
radio and he releases his end of the wire. As the wire is being rewound on 
the drum the transitman is proceeding to the next range where the process is 
repeated. 

The time required to sound a reservoir, such as Douglas where the total 
length of reservoir including tributaries is about 80 miles, the number of 
ranges is 51 and the ranges vary in length from a few hundred feet to about 
6000 feet, is about one week for the field party. In Kentucky Reservoir on the 
main stream, where the total number of ranges is less (only 38), the total 
length of reservoir is about 230 miles with ranges varying in length up to 
10,000, the time required to carry out the investigation is about 3 weeks. 
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AUTOMATIC VHF RADIO TELEMETERING 
OF HYDROLOGIC DATA! 


James A, Dale2 
(Proc. Paper 1278) 


SYNOPSIS 


With the ever-increasing demand for water in this country for power, irri- 
gation, and industrial use, it has become necessary to ensure that existing 
water-supply facilities are operated as efficiently as possible. The efficient 
operation of reservoirs deriving their water from natural flow is dependent 
on a knowledge of the quantity and distribution of water within the drainage 
area. Although this information can usually be obtained by telephone from 
observers hired to read a rainfall or stream level gage, it is often desirable 
to obtain hydrologic data from isolated areas not served by telephone lines. 
The use of radio as a means of making this information available from these 
areas naturally suggests itself. 

As early as 1937 the Tennessee Valley Authority started work on a radio 
system to transmit hydrologic data from isolated areas. By 1949 the system 
consisted of some 20 rain gages and 20 stream gages. Information from these 
gaging stations was received at nine area offices scattered throughout the 
Valley. These reports, along with reports received from observers in the 
area, were tabulated and telephoned to the Forecasting Section in Knoxville 
daily. 

This system had three drawbacks: 


1) Operation was not completely automatic. 
2) The transmitting frequencies used, which were between 2 and 3 mega- 
cycles, were subject to fading and atmospheric disturbances. 


3) Because of the rapid advancement in the electronic field, by 1949 the 
system was obsolete. 


Also, the Federal Communications Commission had served notice that at 


Note: Discussion open until November 1, 1957. Paper 1278 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
83, No. HY 3, June, 1957. 


1. Presented at Knoxville Convention, American Society of Civil Engineers, 
Knoxville, Tenn., June 8, 1956, 


2. Elec. Engr., Hydr. Operations and Tests Section, Hydr. Data Branch, 
Tennessee Valley Authority, Norris, Tenn. 
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some future date operation of this type would have to be discontinued on these 
frequencies. 

In 1950 TVA decided to replace the old radio gaging system. Design ofa 
new system was started in 1950 by TVA engineers, and installation was com- 
pleted in 1955, 


System Requirements 


The system was designed around the following basic requirements: 


1) Location of stream gage and rain gage stations would be essentially the 
same as in the old system. 

2) Receiving stations would be located at seven area offices. Each office 
would receive only the gaging stations within its area. 

3) Each gaging station was to report every two hours and record auto- 
matically at the receiving station. 

4) Radio frequencies for transmission were to be selected from the 16 fre- 
quencies between 169 and 172 megacycles (Mc) allocated by the Federal Com- 
munications Commission for the transmission of hydrologic data. 


System Layout 


The radio frequencies available for use in the system have essentially 
line-of-sight transmission characteristics, and since gaging station and re- 
ceiver locations were already established it was necessary to provide relay 
stations to make possible the reception of gages at the receivers. Tentative 
relay or repeater station sites were selected with the aid of topographic 
maps, but final locations were established only after signal strengths between 
gaging stations, repeater, and receiver were measured with portable and 
mobile radio equipment. It was found that a total of 13 repeaters would be 
required, which is approximately one repeater for every three gaging sta- 
tions. The locations of all stations in the final system are shown on the map 
in Figure 1. 


System Operation 


In order to explain the operation of the system, a typical stream gage- 
repeater-receiver link will be considered. A clock-operated switch at the 
stream gage initiates the operation of the gaging equipment. An unmodulated 
carrier from the gage transmitter is transmitted and received at a repeater 
station, turning on the repeater transmitter. A repeater carrier is trans- 
mitted to the receiving location to complete the circuit from the gage. Con- 
version of stream level at the gaging station into a numerical sequence of 
events suitable for telemetering is accomplished by a digital converter or 
keyer. This device translates the river stage into three groups of from 0 to 
9 contact closures corresponding to the stage in tens of feet, feet, and tenths 
of a foot. A fourth group of from 0 to 9 contact operations precedes these 
three groups and is preset to identify the gaging station. These four groups 
of contact closures key an audio tone generator, which in turn 
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frequency-modulates the gage transmitter. This audio tone frequency is 
designated as frequency A. The digital converter unit also causes a second 
set of contacts to close between each of the four A tone groups and immedi- 
ately before and after the first and last groups. These contacts key an audio 
tone generator on frequency B, which also modulates the transmitter. The 
audio tone pulses are received at the receiving station via the repeater. 
After demodulation they pass through frequency-selective filter networks and 
actuate the A and B tone relays. The A tone relay contacts are connected in 
turn through a stepping switch to each of four solenoid-operated print wheels 
in the recorder. The B tone relay contacts actuate the step coil of the step- 
ping switch. Thus it may be seen that the four groups of A tone pulses are 
each read out in turn by the four print wheels, and the resulting record gives 
the station number and the stage. Included in the print-out section of the re- 
corder is a date and time stamp which prints opposite the stage recorder. 

The chart in Figure 2 shows the time sequence of various operations oc- 
curring during one gage transmission. It will be noted that information from 
the gaging station is sent out twice during one transmission. This is bene- 
ficial in verifying the accuracy of the recorded data. 

The example illustrated in Figure 2 is station 4 broadcasting a stage of 
10.8 feet. At zero time the clock switch initiates equipment operation and 
the transmitter carrier is turned on. Six seconds later the B tone keyer con- 
tacts cause the carrier to be modulated by the B tone. Duration of this 
modulation is about 1/4 second. At the receiver this modulation pulse causes 
the B tone relay to operate, which in turn steps the stepping switch in the re- 
corder off the normal “no contact” position to position 1. Immediately after 
the occurrence of the first B tone pulse the keyer causes a series of four A 
tone pulses to be transmitted. These are received at the receiver and 
actuate the A tone relay. The contacts of this relay are connected through 
the position 1 contacts of the recorder step switch to the solenoid of the first 
print wheel, and the print wheel is stepped around to position 4. The second 
B tone impulse is transmitted at 0 + 12 seconds and causes the recorder step 
switch to move to position 2. The second print wheel is then stepped around 
to position 1 by the following single A tone pulse. At 0 + 18 seconds the third 
B tone pulse is transmitted. Following this there are no A tone pulses so the 
third print wheel is not moved. After the fourth B tone pulse at 0 + 24 sec- 
onds a series of eight A pulses steps up the fourth print wheel to position 8. 
The fifth B tone pulse at 0 + 30 seconds causes the recorder to print and 
clear. The recording cycle is then repeated, and after the tenth B tone pulse 
the transmitter is turned off. The resulting printed record is 4108 indicating 
that station 4 has a stage of 10.8 feet. A typical recorder tape showing the 
record of ten stations in a two-hour period is illustrated in Figure 3. 


Power for Station Operation 


With the exception of two rain gages, all stations are served by commer- 
cial power. In the event of power failure, emergency power is automatically 
provided at repeater stations by a gasoline-driven generator having a fuel 
capacity sufficient for five days of operation. When the repeater is on 
emergency power a tone is sent out from the repeater at the end of each 
transmission. This tone operates a relay at the receiver which in turn lights 
a light on the receiver panel. The light stays on until released manually but 
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Figure 3.--RECORDER TAPE SHOWING RECORD FOR TEN STATIONS 


will come on again on the next repeater broadcast if the main power is still 
off. At gaging stations, ten-day emergency power operation is provided by a 
60-ampere-hour battery. Normally a gaging station transmits the required 
information twice during each broadcast. If the gage is on emergency power, 
the second half of the broadcast appears on the recorder tape as four zeros. 
No emergency power system is employed at the receiving stations. 

Two rain gage stations in the system are in such remote locations that 4 
the cost of bringing in commercial power is prohibitive. One of these, he 
Clingmans Dome, is near the middle of the Great Smoky Mountains National 1 
Park. The other station, Haw Knob, is in an area populated only by wild a” 
boars. 

The equipment at these two stations is basically the same as that used at 
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other stations. Power is provided by two 100-ampere 6-volt heavy-duty 
truck batteries. These batteries will run the station approximately 25 days. 
A gasoline-driven generator is provided to charge the batteries. These sta- 
tions are normally visited every two weeks, and at this time a 3-hour supply 
of gasoline is added to the generator. The generator is then left running and 
it will, of course, stop when the fuel supply is exhausted. With additional 
battery capacity and a longer charging time it should be possible to get at 
least 60 days of gage operation. 


Equipment 


In the interest of economy, the system was designed to include standard 
available components, such as radio transmitters and receivers, when possi- 
ble. Particular attention was given to the standardization of parts and sub- 


assemblies to facilitate servicing and minimize the stock of spare parts 
required. 


Stream Gage 


The components of the stream gaging stations are listed below and shown 
in Figure 4. 


1) Keyer or digital converter 
2) Clock-operated switches 
3) Transmitter assembly 

4) Emergency power battery 


The keyer was developed at TVA’s Hydraulic Laboratory and is manu- 
factured by the Instruments Corporation. Two types are available, one hav- 
ing a range of 00.00 to $9.99 feet and one with a range of 00.0 to 99.9 feet. 

The gaging station is turned on once every two hours by a clock-operated 
switch. Two types of clock switches are used at each gaging station served 
by commercial power. One is operated by a synchronous electric motor and 
the other by an electrically wound spring-driven escapement. The more ac- 
curate synchronous clock normally operates the station. When a power 
outage occurs the spring-driven clock is automatically switched in to take 
over the control function. 


The transmitter assembly consists of six panels which are, from top to 
bottom, identified as follows: 


Tone generator panel 
Relay panel 

Battery charging panel 
Power supply 
Transmitter 
Converter panel 


The audio tones are produced by two General Electric tone generators 
mounted on the tone generator panel. These units are plug-in, making them 
easily replaceable. 


Relays necessary for keyer operation and control functions are located on 
the relay panel. 


The charging panel provides a trickle charge to the emergency power 
battery. 
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Figure 4.--STREAM GAGE KEYER, TRANSMITTER, CLOCK SWITCH, AND BATTERIES 


Sample of gaging station equipment as exhibited at Knoxville Convention 
of ASCE. Gage well is simulated by a transparent water jar to show the 
float. 
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All vacuum tube plate and filament voltages are supplied by the power sup- 
ply panel. 

The radio transmitter is a standard General Electric 6-watt unit. 

On the bottom panel is mounted a converter for converting the 12-volt 
battery power to 115 volts AC. This unit is automatically switched on when 
the commercial power fails and supplies 115-volt AC to all units. 

The emergency power batteries are of the nickel-cadmium type and were 
chosen because of their ruggedness and ability to stand an overcharge. A 
typical stream gaging station is shown in Figure 5. 


Rain Gage 


The rain gage station equipment is identical to that used at the stream 
gage except for the method used to position the keyer. Rainfall is collected 
in a standard Instruments Corporation 20-inch capacity weighing type rain 
gage. The weighing mechanism operates a low-torque potentiometer which in 
turn positions the keyer through a servo system. The rain gage station 
broadcasts the catch in the collector to the nearest tenth of an inch. The 
rain gage collector, keyer, and servo system are shown in Figure 6 and a 
typical rain gaging station is shown in Figure 7. 


Repeater 


The repeater consists of a standard General Electric transmitter and re- 
ceiver, shown in Figure 8. The bottom panel was modified by TVA to include 
a tone generator for indication of emergency power operation. Accessories 
at the repeater station include a 400-watt gasoline-driven generator for 
emergency power, automatic line transfer switch, thermostatically controlled 
heater and ventilating fan, and voltage regulator. A typical repeater station 
installation is shown in Figures 9 and 10. 


Receiving Station 


The receiving station equipment in Figure 11 consists of three units: 
1) Receiver 

2) Recorder 

3) Tone switch unit 


The receiver is a standard General Electric unit with no modification. 
The printing recorder is manufactured by Streeter-Amet and is modified by 
TVA by the addition of a stepping switch and two relay sub-assemblies. 
Suitable mounting and auxiliary equipment are provided by TVA. 


Antennas 


Transmitting antennas used at both repeater and gaging station are of the 
corner reflector type having a power gain of 5 and a 60-degree beam width. 
Non-directional ground plane antennas are used on all receivers. 


System Versatility 


It was mentioned previously that the rain gage collector operated a low- 
torque potentiometer which in turn positioned the rain gage keyer. Any type 
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Figure 5.--STREAM GAGING STATION WITH RADIO TRANSMITTER 


Panel of transmitter assembly is swung open, showing equipment on back 
of panelboard. Antenna is the corner reflector type. 
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Figure 6.--RAIN GAGE KEYER, COLLECTOR, AND SERVO SYSTEM 


The servo system, shown on the exhibition table, replaces the recording 
mechanism in a standard recording gage. 
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Figure 11.--RECEIVER AND PRINTING RECORDER 
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of transducer capable of positioning this potentiometer may be used as a 
sensing element, thus making it possible to transmit other types of informa- 
tion such as pressure, temperature, and relative humidity. 

With no modifications to the repeater stations it is also possible to incor- 
porate one or more 2-way radio telephone circuits. For example, by adding 
a transmitter to the receiver station at Bryson City, 2-way communication 
could be carried on between Bryson City and cars equipped with standard 
mobile radio units. The distance of communication is determined by the 
coverage of the repeater. 


System Performance and Cost 


Probably the two questions most often asked about the system are: 
(1) How much did it cost? and (2) How reliable is it? 
Our actual equipment cost for the four types of stations is as follows: 


Rain Gage - $2,400 
Stream Gage - $1,800 
Repeater - $2,600 
Receiver - $2,100 


All preliminary design costs are included in these figures, as well as 
costs of construction, assembly, and test. Overhead costs are not included. 
Stream gages are installed in the regular stream gage houses so that 
there is no additional cost for housing. Rain gages include the cost of the 
equipment house, collector and servo system, which combined result in a 


higher cost for the rain gage over the stream gage. 

As a performance check, a record of monthly reporting efficiencies of all 
gaging stations is kept. Efficiencies are based on one report every two hours 
from each station or 12 reports a day. The average system efficiency for 
February 1956 was 96.0 percent; for March, 96.8 percent; and for April, 

97.6 percent. 
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Discussion of 
“THE PROBLEM OF RESERVOIR CAPACITY FOR LONG-TERM STORAGE” 


by A. Fathy and Aly S. Shukry 
(Proc. Paper 1082) 


HAROLD EDWIN HURST.*—There are two principal problems involved in 
connection with reservoir capacity for long-term storage. The first is what 
are the characteristics of natural phenomena such as river-discharges on 
which the amount of storage for different conditions depends, and what are the 
relations between them. The second is to devise methods for the regulation 
of the storage when it has been provided. An answer to the first was given by 
the writer in a paper to ASCE, (1) while the second was discussed by him ina 
paper presented to the Institution of Civil Engineers London in March 1956.7 
However a fresh approach to a problem is often valuable. It is unfortunate 
that the second paper was not available to the authors. 

The relations found by the writer are given in equation 2 for random 
events, and in equations 5, 6 and 8 for some natural phenomena. Eq. 2 was 
found theoretically and was later verified by W. Feller using a different 
method.® In the case of natural phenomena a mathematical derivation is not 
at present possible, though some work now in progress may perhaps lead to 
one. In order therefore to arrive at equation 5 seventy-five phenomena were 
analysed and for this purpose their observations were divided into 690 sets in 
which N ranged from 30 to 2000. These were later extended by the addition of 
the discharges of 21 rivers. The phenomena were divided into classes of 
similar members, rainfalls, river statistics, varves etc. For any particular 
class the results from the individual sets of observations were arranged in 
groups according to values of N. For each group mean values of log R/o and 
log N were calculated and plotted against each other. This is shown in Fig. 14 
which is fig. 4 of the A.S.C.E. paper(1) which is here reproduced. From this 
it is clear that the points fall closely on straight lines, which are represented 


by 
log = K log 


a. Scientific Consultant to the Ministry of Public Works, Egypt; formerly 
Director- Gen., Physical Dept., Ministry of Public Works, Egypt. 

7. Methods of using Long-Term Storage in Reservoirs. Hurst. Proceedings 
of the Institution of Civil Engineers. London. Part 1. vol. 5. pp. 519-590 
September 1956. 

. The asymptotic distribution of the range of sums of independent random 
variables. Annals of Math. Statistics. Vol. 22. 1951, p. 427. 
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The mean values of K for each of the classes of phenomena varied from 0.69 
to 0.80. The 690 individual values of K are normally distributed with a mean 
value 0.73 and a standard deviation 0.092. The reasons for adopting the sim- 
pler formula with one parameter K instead of two are given in the discussion 
on the paper. The authors object to this on the grounds that “the consistency 
apparently gained through that extrapolation is the product of the procedure 
itself and not of any uniformity in the character of the various phenomena 
examined.” However the uniformity is plainly to be seen from fig. 14. It 
should be pointed out that the dependence of R on N implied in the above equa- 
tion is a stochastic one, and not, as the authors seem to think it ought to be, a 
functional one. 

The authors’ new line of approach is to express R (Equation 10) in terms 
of Mg a hypothetical absolute mean, N, the length of the period covered by the 
range of the mass-curve of the phenomenon considered, which they call the 
“critical period,” and 6, and 6¢, where 6,9 and 6¢ are the ratios of the de- 
partures from the absolute mean of the mean of the N available observations 
of the phenomenon and the mean of the Nc observations covering the critical 
period, to the absolute mean. The introduction of an absolute mean, Mg, 
which is described as “the most probable value of the mean for any period,” 
but actually does not exist, only serves to obscure the fact that 


R = Nce(Mo- Me) 


The authors then proceed to adopt a standard value of 1/2 for N-/N and pro- 
pose to take a value of Mj-Mc of 3.5 times the probable error of Mo as a 
maximum which will not often be exceeded. But Mc is the extreme value of 
the mean of N, consecutive observations taken from the sample N, so that 
Mo-Mc¢ may easily exceed 3.5 times the probable departure of the mean Mo. 
Substituting these values in equation a we get 


R=05N35.0670NN = 


The coefficient 1.18 is less than the mean 1.25 found for random events. It 
does not correspond with the 1.66 which results from the complicated opera- 
tions carried out by the authors. It can be shown however that the standard 
deviation of R for random events is 0.270YN.8 If this is multiplied by 2.36, 
as proposed by the authors, we get from equation 2 


R = 1°89 


It must be remembered that this only applies to random events, and that for 
natural events R increases more rapidly than YN (see fig. 14). The mistake 
made by the authors is that they fix Nc as 1/2N and then ignore its variation, 
which is by no means negligible. For example in the case of Aswan discharge, 
1870-1956, the critical period started with 1898 and is still continuing. 

V/ith regard to the storage S required to meet the maximum deficit when 
the draft is less than the mean for the period, for which average relations 
were given by the writer, the authors remark “Thus there are no grounds for 
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the supposition of a fixed relation between S and R.” It is only necessary to 
look at figures 8 and 9 of (1) which represent mean results from 45 different 
natural phenomena to see that there are very good grounds in support of a 
stable mean relation between S and R. 

The preceding relates to the statistical description of natural phenomena, 
and it may be emphasised that there is only one method of investigation, and 
that is the analysis of a large number of phenomena with records covering 
long periods. Instead of this the authors make R a function of 4 quantities 
(equation 10) and then proceed to make hypotheses about all of them and base 
their conclusions on these. 

Incidentally it may be remarked, that the writer’s practice of using cumu- 
lative sums of departures from the mean (or a base near the mean) is much 
more convenient than summing the actual observations, since it is less 
laborious, is largely self checking, and produces a much more compact graph 
when the cumulative sums are plotted. 

As the authors say the decision on the amount of storage to be provided de- 
pends on the special conditions of the project, site, water requirements, cost, 
and, in arid countries, evaporation losses; but other things being equal the 
larger the storage the better. 

Having provided the storage, the very important question of how to use it 
arises. This has been discussed by the writer in (7) where the method 
adopted was to experiment with 50 natural phenomena most of them having 
records covering 100 years or more, taken mainly from river discharges and 
levels, rainfall, temperature and varves. 

The first object was to find a regulation which would make good use of the 
water and minimise the risk of emptying the reservoir. Nine regulations 
were tried on the assumption that the first 30 years records were known from 
which the capacity of the reservoir R was calculated by equation 6. In doing 
this o was increased by 10% to allow for the shortness of the known period, 
and N was taken as 100. The reservoir started half full, and the draft was the 
mean of the proceding 10 years which was reduced or increased on a sliding 
scale by a proportion depending on the amount of water in the reservoir. The 
draft was changed every year. The result of this regulation was that in 14% 
of the cases a sliding scale depending on the content was needed to prevent 
the reservoir from emptying, and in only about 10% of the cases would this 
have made more than a trifling reduction of the draft from the 10 year mean. 
Similar investigations were made for flood protection. 

As the author’s paper is not based on the analysis of data relating to a 
large sample of natural phenomena, but on assumptions which lead to results 
not in accordance with the known characteristics of natural phenomena, the 
writer considers it unnecessary to discuss the applications of these results 
to practical problems. 

The following references relate to random variables. 


9. E. H. Lloyd and A. A. Anis. “On the range of partial sums of a finite num- 
ber of independent normal variates.” Biometrika, vol. 40, June 1953. 


10. A. A. Anis. “The variance of the maximum of partial sums of a finite 
number of independent normal variates.” Biometrika, vol. 42, June 1955. 


11. A. A. Anis. “On the moments of the maximum of partial sums of a finite 
number of independent normal variates.” Biometrika. 
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Discussion of 
“FLUID RESISTANCE TO CYLINDERS IN ACCELERATED MOTION” 


by S. Russell Keim 
(Proc. Paper 1113) 


LOUIS W. WOLF.(19)—Mr. Keim has presented useful information on the 
drag exerted on an immersed object in unsteady motion. There is a regret- 
table lack of information on this topic and more is certainly needed. Unfor- 
tunately, the approach based on a single coefficient cannot be readily extended 
to other modes of motion. Because the total resistance in unsteady flows has 
two distinct parts, one due to velocity and the other to acceleration, it seems 
misleading to combine it into a general coefficient, particularly since valu- 
able information concerning the second part is available. 

The unidirectional motion used provides values of the overall resistance 
for cylinders accelerating with the application of a constant force; this motion 
is not well suited to the measurement of the components of resistance or the 
understanding of the mechanism of the formation of the wake. The velocity- 
dependent drag and acceleration-dependent inertial forces always appear to- 
gether and can be separated only by application of a knowledge of the variation 
of the wake or by causing them to vary in such a manner that they can be dis- 
tinguished. Keulegan and Carpenter(20) and McNown(2!) separated the drag 
and inertial effects for oscillating flows. Because the drag is always in phase 
with the velocity and the inertial force in phase with the acceleration, the ef- 
fects were distinguishable. McNown, using Riabouchinsky’s model for the 
wake, has predicted a relationship between Cp and k for flow past a flat plate. 
Good correlation was achieved between the predicted and experimental values. 

A complication in the description of the resistance is the effect of separa- 
tion and subsequent formation of a wake; the size and shape of the wake affects 
both the coefficient of virtual mass and the coefficient of drag. In the afore- 
mentioned studies of flat plates and cylinders, the coefficient of virtual mass 
was shown to increase above the classical values with the increasing size of 
wake. Values of the coefficient of virtual mass of three times the classical 
value were observed. 

Similarly, the pressure behind any object moving in a fluid is dependent 
upon the wake pattern. For small wakes with a small curvature the pressure 
is considerably below the ambient whereas for large wakes with relatively 
large curvature the pressure is higher. Thus the drag coefficient is large for 
a small wake, and approaches the steady state value for a large wake. Small 
variations of this pattern occur when eddies are shed. Drag coefficients as 
large as ten (five times the value for steady flow) have been observed in ex- 
periments on flat plates. 


19. Teaching Fellow, Univ. of Michigan, Ann Arbor, Mich. 
20. Keulegan, G. H., and Carpenter, L. H., Forces on Cylinders and Plates in 
an Oscillating Fluid. NBS Report 4821, Sept. 1956. 


21. McNown, J. S., Drag in Unsteady Flow; Proc. 9th Int. Cong. Appl. Mech., 
Brussels, 1956 (publ. pend.) 
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The force on any body can be expressed as the sum of the drag CpP A v2 


and the inertial force p Vka due to the virtual mass of the fluid, Cp being the 
coefficient of drag and k being the coefficient of virtual mass. Hence, it can 
be seen that the overall resistance coefficient C, defined by the author, can 
be written as 


C +Fk(43) 


For large accelerations the virtual mass term is dominant, the curves in 
Figure 4 being asymptotic to a line having a unit slope and an intercept value 
of 7k for da =1. For small accelerations the drag term is dominant, 

with the curves asymptotic to a horizontal line representing the steady state 
drag. If the coefficients of drag and of virtual mass can be reported sepa- 
rately, a better understanding of the resistance phenomena would be achieved. 

Values of k indicated by the asymptotes to Keim’s experimental curves are 
consistently less than the values obtained by Stelson and Mavis!1 for cylin- 
ders, and do not vary uniformly with the length to diameter ratio L/d. Stelson 
and Mavis report virtual mass coefficients of .98, .97, and .90 for L/d ratios 
of 30, 15, and 5 whereas the author’s curves indicate .73, .90, and .50, respec- 
tively. A higher rather than a lower value of k would be expected since the 
formation of a wake increases the virtual mass. Because of the extreme 
complexity of the problem, limiting values such as these must be carefully 
assessed. 


11. Stelson, J. M. and Mavis, F. T. Virtual Mass and Acceleration in Fluids, 
Proc. ASCE, Vol. 81, No. 670, 1955. 
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Discussion of 
“FREQUENCY ANALYSES OF STREAMFLOW DATA” 


by David K. Todd 
(Proc. Paper 1166) 


H. ALDEN FOSTER,* M. ASCE.—There has been considerable confusion 
for many years in the use of the term “frequency” in hydrology. As generally 
used, if statistics are arranged so as to show the number of times, or fre- 
quency with which, an event happens in a particular way, the arrangement is 
called a “frequency distribution.” The curve or formula used to describe 
such distributions is called a “frequency curve.” (Ref. 14, pg. 145). The 
author illustrates such a curve correctly in Fig. 3 where he also shows a 
“cumulative frequency curve,” which is the integral of the original frequency 
curve. 

As the frequency curve shows how often a flood of a given magnitude (or 
between certain limiting magnitudes) will occur on the average in a 100-year 
period, the cumulative frequency curve shows how often a flood greater (or 
less) than a certain magnitude will occur according to the same 100-year 
record. This is also an indication of the “probability” of occurrence of such 
a flood, according to the record; therefor it is proper to refer to the cumula- 
tive frequency curve as a “flood-probability curve.” This is also in con- 
formance with standard nomenclature in mathematics where the term 
“Probability Integral” is commonly used. 

It has become rather common practice in discussing the probability of 
occurrence of floods to speak of “flood frequency” as indicating the average 
time interval between successive occurrences of a flood of a given magni- 
tude,—such as once in 100 years. A better term for this purpose would be 
“recurrence interval” or “return period” of such a flood, and to confine the 
term “frequency” to the basic meaning as used by statisticians. 

Mr. Todd points out that a number of different mathematical formulas have 
been proposed for the computation of a theoretical frequency curve to repre- 
sent a particular set of hydrological data; and that computational procedures 
have been developed that make solution of the formulas relatively simple. A 
word of caution may be introduced here. Mechanical simplicity of calculation 
based on complicated formulas is no guarantee of the precision or reliability 
of the results. Statistical methods can never give precise answers, but only 
an indication of the probable values that can justifiably be deduced from the 
recorded data. The writer pointed this out in 1924 (Ref. 14, pg. 172) by the 
statement ‘the theoretical curves should be used only as a guide in the study 
of a given record.” The author indicates the fallacy and danger of trying to 
“estimate” a 10,000-year flood from a record of only 20 years duration. 
Great caution should be used in applying flood-probability analyses to the de- 
sign of the spillway for a dam, even though the structure is of minor im- 
portance. The methods have a proper and useful application, however, in 
certain fields, particularly in estimating Average Annual Flood Benefits for a 


* Prin. Associate, Parsons, Brinckerhoff, Hall & Macdonald, New York, N. Y. 
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flood-control project or in trying to determine suitable premiums for flood- 
insurance policies. (Add ref. 28: H. Alden Foster, 1954; “Flood Insurance,” 
Proc. ASCE, Vol. 80, Sep. No. 483). 

The relative advantages of using the “annual flood method” vs. the “basic 
stage method” or the partial duration series are discussed in Ref. 22; they 
were also explained by the writer in Ref. 20, pg. 54-57. 

The frequency analysis of all flows, involving preparation of a Daily-Flow 
Duration Curve, is on a somewhat different basis from the use of theoretical 
frequency curves for determination of flood probabilities. Such a Duration 
Curve should more properly be considered as a convenient tool for study of 
the data, just as the hydrograph or mass curve is customarily treated. (Ref. 
15) Some attempts have been made, however, to plot flow-duration curves by 
mathematical formulas, notably by Lane and Lei (Ref. 21) who represented 
the daily-flow duration by the logarithmic-probability curve. 
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Discussion of 
“THE ESTIMATION OF THE FREQUENCY OF RARE FLOODS” 


by Benjamin A. Whisler and Charles J. Smith 
(Proc. Paper 1200) 


GORDON R. WILLIAMs, ! M. ASCE.—In spite of persistent research over 
more than 30 years, the problem of estimating the frequency of floods, par- 
ticularly rare floods, is far from solved and may have to wait for many more 
decades of records before a reliable answer can be obtained. 

The authors’ conclusions that the frequency of floods does not follow the 
mathematical theory of probabilities, appears to be sound from a hydrologic 
basis alone. It is well known that there are complex cyclical trends affecting 
flood frequency, beginning with periods as short as a year and extending over 
several decades or more. Most localities are subject to seasonal trends, thus 
introducing a short-term rhythmic effect which has little relation to pure 
chance. 

This seasonal effect appears to cause doubts as to the logic in the use of 
one peak discharge for each month of record. Such data from a dry-season 
record will introduce many values which are peaks in the river hydrograph 
but are in no sense flood flows. In fact, the use of monthly peaks causes the 
method to approach the flow-duration curve in which the discharges compiled 
are often interdependent and do not represent isolated hydrologic events. 

More and more evidence is being obtained that the frequency of rare floods 
cannot be obtained by extrapolating the well-defined mean curve for frequent 
floods. This is another conclusion of the paper with which this writer is in 
accord, but he does have misgivings about the hydrologic soundness of the 
method by which the conclusion is reached. 

There are at least two reasons why there should be two frequency curves, 
one is meteorologic and the other is hydraulic. The meteorologic reason is 
that the great storms causing rare floods are possibly different physical 
phenomenon from the storms causing minor flood rises. For example, rain- 
fall-depth-duration-frequency curves in the West Indies show that the trade- 
wind showers define one set of depth-duration-frequency curves while the 
hurricane storms define another set. Rainfall depths from the former curves 
indicate that the depths from the latter curves are unlikely te occur, even 
though both curves are defined over a 40 to 50 year period. The hurricane 
depth-frequency curves have greater slopes when the depths are plotted as 
ordinates. The floods resulting from these two types of storms also fall into 
two different frequency distributions. It is entirely possible that rainfall 
depths from rare storms inthe United States define their own frequency 
curves and that they are not as rare as the extrapolated frequency distribution 
of lesser storms would indicate. 

The hydraulic reason for the different flood frequency relations mentioned 
above is the effect of channel and valley storage. In the minor flood the chan- 
nel storage may exert a modifying influence that is essentially proportional to 


1. Prof. of Hydr. Eng., Massachusetts Inst. of Technology, Cambridge, Mass. 
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the flood peak or flood volume. The rare flood, on the other hand, occupies 
all surcharge storage in swamps and ponds and covers the river valleys from 
hill to hill. The increment of natural storage per increment of flood volume 
becomes materially less and the flood rolls on with apparently little modifica- 
tion. It is these effects that have led to man-made encroachments in valleys 
with later disastrous results when the rare flood takes place. 

The use of 24-hour flows may result in unsatisfactory floodpeak estimates 
in the smaller river basins of less than about 1000 square miles. In even 
smaller basins of less than 100 square miles the duration of the entire flood 
hydrograph may be less than 24 hours, and yet an important flash flood peak 
may have occurred in that interval of time. The 24-hour averages may give 
an erroneous impression of the range of flood magnitudes with respect to fre- 
quency and lead to serious underdesign if hydraulic capacities are to be based 
on frequency criteria. 

Much of the research on flood frequency through the years has attempted 
to apply precise mathematical laws to a phenomenon which is basically em- 
pirical and to a varying degree, cyclical. Furthermore, flood magnitudes are 
affected by surface conditions in each particular basin, some of which are 
permanent and some of which seasonal (cyclical) and by storage effects which 
may vary widely with the degree of basin inundation. Only from a long record 
can be obtain a reliable discharge-frequency distribution. The true flood pro- 
ducing potentialities of a basin for use in immediate design studies can only 
be derived from more complex storm studies translated into flood flows by 
the unit hydrograph method. 
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